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After reviewing the ideal Bose-Einstein gas in a box and in a harmonic trap, I discuss the effect 
of interactions on the formation of a Bose-Einstein condensate (BEC), along with the dynamics 
of small-amplitude perturbations (the Bogoliubov equations). When the condensate rotates with 
angular velocity Q, one or several vortices nucleate, with many observable consequences. With 
more rapid rotation, the vortices form a dense triangular array, and the collective behavior of 
these vortices has additional experimental implications. For Q near the radial trap frequency 
Lj±, the lowest-Landau-level approximation becomes applicable, providing a simple picture of 
such rapidly rotating condensates. Eventually, as f2 — > u>j_, the rotating dilute gas is expected to 
undergo a quantum phase transition from a superfluid to various highly correlated (nonsupcrfluid) 
states analogous to those familiar from the fractional quantum Hall effect for electrons in a strong 
perpendicular magnetic field. 
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I. INTRODUCTION 

The remarkable creation of Bose-Einstein con- 
densates (BECs) in cold dilute alkalai-metal gases 
(lAnderson et aLl . Tl995l iBradlev ei~ai\ , Il995t iDavis et all , 
Il995l ) has produced a wholly new exciting field that con- 
tinues to thrive. Typically, these systems have a macro- 
scopic condensate wave function ty(r,t) that character- 
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izes the static and dynamic behavior of the BEC [for 
a gen e ral background , see, for example, Dalfovo et all 
(ll999l):unguscio et all Jl999h:lPethick and Smith! (|2002t): 
iPitaevskii and Stringaril ( 20031 )]. In the usual limit of 
a dilute gas, H> obeys a nonlinear Schrodi nger equation 
known as the G ross-Pitaevskii equation ( Grossl 1 196lt 
|Pitaevskiil . [l96l . 

Subsequent experiments formed rotati ng condensates , 
initial l y with a small number of vortices ( Madison et all 



l2000at iMatthews eToH. ll999aT). a nd th e n with large 



vortex arrays (|Abo -Shaeer et all . 1200 it lHalian et all 
l2001bl). As in the case of rotating superfluid 
4 He (jPonnellvl . Il99ll ) , the vortices have quantized circu- 
lation, arising from the single-valued nature of the con- 
densate wave function 'J'. Nevertheless, the new dilute- 
gas BECs differ considerably from superfluid 4 He because 
the Gross-Pitaevskii equation here provides a remark- 
ably detailed description of the physics, allowing a careful 
comparison of theory and experiment. 

This article focuses on the physics of quantized vortices 
in ultracold dilute trapped quantum gases. I first review 
the dilute Bose-Einstein gas, both in a box where the 
density is uniform and in harmonic traps where the den- 
sity is nonuniform (Sec. II). The presence of (repulsive) 
interactions have a dramatic effect, leading to (1) a col- 
lective mode with a linear long-wavelength dispersion re- 
lation and (2) superfluidity of macroscopic flow. Section 
III treats the behavior of one vortex (or a small number 
of vortices) in a rotating condensate. For more rapid ro- 
tations, the vortices form a regular array that has many 
analogies with rotating 4 He and with the flux-line lattice 
in type-II superconductors (Sec. IV). As the rotation rate 
f2 continues to increase and approaches the radial trap- 
ping frequency uj± , the effective radial trapping potential 
weakens, and condensate expands. In this limit, the sys- 
tem enters a new two-dimensional superfluid regime that 
has close analogies with the lowest-Landau-level descrip- 
tion of an electron in a uniform magnetic field (Sec. V). 
Ultimately, in the limit CI/oj± > 0.999, the rotating Bose 
gas is predicted to make a quantum phase transition to 
one of several highly correlated ground states that are 
not superfluid (Sec. VI). This regime involves a sequence 
of many-body states that are similar to those developed 
for the quantum Hall effect (a two-dimensional electron 
gas in a strong magnetic field). 

Thr oughout this article, I follow the lead of lBloch et all 
(|2007l ). focusing on models and predictions with exper- 
imental basis or confirmation. As a result, I treat only 
briefly two areas that have been studied in some exper- 
imental detail as nonrotating systems: (I) dipole con- 
densates that interact with long-range electric or mag- 
netic di pole forces [s e e, for example, iGriesmaier et all 
Lahave et all (|2007t ) and re ferences ther e in] (2 ) 



netic aipc 
(|2005l) : E 



spinor conden s ates [ se e, for example, Sadler et all ( 2006h ; 
IStenger et all (|l998t) : IVengalattore et all (|2007f) and ref- 
erences therein]. In both cases, however, rotation and the 
corresponding vortex structures remain an unexplored 
experimental capability (see, however, Sec. VII. B for an 



introductory discussion) . 



II. PHYSICS OF BOSE-EINSTEIN CONDENSATES IN 
DILUTE TRAPPED GASES 

The general topic of Bose-Einstein condensates (BEC) 
in dilute trapped gases has been the f ocus of severa l 
reviews and books (Dalfovo et all, 19991: Inguscio et al. . 



19991: iPethick and Smith! . l2002t IPitaevskii and Stringari 
The current article will emphasize those as- 
pects most relevant to rotating tra pped gas es and 
the associated quan t ized v ortices (lAftalionl. 2006; 
Fetter and Svidzinskvl. 1200 it iKasamatsu and Tsubota 
' l2007l)~ 



20071 : Parker et al. . 2007h . In particular. I will emphasize 



the regime of regular vortex arrays, both the "mean- field 
Thomas- Fermi" limit when the vortex cores remain small 
and the "lowest-Landau-level" limit w hen the cores are 
comp arable to the intervortex spacing ( Bavml . [20051 : iHol . 
120011) . 

To review the basic facts of Bose-Einstein condensa- 
tion, let us consider an ideal (noninteracting) uniform 
three-dimensional gas with N particles in a cubical box 
(volume V = L 3 ) with mean density n = N/V. If the 
gas has a temperature T, the mean momentum per parti- 
cle is pt ~ y/MhsT, ignoring numerical factors of order 
unity, where M is the atomic mass and ks is Boltzmann's 
constant. The de Broglie relation then yields the ther- 
mal wavelength At ~ h/px ~ hj \JMksT ■ The other 
relevant length is the interparticle spacing n" 1 / 3 . The 
short- wavelength limit At <C n -1 / 3 holds when H — ► or 
T is large; it describes the classical regime when quantum 
interference and diffraction are negligible (the analog of 
ray optics). As T falls, however, this inequality eventu- 
ally fails when At ~ n^ 1 ^, or, equivalently, when the 
"phase-space density" defined as nXj, becomes of order 
unity. This condition signals the onset of quantum de- 
generacy; whenever nA^ > 1, quantum mechanics plays 
a central role in the physics. 

The preceding description relies on an atomic-physics 
perspective, but a condensed-matter view yields a sim- 
ilar picture. In an ideal gas, each atom is effectively 
confined to a volume n~ l by its neighbors, leading to a 



zero-point energy e 



zp 



2 / 3 /M. At high T or low 
density (fcsT ^ e zp ), the system behaves classically, but 
when ksT sa e zp (either because of reduced temperature 
or increased density), quantum effects again become cru- 
cial. These two characterizations of the onset of quantum 
degeneracy are essentially equivalent. 

For ideal bosons, the critical temperature T c for the 
onset of quantum degeneracy leads to Bose-Einstein con- 
densation, whereas for ideal Fermi gases, the similar on- 
set temperature is known as the Fermi temperature Tp. 



For liquid helium with 



10 cm , the appropri- 



ate transition temperature is of order 1 K for both iso- 
topes ( 4 He is a boson and 3 He is a fermion). A typical 
bosonic dilute alkalai-metal gas (such as 7 Li, 23 Na, and 
87 Rb) has a much lower T c of order 100-1000 nK, because 
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of the larger atomic mass and reduced number density 
(ri ps 10 13 cm -3 is a typical value for dilute cold gases). 
This article will focus on such bosonic systems. 

The presence of Bose-Einstein condensation means 
that a single quantum state has macroscopic occupation. 
This particular quantum state acts like a particle reser- 
voir that can absorb or emit excited particles with negli- 
gible change in its own properties. Thus it is natural to 
use the grand canonical ensemble in which the system of 
interest is assumed to be in equilibrium with a reservoir 
at temperature T and chemical potential /i that deter- 
mine the mean total energy and mean total number of 
particles. 



A. Ideal Bose gas 

Consider an ideal Bose gas in an external trap potential 
Vtr, with a complete set of quantum- mechanical single- 
particle energies ej, and assume that the system is in 
equilibrium at temperature T and chemical potential /i. 
In the grand canonical ensemble, the mean occupation of 
the jth state is 



1 



exp [0 (ej - fi)] - 1 



(2.1) 



where [3 = 
temperature. 



(k B TY 
Here, 



is proportional to the inverse 



/(e) = {exp[/3(e-M)] - 1}- 



is 



the usual Bose-Einstein distribution function (it de- 
pends explicitly on T and At), as s hown , for exam- 
ple, in_Sec^_5_c£l^tterandmi^y (|2003l ). Chap. V 
of iLifshitz and Pitaevskiil (|l980al)~ Part 1, or Chap. 7 
of iPathrial (|1996h . A detailed analysis shows that the 
mean total number of particles N(T, /i) and the mean 
total energy E(T, /i) are given by 
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(2.2) 
(2.3) 



In principle, the first relation Eq. (|2.2[) can be inverted to 
give n(T,N), and substitution into Eq. (|2.3[) then yields 
the mean total energy E(T, N) as a function of temper- 
ature and the total number of particles. 

It is useful to introduce the "density of states" 



.9 



so that Eqs. (f2~2")l and (|2~3")1 reduce to 

N(T,fi) = f deg(e)f(e), 



E(T,fj) 



deg{e)ef(e). 



(2.4) 



(2.5) 
(2.6) 



Although g(e) is formally singular, a smoothed (coarse- 
grained) version would be well-defined. 



In the classical limit, the chemical potential /i c i = 
fcsTln(?T,A 3 r ) is large and negative since nX^ < 1. As the 
temperature decreases (or the density increases), how- 
ever, n increases toward positive values. The critical 
temperature T c for the onset of BEC is defined implic- 
itly through the relation /j,(T c , N) = eo, where eo is the 
ground-state energy of the single-particle Hamiltonian 
with the potential Vtr- For T < T c , the chemical po- 
tential remains fixed at this value, and the ground state 
has a macroscopic occupation Nq(T), whose temperature 
dependence follows from the conservation of total parti- 
cles N = N (T) + N'(T), with 



N'(T) = de — , 
K ' Je exp[/3(e 



9(e) 



eo)] ~ 1 



(2.7) 



defining the total number of particles not in the conden- 
sate. If <?(eo) vanishes, then this integral is finite and 
well defined, and the number of noncondensed particles 
decreases with decreasing temperature. In contrast, if 
(/(eo) is finite or singular, then this integral for the num- 
ber of noncondensed particles diverges, and the system 
cannot form a BEC. It is helpful to examine two specific 
examples. 



1. Ideal BEC in a box with periodic boundary conditions 

The textbook example of BEC is an ideal gas in a 
three-dimensional cubical box of volume V = L 3 with 
periodic boundary conditions. The single-particle eigen- 
functions are plane waves "0k ( r ) = V r_1 / 2 exp(ifc • r) 
with single-particle energy efc = h 2 k 2 / (2M), where fe = 
(2ir/ L)(n x , n y , n z ) with rij any integer. It is not hard to 
obtain the corresponding density of states 



9(e) 



V 



4tt 2 \ h 2 



(2.8) 



In a box with periodic boundary conditions, the lowest 
single-particle energy is eo = 0, which corresponds to 
k = (a uniform state with density l^ol 2 = V" 1 ). A 
standard analysis with the conventional definition Ay = 
2-Kh 2 /(Mk B T) yields the condition 



= C(1) 



2.612, 



(2.9) 



which is the critical phase-space density for the onset of 
BEC. Alternatively, the critical temperature T c is given 
by 



2tt 



h 2 n 2 / 3 



[C(|)] 2/3 



M 



3.31 ■ 



H 2 n 2 / 3 



M 



(2.10) 



as anticipated from the previous qualitative discussion. 
Below T c , the (uniform) condensate with zero momen- 
tum has a macroscopic occupation number, which is the 
signature of BEC. In particular, it is easy to verify that 



N'(T) 
N 



3/2 



(2.11) 
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for T < T c , so that the macroscopic occupation of the 
uniform ground state has the temperature dependence 



N (T) 



T 

V ' " ( T r 



3/2 



(2.12) 



It is worth generalizing these results to a uniform ideal 
Bose gas in a d-dimensional hypercubical box; the corre- 
sponding density of states g(e) is proportional to e d / 2 ~ x . 
This alteration means that the integral in Eq. fl2.7|) di- 
verges for d < 2. As a result, the limit /i —> eq requires 
a more careful analysis, leading to the conclusion that 
such a uniform gas in two or one dimension cannot form a 
BEC; for a detailed treatment of the two-dimensional un i- 
form Bose gas, see Chap. 2 of iPethick and Smith! (|2002fh 



2. Ideal BEC in a harmonic trap 

In three dimensions, an anisotropic harmonic trap po- 
tential 

V tI (r) = \M (oj 2 x x 2 + uy + ulz 2 ) (2.13) 

yields single-particle energies that are labeled by a triplet 
of non- negative integers n x ,n yi n z 



h{n x L0 x + n y Ljy + n z uj z ) + e , (2.14) 



where the second term is the zero-point energy eo — 
ifi. [u> x + ojy + u) z ). The ground-state wave function is 
a product of three one-dimensional Gaussians with mean 
widths dj = \Jhj (Muij) (here j — x, y, or z). If the 
sums in the density of states are approximated by inte- 
grals (which holds for large e) , the corresponding density 
of states is 



9(e) 



2h\ 



(2.15) 



where defines a geometric-mean trap fre- 

quency. 

The onset of BEC in a harmonic trap occurs at /i = eo, 
and an elementary analysis yields the transition temper- 
ature 

k B T c = [C(3)P 1/3 huoN 1 / 3 sa OMhujoN 1 ^ 3 , (2.16) 

where £(3) w 1.212. Below T c , a macroscopic number of 
particles occupies the anisotropic Gaussian ground state 
of this harmonic trap, which has a macroscopic occupa- 
tion with temperature dependence 



No(T) (T 



(2.17) 



Images of this temperature-dependent Gaussian conden- 
sate density profile provided th e first clear evid e nce fo r 
the creation of BEC in 87 Rb (lAnderson et all , 1 19951 ). 
Typical traps have do — y/ti/(Mu)o) ~ a few (an and 



N ~ 10 6 , confirming the previous value T c « 100-1000 
nK (depending on the atomic mass). 

For a harmonic trap in d dimensions, the density of 
states g(e) is proportional to e d_1 , which differs from 
the d dependence for a hypercubical box. As a re- 
sult, a two-dimensional Bose gas in a trap can form a 
BEC with a finite transition temperature T c , in con- 
trast to a uniform two-dimensional gas in a box. This 
behavior is especially important in the limit of rapidly 
rotating Bose-Einstein condensate, when the centrifu- 
gal forces flatten the atomic system, producing an ef- 
fectively two-dimensiona l trapped gas. Chapter 2 of 
IPethick and Smith! (|2002l) discusses this situation in some 
detail, especially the behavior in a one-dimensional trap, 
where the order of limits N — ► oo and T — > becomes 
crucial. 



B. Inclusion of interparticle interactions 



Roughly sixty years ago. lBogohubovj ( 1947?) introduced 
what is now seen as the essential physical approxima- 
tion for a dilute Bose gas with repulsive interactions. To 
understand his idea, note that the repulsive interparti- 
cle interaction can be characterized by the positive s- 
wave scattering length a, which is typically a few nm for 
th e dilute alkalai-metal ato ms of interest [see Chap. 5 
of IPethick and Smith! ( 2002)]. In the low-density limit 



(a -C n -1 / 3 ), the "gas parameter" na 3 is small and pro- 
vides a natural expansion parameter. Bogoliubov notes 
that the ground state of the dilute gas will then be close 
to that of the ideal gas. Thus the number of particles 
Nq in the condensate will be close to the total number 
N, and the difference can be ignored in first approxima- 
tion. The macroscopic occupation of the single-particle 
mode for the condensate means that the creation and 
annihilation operators for this particular mode can be 
approximated as "classical fields" similar to the electric 
fields for a laser mode. Bogoliubov's seminal physical pic- 
ture underlies almost all the subsequent developments in 
the field and is basic for understanding the behavior of a 
low-temperature BEC. 

For a uniform Bose gas at rest, the condensate is 
the zero-momentum state, and the excited states are 
labeled by wavenum ber fc ^ 0, as discussed, for ex- 
ample, in Sec. 25 o f Lifshit zTand Pitaevskiil ( 1980bl ). in 
Sec. 7 2 of IPethick and Smith! (l2002t) . and Sees. 18 and 
35 of [Fetter and Waleckal (|2003t h In the present case 
of a trapped condensate, however, the analysis is some- 
what more intricate and relies on the macroscopic con- 
densate wave function ^ (sometimes called the "order 
parameter"). This condensate wave function describes 
the single-particle mode that has macroscopic occupa- 
tion, with its squared absolute value ^(r)! 2 = n(r) giv- 
ing the nonuniform condensate particle density n(r). For 
a dilute gas at low temperature with Nq f=a N, the nor- 
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malization requires 

J dV\^{r)\ 2 = Nq ~ N. 



(2.18) 



In the limit ol a nearly ideal Bose gas at T = K, the 
spatial form of this condensate wave function follows by 
minimizing the total energy E, along with the constraint 
that the total number of particles N is conserved. Equiv- 
alently, one can simply minimize the grand-canonical 
thermodynamic potential E — fj,N, where fi is the chem- 
ical potential. To understand this approach, I shall gen- 
erally rely on physical arguments, but many-body quan- 
tum field the ory provides a rigorous foundation for the 
same results ( Fetterl . Il972l Il996t iHohenberg and Martini . 
Il965f l. In a trap with confining potential V tr, the Gross- 



Pitaevskii (|Grossl . Il96lt iPitaevskiil . Il96lh energy 
tional for the condensate has the form 



func- 



GP 




kinetic 



interaction/ 

(2.19) 

containing the kinetic energy, the trap energy [see 
Eq. (|2.13p ]. and the interaction energy, respectively. The 
first two terms (those quadratic in ^) are just the one- 
body energy for an ideal Bose gas in a (usually) har- 
monic trap. In contrast, the quartic term (the two-body 
energy) describes the effect of interactions, where the in- 
terparticle potential has been approximated by a short- 
range "pseudopotential" V(r — r') ~ gS^'ir — r'), with 
g = 4irah 2 /M a coupling constant fixed by the s-wave 
scattering length a [this result follows from standard two- 
body quantu m-mechanical scattering th eory, for exam- 
ple, Chap. 5 of iPethick and Smith! {2002) or Sees. 11 and 
35 of lFetter and Waleckal (|2003l )]. 

Comparison of the kinetic and trap energies for a har- 
monic potential yields the familiar oscillator length 



do = 



Mu 



(2.20) 



that characterizes the mean size of the noninteracting 
condensate, with luo = (uj x u) y uj z ) 1 / 3 . Similarly, compari- 
son of the kinetic energy and the interaction energy for 
a uniform condensate yields the "healing length" 



1 



y/2Mgn y/Siran' 



(2.21) 



which characterizes the length scale over which a local- 
ized alteration in the condensate density heals back to its 
uniform value n. In this context, recall the Bogoliubov 
approximation of small quantum depletion N' <C N, 
which requires na 3 <C 1. As a result, Eq. (|2.21[) has 
the important corollary 



1 



87rn 1 / 3 a 



» 1, 



(2.22) 



namely the healing length must exceed the interparticle 
spacing in any GP description. 

Variation of E GP [$i] in Eq. (|2.19|) with respect to ** 
at fixed normalization yiel ds the time- independent Gross - 
Pitaevskii (GP) equation (jGrossl . 1 19611: IPitaevskiil ll96lD 



h 2 v 2 

' 2M 



+ V tr (r)+g\^(r)f 



W(r) = fx^(r), (2.23) 



Here the chemical potential \i can either be consid- 
ered a Lagrange multiplier or a parameter in the zero- 
temperature grand-canonical thermodynamic potential 
E — fiN. In the present context, this GP equation is es- 
sentially a nonlinear Schrodinger equation that includes 
an additional quadratic self coupling g| , 3 / | 2 ; such a term 
can be interpreted as a Hartree potential Vh(t) = gn(r) 
that represents the interaction with the local nonuniform 
condensate density. 

Compared to a uniform BEC, a trapped condensate 
involves an additional characteristic length do (the non- 
interacting co ndensate size) . A simple scaling argument 
[see Sec. 6.2 of lPethick and Smithl ( 20021 )] yields a new di- 
mensionless parameter Na/do that characterizes the im- 
portance of the interaction in a trapped condensate. As 
a result, the three terms in Eq. (|2. 19[) can have very dif- 
ferent magnitudes depending on the parameters. In the 
usual situation (N m 10 6 , a w a few nm, and do ~ a 
few /im), this dimensionless parameter Na/do is large, 
and th e resulting regime is known as the "Thomas- Fermi 
limit" (|Bavm and Pethickl [1996 s ). 1 In this case, the re- 
pulsive interactions dominate and expand the condensate 
to a mean radius Ro that greatly exceeds the mean oscil- 
lator length do (a factor of 10 is typical). This expansion 
dramatically reduces the radial gradient of the density, 
and the associated kinetic energy thus becomes negligi- 
ble relative to the trap energy and the interaction energy. 
Equation (|2.19[) then reduces to 



E TF [V] » J ^(y tr |*| 2 + ig|*| 



(2.24) 



which involves only | , I'| 2 and | , I / | 4 . Minimization with 
respect to l^l 2 at fixed normalization gives the Thomas- 
Fcrmi (TF) approximation 



V tr (r)+g\^{r)\ 2 = ^ 



(2.25) 



which also follows by omitting the kinetic energy term in 
Eq. (|2.23p . This algebraic equation can be solved directly 
for the equilibrium density 

n ( r ) = - VtI (r)] , (2.26) 



1 As seen below, the limit Na/do 2> 1 means that the condensate 
density is effectively uniform on relevant length scales, which 
explains the name "Thomas-Fermi." 
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where 6{x) is the unit positive step function. Although 
this expression holds for any reasonable trap potential, 
the usual case of a quadratic harmonic trap (|2.13|) yields 
the simple and explicit expression 



n(r) = n(0) 1 



x 



z 

m 



(2.27) 



where the right-hand side is positive and zero other- 
wise. Here n(0) = fi/g is the central density and R 2 = 
2/x/(Mw|) are the squared condensate radii in the three 
coordinate directions. This TF density has a parabolic 
cross section and fills the interior of an ellipsoid. 

Application of the normalization condition (|2.18p to 
the TF density (T2T27)) yields N = 87rri(0)i^/15, where 
Rq = R x R y R z depends on the chemical potential \i. 
Some algebra shows that n(0) = fi/g — R 2 , / (Snad^) , and 
a combination of these results yields the important di- 
mensionless relation 



Na 



r 5 



(2.28) 



which is large in the present TF limit. Correspondingly, 
the TF chemical potential becomes 



fiTF 



d 2 ' 



(2.29) 



so that /itf 3> fiujo in the TF limit. Since /xtf is propor- 
tional to TV 2 / 5 , the thermodynamic relation fi = dE/dN 
yields the TF energy for a trapped condensate 



f MTFiV, 



(2.30) 



which also follows by direct integration of Eq. (|2.24[) for 
the TF density profile. 

It is conventional to use the central density n(0) to 
define the healing length [Eq. (|2.21[) ] in a nonuniform 
trapped condensate, which gives 



87r<m(0) 



(2.31) 



These relations lead to the important conclusion that 
£ 2 = <Iq/Rq; equivalently, 



1 

do 



dp 

Ro 



(2.32) 



This TF limit holds when the mean condensate radius 
Rq is large compared to the mean oscillator length dp. 
Thus the TF oscillator length do is the geometric mean 
of £ and Rq, and Eq. (|2.32p then yields a clear separation 
of TF length scales, with £ <C d <C R . 

Equation (|2.23p has the expected time-dependent gen- 
eralization (known as the time-dependent GP equation) 



ih 



dt 



ft 2 V 2 
" 2M 



+ V tI (r)+g\y(r 7 t)\' 



*(r,t). 
(2.33) 



where * now depends on t as well as on r. Al- 
though I shall not discuss the field-theoretical deriva- 
tion of the time-dependent GP equation, it is worth 
noting that the condensate wave function is interpreted 
as an ensemble average ^(r,t) = (ijj(r,t)) of the 
field operator ijj(r,t) p 
tinf 



(IFetter 



iHohenberg and Martini 1965 



19721 [19961: iGrosd 
Pitaevskiil ["l96lh . 



Com- 
parison of Eqs. (|2.23p and (|2.33p implies that a station- 
ary solution has the time dependence exp(— i/j,t/h). This 
result also follows directly from the field-theoretical def- 
inition because the states on the left-hand side of the 
ensemble average have one less particle than those on 
the right-hand side. 

This nonlinear field equation (|2.33p can be recast in an 
intuitive hydrodynamic form by writing the condensate 
wave function 



tf(r,t) = |tt(r,i)|exp[iS(r,t)] 



(2.34) 



in terms of its magnitude \^\ and phase S. In this 
picture, the condensate (particle) density is n(r, t) = 
|\P(r,t)| 2 , whereas the particle-current density becomes 



2Mi 



1*1 



M 



(2.35) 



The usual hydrodynamic relation j = nv identifies the 
local velocity as 



v{r,t) = —VS(r,t) = V$(r,t), 



(2.36) 



w here <i> = HS/M is the velocity potential, as proposed 
bv lFevnmanl ( 1955h in his discussion of the quantum me- 
chanics of vortices in superfluid He II. Note that V X v 
vanishes wherever S is not singular, so that the "super- 
fluid v elocity" v is irrotational, as assumed by lLandaul 
(1941) in his phenomenological two-fluid hydrodynamics 
for superfluid He II. For more complete discussions, see , 
for example, Chap . Ill of Lifshitz and Pitaevsk ii (1980b) 
and Chap. XVI of lLandau and Lifshitzl (|1987[ ) 

The representation of the velocity v in terms of the 
gradient of the phase has one major implication. Con- 
sider the "circulation" k, defined at a given instant of 
time in terms of a line integral n ~ <jL dl ■ v around a 
closed path C. Substitution of Eq. (|2.36p yields 



K =±t d l.VS=^-AS\e, 



M 



M 



(2.37) 



where AS'lc is the net change in the phase on following 
C. The single- valuedness of the condensate wave func- 
tion requires that this quantity must be an integral mul- 
tiple of 27r, which shows that the circulation in a dilute 
BEC must be q uanti zed in units of 2irh/M (Fevnman, 
119551: lOnsaeerL Il949h . For accounts of this remark- 
able q uantum-mechan ical relation, see Sec. 23 of lLondonl 
(119541) and Sec. 2.3 of iDonneilvl (|l99lh . 

Substitute Eq. (|2.34p into the time-dependent GP 
equation (|2.33p . The imaginary part yields the familiar 
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continuity equation for a compressible fluid 

+ V • (nv) = 0. (2.38) 

at 

Correspondingly, the real part provides a generalized 
Bernoulli equation for this quantum fluid 



1 



-Mv 2 + V tl 



H 2 n <9<f> 

M 7 ^ 2 V* + 9 n + M- = 0, (2.39) 



where the explicitly quantum-mechanical term contain- 
ing h 2 is sometimes called a "quantum pressure." The 
full structure of this Bernoulli equation (including the 
quantum pressure) can be understood as that for an 
irrotational compressible i sentropic fluid [see Sec. 2 of 
lLandau and Lifshitd (If 987 )1 w i th an enthalpy density ob- 
tained from Eq. ([239]) (jFettel |2002|) . Such an isentropic 
picture is clearly appropriate for a low-temperature su- 
pcrfluid. 

The existence of a Bernoulli equation for the time- 
dependent GP equation has one basic consequence in 
the present context of vortices in dilute quantum gases. 
All the classical results about vortex dynamics in an ir- 
rotational fluid [such as Kelvin's theorem on the con- 
ser vation of circulation , prove d, for example, in Sec. 8 
of lLandau and Lifshitzl ( 19871) ] automatically hold for 
any vortex configuration in a dilute BEC. Nevertheless, 
the presence of a nonuniform trapping potential Vt T (r) 
and the resulting nonuniform density in Eq. (|2.26|) sig- 
nificantly affect the dynamics of vortices, so that quali- 
tative classical pictures based on uniform incompressible 
fluids may not always apply. 



C. Bogoliubov equations 



The time-dependent GP equation (|2.33[) describes the 
dynamics of the condensate, and it is natural to consider 
the linearized behavior of small perturbations around 
such solutions ^(r,t). Alth ough a fu l ly quantum- 
mechanical analysis is feasible ( Bogoliubov], Il947t iFetterl . 
If972t IPitaevskiil. Il96ll). I here rely on a simpler classi- 



cal approach (iDalfoyo et aH , If999 ) [see also Sec. 7.2 of 
iPethick and Smith! (|200^) ]~ that yields the same eigen- 
value equations. Assume a stationary nonuniform con- 
densate with wave function ^>(r) e ~ l > lt / h . Since the 
nonlinear Hartree term involves | V I'| 2 , the perturbations 
must, for consistency, include e Tlu>t with both signs of 
the frequency. Specifically, take 



#(r,t) = eT i ^l % [tf(r) + u(r). 



v*{r)e 1 ' 



(2.40) 

where the Bogoliubov amplitudes u{r) and v(r) are 
treated as small. 

Substitute (|2~40| into Eq. |2~33| and collect first-order 
terms proportional separately to e TluJt . The result is a 
pair of linear eigenvalue equations 



Cuj — g{^!) 2 Vj — HujjUj, 
Cvj-g(**fuj = 



(2.41) 
(2.42) 



where £ = -ft 2 V 2 /(2Af) + V tr -^ + 2 i g|*| 2 is a Hermitian 
operator. These coupled equa tions are known a s the Bo- 
goliubov equations, although IPitaevskiil ( 196lh was the 
first to consider the application to nonuniform systems 
[sometimes they are called the Bogoliubov-de Gcnncs 
equations because of their close correspondence to similar 
equations in the theory of supe rconductivity, as discussed 
i n Cha p. 5 of lde Gennesl (|l966l ) and Sec. 10.1 of lTinkhaml 
(119961) ]. In many ways, the Bogoliubov equations are 
analogous to a nonrelativistic version of the Dirac equa- 
tion, with u and v as the particle and hole amplitudes, 
including the (+, — ) metric seen in the minus sign on the 
right side of Eq. (pT42|) compa red to t hat o f Eq. (f2~4Tj) 
[see, for exa mple, Chap. 12 o f lAbersI (20041) Sec. 52 of 
ISchifj (|1968| ). and Chap. 20 of IShankarl (|l994h ]. 

For a uniform condensate, the appropriate eigenfunc- 
tions are plane waves labeled by k. The Bogoliubov en- 
ergy eigenvalue is easily determined to be 



Eh — 



gnh 2 k 2 
M 



h 2 k 2 

~2M 



-l 1/2 



(2.43) 



where n is the condensate density and /i ~ gn for this 
uniform gas. For long wavelengths (kt; C 1), where £ is 
the healing length from Eq. (|2.2ip . the Bogoliubov energy 
reduces to a phonon spectrum 



E k w hsk 



(2.44) 



with s = gn/M the speed of compressional sound. For 
short wavelengths (fc£ 3> 1), in contrast, the Bogoliubov 
spectrum becomes that of a free particle. The existence 
of a linear (phonon) spectrum at l ong wav e length s is cru- 
cial for superfluidity because the ILandaul ( 194lh critical 
velocity v c for onset of dissipation is the minimum value 
of Ek/(hk) with respect to k. In the present case, the Bo- 
goliubov excitation spectrum (|2~43"1) yields v c = s. Note 
that s cx yfg is real and positive only because of the repul- 
sive interactions (g > 0). Furthermore, v c = s vanishes 
for an ideal Bose gas with 5 = 0. Hence the ideal Bose 
gas is not strictly superfluid even at zero temperature, 
although it does indeed have a condensate with k = 0. 
In this well-defined sense, the existence of superfluidity 
in a uniform dilute Bose gas arises from the repulsive 
interactions. 

Manipulation of the Bogoliubov equations shows that 



■JdV(\t 



is real. In most cases, the integral 



is nonzero, so that u)j itself is real. To understand the ba- 
sic issue involved in the normalization of the Bogoliubov 
eigenf unctions, it helps briefly to recall the usual parti- 



from 



cle creation and annihilation operators and aj 
quantum field theory. These particle and hole operators 
obey the familiar Bose-Einstein commutation relations 
[aj-,at,] = Sjj'. The solution of the Bogoliubov equa- 
tions constitutes a canonical transformation to a new set 
of "quasiparticle" operators and aj that are linear 
combinations of the original particle and hole operators. 
The situation is very similar to that for BCS supercon- 
ductors, as shown, for example, in Sec. 5-1 of lde Gennesl 
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(|l966l ) or Sec. 3.5 of lTinkhaml (| 19961) . Almost by defi- 
nition, any canonical transformation preserves the com- 
mutation relations, so that the Bogoliubov quasiparticle 
operators must obey the same Bosc-Einstein commuta- 
tion relations [ay, at,] = Sj t j>. A detailed analysis shows 
that physically relevant Bogoliubov eigenfunctions must 
satisfy the following positive normalization condition 



dV{\i 



1. 



(2.45) 



Under this canonical transformation, the Bogoliobov 
approximate grand-canonical Hamiltonian assumes the 
very simple form 



For some special condensed systems (typically when 
BEC occurs in an excited single-particle state), the lin- 
earized Bogoliubov equations also have zero-norm solu- 
tions with J dV (\u\ 2 — \v\ 2 ) = 0. In these cases, the fre- 
quencies can be (and usually are) complex. Such systems 
are said to exhibit a "dynamical instability," and the un- 
stable mode(s) will grow exponentially until various non- 
linearities limit the evolution. Such dynamically unstable 
behavi or occurs in a condensate with a d oubly quantized 
vortex ( Pu et all 11999^ Shin et a71 . l2004f ) in certain one- 
dimensional solitons jlFeder et all \200& iFedichev et all 



tical lattices (|Wu and N iu. 2001) 



1999; lMurvshe v et aZ.l.ll999j). and in one-dimensional op- 



H — fiN 1=3 const. 



(2.46) 



where the sum runs over all states with positive normal- 
ization. 2 It is clear that the sign of the eigenfrequencies is 
crucial for the stability of the system. If all the frequen- 
cies u>j are positive, then creation of quasiparticles raises 
the energy, and the system is manifestly stable. If, how- 
ever, any of the positive-norm solutions has a negative 
frequency, then creation of those particular quasiparti- 
cles lowers the energy, and the system is said to display 
an "energy instability." 

A simple example of this situation is a uniform con- 
densate that moves with velocity v = Hq/M. The con- 
densate wave function now has a nonuniform complex 
structure \P(t") = y/ne lq r , and the solution with positive 
normalization has the excitat ion energy (jFetterl . 1 19721: 
iPollockL fl967t iRokhsarl fl997cJ) 



E k (v) = hk ■ v + E k , 



(2.47) 



where Ek is the Bogoliubov energy (|2.43[) for a stationary 
condensate. For small v, this excitation energy is positive 
for any k. When v exceeds the speed of sound s, however, 
then this eigenvalue will be negative in certain directions, 
which represents the onset of the Landau insta bility as- 
socia ted with spontaneous creation of phonons (Landau], 
Il94ll ). Although the energy instability implies that the 
system is formally unstable, the absence of dissipation 
means that no mechanism exists to support the growth 
of the unstable process. 



III. PHYSICS OF ONE VORTEX (A FEW VORTICES) IN 
A TRAP 

When the GP equation was developed in 1961 to pro- 
vide a tractable model for a vortex line in a BEC, rotat- 
ing superfluid 4 He was the only system of experimental 
interest (even though 4 He is a dense highly correlated 
system). Since 1995, however, dilute ultracold trapped 
alkalai-metal gases have provided remarkable physical 
systems that indeed fit the GP picture in considerable 
detail. 



A. One vortex in unbounded condensate 



iGrossI (|!96lh and lPitaevskil (|l9fill ) introduced the GP 
equation specifically to describe a single straight vortex 
in an otherwise uniform dilute Bose-Einstein condensate 
with bulk particle density n. The condensate wave func- 
tion has the form W(r) = \/nx( r )i with 



x (r) = e**/(r/0 



(3.1) 



where r, <f> are two-dimensional plane-polar coordinates 
and / — > 1 for large r. Equation (|2.36p gives the circu- 
lating velocity of a singly quantized long straight vortex 
line 



v(r) = — 0; 



(3.2) 



2 The criterion for inclusion in this sum has been a source of con- 
fusion, because every acceptable solution of the Bogoiubov equa- 
tion {uj,Vj} with frequency ujj and positive normalization has 
a second "antiparticle" solution {v*,u*} with frequency — luj 
and negative normalization. The original study of these rela- 
tions (Fetter, 1972) incorrectly assumed that the positive-norm 
solution always had positive energy. This result indeed describes 
the simplest case of a uniform stationary condensate with plane- 
wave eigenfunctions. In contrast, nonuniform condensates, such 
as one in uniform motion ( discussed here) or one containing a 
vortex [see Sec. III.B.2 and IRokhsarl jl997alH can have physical 
states with positive normali zation and negative frequency . For a 
more general discussion, see IFetter and Svidzinskvl J200lh . 



this hydrodynamic flow has circular streamlines and a 
magnitude that diverges as r — > 0. The circulation 
k = § dl ■ v = 2irh/M can be transformed with Stokes's 
theorem as J dS ■ V X v = 2irh/M, implying a singular 
localized vorticity at the center of the vortex core 



Vxv = — 6W(r)z. 



(3.3) 



The kinetic energy per unit length of the vortex is 



2M 



d 2 r |V*| 2 



h 2 n 
2M 



d 2 r 



, (3.4) 
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where the first term arises from the density vari- 
ation near the vortex core and the second term 
[= \M J d 2 rn\f(r)\ 2 v 2 ] is the kinetic energy of the cir- 
culating flow. With this kinetic- energy functional, the 
Euler-Lagrange equation yields the nonlinear GP equa- 
tion for the radial amplitude f(r). The resulting centrifu- 
gal barrier forces the amplitude f(r) to vanish linearly for 
r < £, which characterizes the vortex core, and f(r) ~ 1 
for r ^> £. In the present dilute limit, Eq. (|2.22|) shows 
that the vortex core is significantly larger than the inter- 
particle spacing. In addition, the particle-current density 
|j(r)| = n\f(r)\ 2 h/ (Mr) vanishes both far from the vor- 
tex and at the center of the vortex, with a maximum near 
r « £. Finally, the spe ed of s ound in a dilute Bose gas 
can be rewritten s = yfgn/M = %/{\/2M£), so that the 
circulating flow in Eq. (|3.2|) becomes supersonic for r < £; 
in this sense, local acoustic cavitation can be considered 
the source of the vortex core. 

In terms of the dimensionless scaled variable 
x = r/£, the variational approximation f vaI (x) = 
x/y/x 2 + 2 (jFetterl . fl969l ) provides a good fit to the nu- 
merical solution of the radi al GP equation [see , for e x 



of iLifshitz and Pitaevskiil (|l980af) ]. 3 



ample, Fig. 9.1 in Sec. 9.2 of iPethick and Smith! ||2002h ]. 
The additional GP energy per unit length of vortex is 
E v ps (nh 2 n/M) ln(1.46i?/£), where R is a large-distance 
cutoff. Apart from the additive numerical constant, this 
result is essentially the integral of ^Mv 2 n over a large 
ci rcle of radius R with a short-d istance cutoff of order 
£ ( Ginzburg and Pitaevskiil Il958l ). 



B. One vortex in a trapped condensate 

A classical viscous fluid in a container that rotates with 
angular velocity fi rapidly acquires the same angular ve- 
locity because of the microscopically rough walls. From 
a theoretical view, the laboratory is no longer the ap- 
propriate reference frame because the rotating walls are 
a time-dependent external potential that can do work 
on the system. Thus it is essential to transform to the 
frame rotating with the walls, in which case the exter- 
nal potential becomes stationary. Consider an originally 
static Hamiltonian H(r,p) with an external trap poten- 
tial Vt v (x,y, z) that in general is anisotropic about O. 
When the potential rotates with angular velocity 17, a 
direct transformation of the Lagrangian from the labo- 
ratory frame to the rotating frame eventually yields the 
Hamiltonian H'(r' ,p') in the rotating frame 



H'(r',p')=H(r',p')-n-L(r',p') 



(3.5) 



where r' and p' are the coordinates and canonical mo- 
menta in the rotating frame. Here, the right-hand side 
involves the nonrotating Hamiltonian H(r',p') and an- 
gular momentum L(r',p') = r' X p', with the origi- 
nal variables r eplaced by those in t he ro tating frame 
[see Sec. 39 of lLandau and Lifsh"it3 ( 1960h and Sec. 34 



1. Gross-Pitaevskii energy in a trap 

In the context of the GP equation, this transformation 
to a rotating frame yields a modified energy functional 



dV 



2M 



;g\v(r)\ z 



|V*(r)| 2 -(-y tr (r)|*(r)| 2 (3.6) 
-SI- J dVZ*(r)r x p*(r), 



where all variables here are in the rotating frame (writ- 
ten without primes for notational simplicity) [compare 
Eq. (|3.5p ]. The minimization of £"[\&] now explicitly in- 
volves the angular velocity, and the last term (— fi • L) 
clearly favors states with nonzero (positive) angular mo- 
mentum. 

Most experiments on vortices in rotating dilute BECs 
occur in the TF regime. Thus I shall concentrate on this 
situation, which has several simplifying features. Specifi- 
cally, the dominant repulsive interactions imply that the 
mean condensate radius Rq is large relative to the original 
mean oscillator length do. Equation (|2.32[) then implies 
that the vortex core size (~ £) is smaller than do by the 
same ratio do/Ro- Thus the presence of a few vortices 
does not significantly affect the number density, which 
can be taken as the unperturbed TF shape in Eq. I|2.27[) . 
In most cases, the condensates are axisymmetric with 
radial and axial trap frequencies ui±_ and oj z , so that the 
nonrotating TF density depends only on r = \/ x 2 + y 2 
and z. 

If an axisymmetric condensate contains a vortex 
at th e center, the genera l wav e function h a s the 
form (iDalfovo and Stringaril . Il996l ; iLundh et all Il997t 
ISinhatll997fi ~ 



#(r) = e i4> |*(r, z)\ 



(3.7) 



with the phase <fi appropriate for a singly quantized vor- 
tex. The hydrodynamic velocity is again v — h<f>/{Mr), 
as in Eq. (|3.2[) for a singly quantized vortex in an asymp- 
totically uniform fluid. The resulting centrifugal barrier 
near the axis of symmetry creates an axial node in the 
condensate wave function; thus the presence of one vor- 
tex fundamentally alters the topology of the condensate 
density from ellipsoidal to toroidal because of the hole 
along the vortex core. 

It is important to generalize this wave function to 
an off-center vortex. For simplicity, consider a three- 
dimensional disk-shaped condensate (u> z 3> u>±) with ra- 
dial condensate dimension R±. In this case, a vortex 



3 Note that H(r' , p') is not simply the original Hamiltonian in the 
laboratory frame because of the primed variables. 
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at ro = (xo,yo) will remain essentially straight, and the 
principal effect of the lateral translation is the altered 
phase 



S(x, y) = arctan 



y-yo 



x - 



(3.8) 



which is the azimuthal angle around the shifted origin ro. 
In the TF regime, the density is taken to be unchanged 
(the logarithmically divergent circulating kinetic energy 
is cut off at the vortex core radius £) . 

Use the energy functional E'[ty] in Eq. (|3.6[) and let 
E' denote energy of a vortex-free state (it is indepen- 
dent of O) . If E[ (ro , O) is the total GP energy of a 
single off-center vortex in the rotating frame, then the 
difference of the two energies A£"(ro,f2) = E[(ro,il) — 
Eq represents the formation energy of the vortex. A 
straig htforward analysis with this model TF wave func- 
tion (jSvidzinskv and Fetter! . l2000al) yields 



AE'(r ,n) 
A£'(0,0) 

where 



= 1- 



3/2 



2 MR 2 , 



hi 



5/2 



(3.10) 



and the detailed form of AE'(0, 0) is unimportant here. 
For fixed fi, this expression gives the energy in the 
rotating frame AE'^Qjfl) as the vortex moves from 
the center of the condensate (ro = 0) to the edge 
(r = R i ) , shown in Fig. [1] for various fixed values 
of f2 [Packard and Sanders! (jl972h make a similar anal- 
ysis for a vortex i n an in compressible uniform supcrfluid; 
ICastin and Duml (jl999) obtain comparable results nu- 
merically for a vortex in a trapped condensate] . 

(a) Curve (a) shows the energy AE'(ro,0) of a vortex 
in a nonrotating condensate as it moves from the 
center to the outer edge. This energy decreases 
monotonically as ro increases, with negative curva- 
ture for ro < i?j_/\/2. If there is no dissipation, 
the distance ro determines the energy, and the vor- 
tex follows a circular trajectory at fixed r . At 
low but finite temperature, however, weak dissipa- 
tion slowly reduces the energy of the vortex, which 
thus moves outward along curve (a), executing a 
spira l trajectory tow ard the edge of the conden- 
sate (|R.okhsail[l997al lbh. 



(b) With increasing fi, the function AE' (r ,il) for 
fixed fl flattens and the (negative) central curva- 
ture decreases. Curve (b) has zero central curva- 
ture, which occurs for 



2 MR\ 



In 



R± 



(3.11) 



This angular velocity O m represents the onset of 
metastability for a vortex near the center of a ro- 
tating condensate. For f2 < f2 m , the negative cen- 
tral curvature means that weak dissipation moves 



the vortex away from the center, and it eventually 
spirals out to the edge. For O > Q m , however, the 
central curvature is now positive, and weak dissi- 
pation will move the vortex back to the center. In 
this latter regime, the surrounding barrier makes 
the vortex locally stable near the trap center, even 
though it is not globally stable because the energy 
at the edge is definitely lower than the energy at 
the center. 



(c) Eventually, at fl c defined in Eq. (|3.10[) , the cen- 
tral energy A£"(0, f2 c ) vanishes [and thus equals 
E'(R±,n c )]- For O > fl c , the central vortex is 
both locally and globally stable. This value f2 c 
is interpreted as the thermodynamic critical angu- 
lar velocity for the creation of a singly quantized 
vortex in the TF limit. 4 As discussed below, ex- 
periments on rotating trapped gases generally find 
a different (larger) critical angular velocity for the 
creation of a vortex, but similar analyses for incom- 
pressible fluids in cylindrical containers do provide 
a reasonable descripti on of vortex formation in ro- 
tating superfluid 4 He (iPackard and Sandersl Il972t 



lYarmchuk and Packard , 
5.1 and 5.2 of iDonnellv 



1991ft ] 



1982) [see also Sees. 2.5, 




4 For completeness, the corresponding thermodynamic critical 
angular velocity for a weakly int eracting Bose-Einste i n gas 
is fl c /uj ± 1 - N a/jV&Kdz) <Butts and Rokhsarl. Il999l ; 

iLinn and Fettel ll99Sft . Note that this result also holds for 
attractive interactions (a < 0), but the resulting inequality 
Q c > implies that the radial confinement disappears before 
the vortex would be stable, in agreement with earl ier conclusions 
jDalfovo and StringariL Il996t IWilkin et aR Il998h . Indeed, at- 
tractive interactions quite generally tend to suppress BEC in fa- 
vor of "fragmented" states wit h macroscopic occupation in more 
than one sing le-particle state ijNozieres and Saint Jamej . 1 19821 : 
iPolloclj . Il967h 
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FIG. 1 Dimensionless scaled energy AE'(r , Q)/AE'(0, 0) 
[Eq. (|3.9[) ] of a singly quantized straight vortex in a rotat- 
ing disk-shaped TF trap as a function of the displacement ro 
from the symmetry axis. Different curves represent different 
fixed values of the angular velocity Q: (a) Q — 0; (b) Q — fl, n , 
defined in Eq. (|3dT|) ; (c) Q = fi c , defined in Eq. (pHO) . 



2. Dynamical motion of a trapped vortex 

Although the GP energy functional Eq. (|3.6p describes 
the equilibrium of a GP condensate with one vortex (or 
more vortices), it also provides a very useful dynamical 
picture of the vortex motion in a trapped condensate. 



a. Lagrangian functional The most direct approach relies 
on the basic observation that the time-dependent GP 
equation (|2.33p is the Euler-Lagrange equation for the 
time-dependent Lagrangian functional 



dV — [ty* — 
2 V dt 



~dt 



E'[9] (3.12) 



under variation of If the condensate wave func- 

tion \& depends on one or more parameters, the resulting 
Lagrangian functional provides approximate Lagrangian 
equations of motion for these parameters. In the context 
of BECs, this approach was first applied with great suc- 
cess to the low-lying normal modes of a trapped station - 
ary (nonrotating) condensate (jPerez-Garcia et"aZl . ll996lh 
Subsequently, its application to a single straight vor- 
tex in a disk-shaped condensate yields an explicit pre- 
diction for the angular precession of such a trapped 
vortex dLundh and Ad l2000t iMcGee and Holland!. 12001]; 



ISvidzinskv and Fetterl . l2000ah . In this case, the position 
ro of the off-center vortex serves as an appropriate pa- 
rameter, and the Lagrangian formalism shows that the 
precession rate is proportional to (the negative of) the 
slope of the appropriate curve in Fig. [1] For a nonrotat- 
ing disk-shaped TF condensate, this analysis gives 



2 ' 



(3.13) 



where Q m = § H ln(Rx/£)/(MFi]_) is the critical rota- 
tion frequency for the onset of metastability given in 
Eq. (|3.1ip . Note that a vortex in a nonrotating con- 
densate precesses in the positive sense, which is the same 
sense as the circulating vortex velocity field. 

In this Lagrangian approach, the factor 1 — TqJR\ in 
the denominator arises from the parabolic radial TF den- 
sity profile. It means that a vortex near the outer bound- 
ary precesses more rapidly than one near the center, al- 
though this simple picture omits bending of the vortex 
and thus fails as ro — * R±. 

It is interesting to compare this TF result [Eq. (|3.13p ] 
to the precession rate for a long straight classical vortex 



in an incompressible fluid inside a circular cylinder of ra- 
dius R. In this latter case, the boundary condition of 
zero radial flow velocity at the radius R requires an ex- 
ternal ima ge vortex, yielding the resu lt [see, for example, 
Sec. 9.4 of iPethick and Smith! (|2002h ] 



1 



£>cl = 



MR 2 1-rl/R 2 ' 



(3.14) 



Although the two expressions for the precession rate 
have the same denominators, the precession here arises 
from the motion induced by the image (and there is no 
large logarithmic factor). Since the TF density neces- 
sarily vanishes at the condensate radius, image vortices 
are generally omitted for trapped condens ates, although 



this question remains partially unre solved ( Anglinl . [2002] 
iKhawaja] [2005 iMason et a/.l . l2006t l. 



b. Matched asymptotic expansion An alternative ap- 
proach to the dynamics of a vortex in a trapped con- 
densate is to study the time-dependent GP equation 
(|2.33p itself; this analysis reli es on the method of matched 
asymptotic expansions dPismen and Rubinsteinl . 1 1991] 
iRubinstein and Pismenl . Il994[ ) [see also Sees. 2.2 and 
5.2 of iPismenl (Il999l)l. In the TF limit, the vortex 



core is small compared to other length scales, and it 
is pos sible to use the separation of length scales as fol- 
lows ( Svidzinskv and Fetterl . l2000bh . Assume that the 
straight vortex (located at ro) moves with a locally uni- 
form velocity V(ro) that lies in the xy plane. Transform 
to a locally comoving reference frame in which the vor- 
tex is stationary. The problem is solved approximately 
in two different regimes: 

1. Near the vortex core, the short-distance solution 
shows that the trap potential exerts a force pro- 
portional to Vj_T4 r evaluated at ro. This solu- 
tion includes both the detailed nonuniform struc- 
ture of the vortex core and the asymptotic region 

|r-r |. 

2. Far from the vortex, the core can be approximated 
as a line singularity and the resulting solution also 
includes the region £ <C |r — t*o|- 

The two solutions must agree in the common region of 
overlap, and a lengthy analysis determines the transla- 
tional velocity V(ro) of the vortex line near the center 
of the three-dimensional disk-shaped trap 

V(r ) = ^ In (^j z X VxKr(ro). (3.15) 

To understand the structure of this result, note that the 
trap exerts a force — Vj_Vt r on the vortex. Since the vor- 
tex itself has intrinsic angular momentum, it acts gyro- 
scopically and hence moves at right angles to this applied 
force. As a result, the motion follows an equipotential 
line along the direction z X Vj_Vtr and thus conserves 
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energy, as expected for the GP equation at zero temper- 
ature. 

In fact, this result is essentially the same as that in 
Eq. (|3.13[) for r$ near the trap center. The gradient of the 
axisymmetric potential Vtr = \M{ui\r 2 + w 2 z 2 ) and the 
squared TF condensate radius R\ — 2fj,/(M(Jj_) readily 
yield the equivalent form for Eq. (|3.15[) 



V(t-q) = H m r0, 



(3.16) 



involving the critical angular velocity Q m for the onset 
of metastability given in Eq. (|3 . 1 1 [) . These expressions 
for the translational velocity V can be generalized to 
treat a rotating anisotropic harmonic trap with uj x ^ 
ujy , and to include the three-dimensional effects of vortex 



curvature ( Svidzinskv and Fette r. 2000a,tJ. 

Although the final result in Eq. (|3 . 1 5[) is simple and 
physical, the details of the method of matched asymp- 
totic expansions are not very transparent. A more 
physical picture starts from the conservation of par- 
ticles V • (nv) — in the locally comoving coordi- 
nate system. In the TF limit, the particle density is 
explicitly known in terms of the trap potential [com- 
pare Eq. (|2.27p ]. Since v oc VS*, this equation can 
be rewritten as Vn • V5 + nV 2 S — 0, which deter- 
mines the phase S of the condensate wave function in 
terms of the known density. A perturbation expan- 
sion with the appropriate azimuthal angle [see Eq. (13. 
as the leading 
tional to In \ r 

ISvidzinskv and FetterT 2000bl ). The gradient of this ad- 
ditional phase evaluated at t*o then fixes the precessional 
velocity of the vortex (|3.15|) . In this latter perspective, 
the precessional motion can be thought to arise from 
the gradient of the condensa te density n(r) evaluated 
at the position of the vortex ( Sheehv and Radzihovskvl 
l2004alTbh 



as the leading term yields an additional phase propo r- 
r n l (ISheehv and Radzihovskvl l2004allb1 ; 



= h_ z X (-Vn) 

V 01 M 2n 



R± 



(3.17) 



where the minus sign appears because the density de- 
creases with increasing r (this detailed result describes 
a two-dimensional condensate). Apart from a numeri- 
cal factor that reflects the difference between two and 
three dimensions, Eqs. (|3.15[) and (|3.17[) describe the 
same physics because the gradient of the TF density Vn 
is proportional to — VVt r [see Eq. (|2.26l) ]. 



c. Bogoliubov equations and stability of a vortex Soon af- 
ter th e creation of BECs in dilute 87 Rb ( Anderson et all 
Il995l ). theorists studied the collective- mode spectrum of 



a nonrotating stationary vortex-free trapped condensate 
using the GP equation (|2.23| for the condensate and the 
Bogoliubov equations (|2.41[) and (12.421) for t he sm all- 
amplitude normal modes. As lEdwards et al] ( 19961 ) re- 
ported, the agreement between the theoretical predic- 
tions and the experi mental observa t ions w as impressive 
[see also Sec. IV. A of lDalfovo efal\ (|l999h ]. 



This success ra p idly led to a similar numerical 
study ( Dodd et all Il997l ) of the collective modes of 
a singly quantized vortex in an otherwise nonrotating 
trapped condensate. If the vortex is on the axis of 
symmetry, the condensate wave function is like that in 
Eq. (|3.7[) . and the conservation of the z component of 
angular momentum greatly simplifies the analysis. Re- 
markably, the spectrum of quasiparticle excitations for 
the vortex turns out to contain a very unusual mode (now 
called "anomalous" ) with positive normalization integral 
[Eq. (|2.45[) ]. negative angular momentum m a = — 1 rel- 
ative to that of the vortex, and negative frequency tu a 
Rokhsarl (|l997bl) has emphasi zed this particular identi- 
fication; see also iFetterl (jl998l )]. As noted in Sec. II.C 
[see Eq. (|2.46|) ]. such a negative eigenfrequency implies 
that a vortex in a nonrotating condensate is formally (en- 
ergetically) unstable, but the vortex cannot spiral out of 
the condensate in the absence of dissipation [compare the 
discussion concerning curve (a) in Fig. [J. Detailed stud- 



ies (|Linn and Fetterl . Tl999l : rSvidzinskv and Fetterl . l2000al ) 
show that io a w — ojj_ in the near-ideal limit and uj a w 
— J7 m in the TF limit, so that u> a remains negative for all 
coupling strengths. 

The physics underlying the anomalous mode is es- 
pecially transparent in the weak-coupling (near-ideal) 
limit, when the condensate with a singly quantized vor- 
tex at the origin has the radial wave function tfti (r, (ft) cx 
(x + iy) exp(— \r 2 ) ~ re 1 ^ exp(— \r 2 ), with one unit of 
angular momentum and an excitation energy fko±. Note 
that differs from the true radial ground-state wave 



function -00 oc exp(- 



which instead has zero angu- 



lar momentum and zero excitation energy. In principle, 
the particles in the macroscopic vortex condensate could 
make a transition to the nonrotating ground state, with 
a change in angular momentum m = — 1 and a change 
in energy —huj±; this transition is just the anomalous 
mode in an ideal Bosc-Einstein gas. Such anomalous 
modes typically occur whenever the condensate occupies 
a single-particle state that is not the true ground state 
[compare the discussion of a moving condensate below 
Eq. 

In the present case of a condensate with a single vor- 
tex, the anomalous mode has two very important conse- 
quences 

1. The linearized solution for the jth mode of the 
time-dependent GP equation in Eq. (|2.40[) imme- 
diately yields the corresponding perturbation in 



the density dr 



In par- 



ticular, a singly quantized vortex with 'J given in 
Eq. (|3.7|) has an anomalous-mode density perturba- 
tion 5n a (r, t) tx exp[i(m a (j) — Lo a t)], with m a = — 1 
and negative excitation frequency uj a = — \uj a \. 
This density perturbation Sn a oc exp[z(|w a |i — (ft)] 
precesses in the positive sense around the symmetry 
axis of the condensate at a rate — u> n = \uJn\- In the 
TF lim it when |c<j a | = f2 ro ( Svidzinskv and Fetterl . 
2000a), this rate a grees with that in Eq. (|3.13[) and 
with experiments ( Anderson et aLl . l2000f) discussed 
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below in Sec. III.D.l. 

2. When the condensate rotates, Eq. (|3.5p shows that 
the frequency uij of the jth mode in the laboratory 
frame shifts to 



rrijVL 



(3.18) 



in the rotating frame, where rrij is the correspond- 
ing angular quantum number. For the anomalous 
mode, this result implies that w (fi) = — \u) a \ + 
which becomes positive and thus energetically sta- 
ble when O exceeds \u> a \. This conclusion confirms 
the critical rotation frequency O m for the onset of 
metastability in the TF regime, as seen in curve (b) 
of Fig. CD 



C. Small vortex arrays in a trap 



In the TF limit, Eq. (|3.8j) shows that the condensate 
wave function has a 27r-phase singularity associated with 
a ■ 



vortex at Tq. More generally, a v ortex is simply a node 
. the condensate wave function dFevnmanl. 11955). In 
the co ntext of dilute ultracold BECs. lButts and Rokhsarl 
( 19991) exp loit this property i n thei r studies of small vor- 
tex arrays [Castin and Duml (|l999|) obtain similar vortex 
configurations with a somewhat different variational ap- 
proach] . 




FIG. 2 Angular momentum I vs. dimensionless angular veloc- 
ity fl/u>± for stable (black lines) and metastable (gray lines) 
states. Here, 7 = (2/n) 1 ^ 2 Na/d z is assumed small. There 
are no stable states with < I < 1, and I = 1 represents 
a single vortex on the axis of symmetry. For larger I, the 
angular momentum changes smoothly because the positions 
of the two or more vortices vary continuously with f2 (for 
each stable branch, the integer denotes the rotational sym- 
metry). Three-dimensional plots of the density for two and 
se ven vortices are remark ably similar to experimental images 
Yarmchuk el al\ lll979]) for superfluid 4 He (see Fig. [Sj and 
' ,n 2000ali for ultracold dilute 87 Rb gas (see 



of 
of 



-- .Madison et al.\ ... . 

Fig. I7]l. From lButts and RokhsaJ ljl999l ). 



To understand their procedure, assume that the con- 
densate is tightly confined in the z direction with a small 
axial oscillator length d z <C d±. The remaining two- 
dimensional isotropic harmonic trap has the one-body 
Hamiltonian H = \ (p 2 /M + Muj\r 2 ). For positive 
angular momentum, the low-lying normalized one-body 
eigenfunctions are 



(ttto!) 1 ^ d A 



exp 
exp 



u 
~1> 

' 2 

oi<t> 



(3.19) 



where u = r/d±, £ = (x + iy)/dj_ — ite 1 ^, and m > 0. 
Note that ip m is also an eigenfunction of the dimension- 
less angular-momentum operator I = —id/ dip with eigen- 
value to. As a simple example of a possible trial state, 
consider the linear combination = CoV>o + Ci^i, which 
vanishes at (0 = (xo+Wo)/d± — ~ c o/ c i (the onr y node). 
The point Co is the position of the vortex, since the phase 
of ^ increases by 2-7T on going around Co m the positive 
sense. 

For a r otating condensate with p ositive angular mo- 
mentum, iButts and Rokhsarl (|l999h use a general lin- 

Y^ m >0 C m^m With 



ear combination of such states \E' 
normalization Y^ m >o l Cm ' 2 = 



The expectation 
value of the angular momentum per particle is hi = 
hJ2m>o m \ c m\ 2 - Furthermore, the expectation value of 
the one-body Hamiltonian (H ) = hu± J2 m >o m \ c m\ 2 + 
huj± = fkuj_(l + 1) is essentially that of the angular mo- 
mentum. Thus all states with the same value of I have 
the same one-body energy, independent of the specific 
coefficients c m . To lift this degeneracy, it is necessary to 
consider the interaction energy E- lnt = \g J ci 3 r|\E'| 4 (as- 
sumed to be small). In the limit of tight z confinement, 
this energy reduces to a two-dimensional integral, which 
they evaluate numerically with the linear-combination 
trial state 'f '. For a given angular momentum /, the 
equilibrium state minimizes the total energy, which also 
means the i nteraction energy because of the degeneracy. 
In this way, IButts and Rokhsarl ( 19991 ) determine the co- 
efficients c m and thus the position of the vortices (the 
nodes of fy) as a function of I 

Given the equilibrium energy E to t{l), the thermody- 
namic relation Ml — dE tot /dl fixes the correspond- 
ing angular velocity fl. E xtensive numerical stud- 
ies ( Butts and RokhsarUl999T ) then determine the result- 
ing physical sequence of vortex states with increasing f2, 
as shown in Fig. [2l The state with I — 1 arises from 
a single vortex on the symmetry axis of the condensate, 
but configurations of two or more vortices vary with I, 
so that the corresponding angular momentum does not 
generally take integer values. 



Similar solutions play a crucial role in the study of rapidly ro- 
tating condensates in Sec. V. 
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D. Experimental and theoretical studies in a trap 

Two basic experimental methods have been most ef- 
fective i n creating vort ic es in a BEC. The first ap- 
proach ( Matthews et all . Il999bh manipulates two hy- 
pcrfinc components of 87 Rb, spinning up one compo- 
nent with a n external coherent el e ctrom agnetic cou- 
pling beam ( Williams and Holland . Il999| ). This ap- 
proach relies heavily on techniques from atomic physics, 
especially the optical Bloch equations for a two-level 
system [see, for exam ple, Chaps. 7-11 of Feynmanl 
(|l965lh Cha p IV of ICohen-Tannoudii et all (I1977T). 
Chaps. 1-3 oflAllen and Eberlvl (Il987f). and Chan. V of 
ICohen-Tannoudji et al. ( 19981)1. The second nearly si- 



multaneous approach ( Madison et all , l2000af ) is similar 
to the "rotating-bucket" method of conventional low- 
temperature physi cs [see, for example, Chaps. 2 and 5 
of lDonnellvi (|l99lD ]. In the present context of BECs, a 
driven rotating deformation of the confining trap serves 
to spin up the whole condensate. 



1. Vortex behavior in two-component BEC mixtures (JILA) 

The search for quantized vortices in dilute BECs 
started soon after the creation of trapped conden- 
sates. Since a mixture of two hyperfine components of 
87 R b play an essent i al role in the first reported vor- 
tex dMatthews et all Il999bh . I first review the physics 



of such systems. The GP energy functional in Eq. (|2.19|) 
is readily generalized for a mixture of two BECs with con- 
densate wave functions and \&2- The total energy con- 
tains (1) the two separate energies ^ [^j] of each compo- 
nent j = 1,2, with appropriate selfinteraction parameters 
gj (I assume repulsive interactions with gj > 0) and (2) 
the interaction energy £i 2 [*i, ^2] = f dV gi 2 \^i\ 2 |* 2 | 2 
between species 1 and 2, where gi2 is the interspecies in- 
teract ion parameter [se e Sec. 12.1 of[Pethick and Smith! 
(|2002h and Sec. 21.11 of lPitaevskii and Stringaril (|2003l )]. 
For the simplest case of two uniform BECs in a box, 
if g^2 > 5i<72, the two components are immiscible 
and separate into nonoverlapping phases. Otherwise, 
if g 2 2 < g\g2i they form two uniform interpenetrating 
BECs [in the context of dil ute BECs, this result a ppar- 
ently was first obtained by Nepomnvashchil ( 19741 ); see 
also (|Colson and Fetterl . Il978l )]. 

For trapped dilute BECs, the situation is more com- 
plicated because the nonuniform trap potentials Vj in- 
duce both kinetic and potential energies that affect 
the o verall energy functional. The first such experi- 
ment (|Mvatt et all , Il997| ) at JILA (Boulder, CO) uses 
two different hyperfine states of 87 Rb. This isotope has 
nuclear spin / = |, which combines with the single 
valence electron (S — 5) to give two hyperfine man- 
ifolds, a lower one with F = 1 and an upper one 
with F — 2. The typical trap has a local minimum 
in the magnetic field, and only some of these hyper- 
fine states are weak-field seeking and hence stable [see 



C haps. 3 and 4 of iPethick and Smith! (|2002ri and Chap. 9 
of iPitaevskii and Stringaril (|2003h for a description of 
magnetic and optical traps]. In the common notation 
\F,mp), the stable weak-field seeking states are [1, — 1), 
|2,l),and|2,2). 

The experiment of iMvatt et all ( 19971 ) cools state 
[1,— 1) evaporatively and the interspecies interaction 
cools state |2,2) through sympathetic cooling. Al- 
though the interaction paramet ers approximately satisfy 
S12 ^ 5iff2 ( Cornell et a/.l . fl998! ) . implying that two such 
uniform BECs would overlap, the trapped condensates 
in fact separate because of the difference s in the two 
trap potentials and interaction parame ters (IMvatt et all , 
I1997D . A careful theoretical analysis (jEsrv et all Il997l ) 
that includes all the relevant effects finds a reasonable fit 
to the experiments. 

To avoid the problem of different trap potentials, sub- 
sequent experiments on these mixtures use the pair of 
states |1) = |1,-1) and |2) = |2, 1) that differ by an 
angular momentum Am = 2. Each component can be 
imaged selectively with appropriately tuned lasers. Fur- 
thermore the state can be quickly changed from |1) to 
1 2) and back by a two-photon transition (microwave at 
w 6.8 GHz and rf at « 2 MHz). To understand the 
physics of this fascinating system, it is helpful to exam- 
ine the coupl ed time-dependen t GP equations for the two 
components ( Hall et all , Il998l ) 



ih 



dt V* 



(3.20) 



T + Vi + V m 
\Ml{t)e- 1 ^ 



\m(t) e iu < 

T + V 2 + V H 2 




where T = — fi, 2 V 2 /(2M) is the kinetic energy, Vj is the 
trap potential, and Vnj is the Hartree energy including 
the interaction with both densities \^\\ 2 and l^l 2 - The 
crucial new feature here is the applied electromagnetic 
(microwave and rf) fields with combined frequency w cm ; 
they produce an off-diagonal coupling that involves the 
phase of each condensate wave function. The effective 
coupling f2(t) is known as the "Rabi" frequency, which 
increases with the strength of the app l ied ele ctromagnetic 
fields [see Ch. 3 of lAllen and Eberlvl (|l987l )]. In general, 
this coupling depends explicitly on time (it vanishes when 
the fields are turned off). Finally, 8 is the detuning be- 
tween the combined frequencies of the driving fields and 
that of the atomic hyperfine transition |1) <-> |2). 

This electromagnetically coupled two-component sys- 
tem provides a remarkable example of topological con- 
trol of a quantum-mechanical state through an external 
parameter. Specifically, when the near-resonant elec- 
tromagnetic coupling Q (t) is turned off, the interac- 
tion between the condensates involves only the densities 
(through Vhj) with no phase information. In this case, 
each condensate has its own complex order parameter 
The associated U (l) symmetry [see, for example, 
Sec. 27 of lSchifj (|l968t) ] has the topology of a circle or a 
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cylinder. For T <C T c , the magnitude \^fj\ is fixed, and 
each complex order parameter has its own phase angle, 
yielding a quantized circulation that remains a topologi- 
cal invariant. 

The topology is very different when the electromag- 
netic coupling is turned on because the coupled con- 
densates in Eq. (|3.20p now have a single two-component 
SU(2) order parameter with the topology of a sphere, like 



that of a particle with spin \ [see Sec. 27 of lSchifj ( 19681) 



and (|Hall et all Il998h ] . Two spherical-polar angles suf- 
fice to characterize the coupled syste m, whose dynam- 
ics ob eys the optical Bloch equations (jAllen and Eberlvl . 
1987), similar to those of nuclear magnetic resonance 
(NMR). Specifically, the transformation of |1) into |2) is 
analogous to a tt pulse in NMR that rotates a spin from 
up to down. If the pulse is twice as long, the resulting 
2ir rotation essentially reproduces the initial |1) state up 
to a sign. 

An elegant experiment ( Hall et all . Il998h shows the 
ability to manipulate the coupled condensates as a single 
SU(2) order parameter. Application of a n/2 pulse to 
a pure state |1) produces a coherent mixture with equal 
population in both states that are spatially separated 
with a small region of overlap. After a variable delay, a 
second 7r/2 pulse measures the relative quantum phase 
between them. A subsequent experiment demonstrates 
the three- dimensional character of the SU(2) topology 
explicitly (jMatthews et all 1 1999at I Williams et aZ.I . l2000h 
as follows: An external coupling field with a spatial gra- 
dient produces a differential twist along the length of the 
condensate. With increasing time, the increasing phase 
difference between the two ends of the condensate relaxes 
periodically by a three-dimensional evolution of the cou- 
pled SU(2) order parameter, leading to a collapse and re- 
currence of the structure. Note that such behavior would 
not occur for two uncoupled U(l) condensate wave func- 
tions. 



La ser rela tion, to 
Laser — 
profile f \ 




(a) 



FIG. 3 (a) Schematic illustration of the technique used to 
create a vortex: An off-resonant laser provides a rotating 
force on the atoms across the condensate with a simultane- 
ous microwave drive of detuning S. (b) Level diagram show- 
ing the microwave transition to the sublevels of the state 
1 2). When the rotation frequency of the laser beam ui is 
close to the detuning S, the perturbation selectively modu- 
lates the transition from the 1 = state to the I = 1 state. 
From iMatthews eTal\ (|l999bl ). 

To create a quantized vortex, IWilliams and Holland! 



(1999) propose applying an external laser beam in addi- 
tion to the external electromagnetic fields that produce 
the two-photon transition. This external laser beam in 
effect adds a rotating gradient to the coupling fields with 
a frequency oj (see Fig. [3]). Such a process starts from the 
ground state |1) and yields the state |2) with one unit 
of angular momentum. Specifically, the theory predicts 
that the angular momentum of the state |2) oscillates and 
increases with time. Turning off the coupling at the ap- 
propriate time leaves the state \2) with a persistent singly 
quantized vortex because the associated order parameter 
for this state |2) now has U(l) symmetry. 




0.5 1,0 1.5 
azlmuthal angla 8 (i) 



FIG. 4 Cosine of the relative phase around the vortex, 
showi ng the sinusoidal varia tion of the azimuthal angle. 
From IMatthews et all (|l999bh . 



In the experiment ( Matthews et all . Il999b| ), the vor- 
tex in state \2) surrounds a nonrotating core consisting of 
state |1). Furthermore, a tt pulse can interchange |1) and 
1 2), so that the vortex can occur in either state. There 
is one important difference between the states |1) and 
|2), because the state |2) with = 1 is not maximally 
aligned. Consequently it can decay by spin relaxation 
with a measured lifetime ~ 1 s. As expected, a vortex 
in state |1) surrounding a core of state |2) is more sta- 
ble because the state |2) core eventually decays, leaving 
a simple one-component vortex in state |1). The experi- 
ment can image either the filled core or the surrounding 
vortex, but these images by themselves do not demon- 
strate the presence of quantized circulation. Fortunately, 
the interference procedure described above allows a clear 
measurement of the relative phase between the two con- 
densates at various positions, demonstrating the sinu- 
soidal variation predicted for a vortex (see Fig. 0|. 

One very important advantage of this two-component 
experiment is the ability to follow the motion of 
the vortex core with no ndestructive time-lapse im- 
ages (| Anderson et all , feOQOt ). Figure \E\ shows the uniform 
precession of the filled vortex core in state |2) around 
the trap center, determining the angular velocity of the 
circular motion. Note [Fig. GJd)] the gradual shrinkage 
of the core radius as state |2) decays through spin- flip 
transitions. The experimenters can also remove the core 
component with a brief resonant laser pulse, leaving a 
one-component vortex whose empty core (of order £) is 
too small to image in the trap. To study the precession 
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FIG. 5 (a) Successive images of a two-component condensate 
with a vortex in state |1) surrounding a core in state |2). (b) 
The image is then fit with a smooth Thomas-Fermi distri- 
bution, (c) The azimuthal angle of the vortex core plotted 
against the time, giving a precession frequency of 1.3(1) Hz. 
(d) The gradual decrease in the core radius r (measured in 
units of the coherenc e lengt h £) because of spin relaxation. 
From lAnderson etail (|200oh . 



of such a one-component vortex, the experimenters first 
image the two-component vortex, then remove the core, 
wait a variable time and finally turn off the trap, imaging 
the expanded vortex core. In this way, they measure the 
precession rate of a one-component empty-core vortex, 
yielding val ues in reasonable agreement w ith theoretical 
predictions (jFetter and Svidzinskvl . 1200 ll ). These data 
show no evidence that the vortex spirals outward, sug- 
gesting that thermal damping is negligible on the time 
scale of ~ 1 s. 



2. Vortex behavior in one-component BECs (ENS and 
subsequently others) 

Most other experimental methods of creating vortices 
in dilute BECs rely on the rotation of an anisotropic 
potential, and I briefly review the physics of such a 
system. Consider a classical incompressible fluid in a 
long cylinder with elliptical cross section \x\ < a and 
\y\ < b. When the container rotates with angular veloc- 
ity = Qz, the moving walls (in the laboratory frame) 
exert forces on the fluid. As a result, the fluid executes 
a purely irrotational motion with velocity i>i rr = V$j rr 
[see Eq. (|2.36[) ], where the associated velocity potential 
<I>i rr = flxy(a 2 — b 2 )/ ( a 2 + b 2 ) is proportional to xy 
[see pp. 86-88 of iLambl (|l945| )]. Remarkably, this irro- 
tational velocity field induces a nonzero angular momen- 
tum per unit length L irT = L sh (a 2 — b 2 ) 2 / (a 2 + b 2 ) 2 , 
where L s b is the corresponding angular momentum per 
unit length for a fluid that executes solid-body rotation 
with t> s b = f2 X r. Note that this irrotational angular mo- 
mentum occurs without quantized vortices. Such vortex- 
free angular momentum increases the critical ang ular ve- 
locity f2 c for vortex formation in superfluid 4 He ( Fetter] , 



11974 ). as lDeConde and Packard! (|l975l ) verify by experi- 
ments on rotating elliptical and rectangular cylinders. 

The same irrotational flow occurs in rotating 
anisotropic trapped BECs, where it appears in the phase 
of the condensate wave function S = M<fr- m /h [see 
Eq. (|2.36[) ] and leads to unexpected phenomena. For 
example, a sudden small rotation of the anisotropic 
trap potential induces an unusual oscillatory motion 
k nown as the "scissors mode" [s ee Sees. 14.3 and 14.4 
of iPitaevskii and String ari (2003) for a detailed descrip- 
tion of this and related experiments] . 




Laser 



FIG. 6 A cigar-shaped condensate is confined by an axisym- 
metric magnetic trap and stirred by an off-center far-detuned 
laser beam. The laser beam propagates along the long axis 
and creates an anisotropic qua drupole potential that rotates 
at an angular velocity fi. From [Madison et al\ (|2000bh . 



The experimental group at the Ecole Normale 
Superieure (ENS) in Paris uses such a direct ap- 
proach (see Fig. [6]) to create one or more vor- 
tices in a one-component 87 Rb cigar-shaped conden- 



apo 

sate ([Madison et all , l200 0a) . Starting with an axisym- 
metric condensate, they apply a toggled off-center stir- 
ring laser beam that slightly deforms the trap potential; 
the plane of deformation then rotates about the long axis 
with an applied angular velocity f2/(27r) < 200 Hz. Phys- 
ically, the laser beam acts on each polarizable atom in the 
BEC with a "dipole" force 6 proportional to the gradient 
of the squared electric field of the laser [for the de tails of 
the stirring procedure, see (|Madison et ali , I^OOObD ] . This 
dipole-force mechanism is the basis for optical trapping 
of atoms, as discussed briefly in Sec. III.E in connec- 
tion with optical l attices [for a more d etailed account, 
see Sec. 4.2 of iPethick and Smith! (|2002l ) and Sec. 9.4 of 
IPitaevskii and Stringari ( 2003f )]. 

Since the vortex core radius ~ £ w 0.2 — 0.4 /im is too 
small for direct visual detection, they turn off the trap 
and allow the condensate (including the vortex core) to 
expand. The initial tight radial confinement means that 
this expansion is mostly in the radial direction, and the 
expanded condensate usually becomes disk-shaped. The 
resulting images of the expanded condensate clearly show 
the presence of one or a few vortices (see Fig. [7]) . These 
latter images are remarkably similar t o those in Fig. EQ 
of vortices in rotating superfluid 4 He (jYarmchuk et all 



6 Here, "dipole" means the induced atomic electric dipole; the ac- 
tual deformation of the originally axisymmetric condensate usu- 
ally has quadrupole symmetry. 
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Il979l ) obtained by trapping charged p articles on t he vo r- 
tex cores [see Sec. 4.5 and Chap. 5 of iDonne El (EMI)]. 
Hodby et al\ (|2002l ) (Oxford) observe similar vortex pat- 
terns with a rotating magnetic deformation of the trap 
potential. Their experimental arrangement allows them 
to explore the physics of vortex nucleation for relatively 
large quadrupole deformations. 
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FIG. 7 Absorption images of a BEC stirred with a laser beam; 
in all images the condensate has about 10 5 atoms, the temper- 
ature is below 80 nK and the rotation rate Q/(2ir) increases 
from 145 Hz for (a) and (c) to 168 Hz for (g). Images (a) 
and (b) show the optical thickness for images (c) and (d), 
with the clear appearance of the vortex core. Images (e), 
(f), a nd (g) show states wi th two, three, and four vortices. 
From IMadison efal\ (|2000af ) . 
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FIG. 8 Photographs of stable vortex configurations in su- 
perfluid 4 He, with u p to four vortices. Adapted from 
lYarmchuk et al\ [| 19791 ). 



Careful experimental st udies ( Chew et all l2000t 
IMadison et aZ] . l200lL [2000al [bh show that the critical an- 
gular velocity fi c for the appearance of the first vortex is 
« 0.7wj_, which is typically much larger than that pre- 
dicted in Eq. (|3.10[) for the relevant TF limit. In fact, 
a rotating axisymmetric condensate has a quadrupole 
instability at an angular velocity f2 = luj_/\/2, which 
agrees well with the measured value for the appear- 



ance of the first vorte x. Sections 14.5 and 14.6 of 

iPitaevskii and Stringaril ( 20031 ) discuss this intriguing be- 
havior in detail. 



3. Normal modes of an axisymmetric condensate with one 
central vortex 

A static condensate without vortices has numerous 
low-lying small-amplitude normal modes characterized 
by coupled density and velocity perturbations. Many 
experimental studies have confirmed the various theo- 
retical predictions in cons i derabl e detail [see pp. 178- 
195 of IPethick and Smlthl (l2002f) and pp. 177-190 of 
IPitaevskii and Stringaril (|2003f)l. A direct tr eatment of 
the hydrodynamic equations (jstringaril . [l996) for a non- 
rotating axisymmetric condensate in the TF limit yields 
the frequency of quadrupole modes with I = 2, to = and 
I = 2, to = ±2. The last two modes with to = ±2 are de- 
generate and together produce an oscillating quadrupole 
distortion of the condensate in the plane perpendicular to 
the original axis of symmetry. In the absence of any per- 
turbation that breaks time-reversal symmetry, the axes 
of the deformation remain fixed in space [see Figs. OJa)]. 
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FIG. 9 Nondestructive images of quadrupole mode of an 
axisymmetric condensate. Condensate (a) is vortex-free, 
whereas condensates in (b) and (c) have straight vortices 
(with opposite sense of circulation) aligned perpendicular to 
the plane of the figure; (d) shows a linear fit to the precession 
angle. From lHalian et al\ (|2001al ). 

When the condensate has a vortex, how- 
ever, its circulating velocity breaks time-reversal 
symmetry and splits the previously de ge nerate 
states with ±m (ISvidzinskv and Fetterl . Il998t 
IZambelli and Stringaril Il998l ). This effect is analo- 
gous to the Zeeman splitting of magnetic sublevels. In 
the usual TF limit, the fractional splitting is small (of 
order \m\d\ / R\) . One dramatic consequence is that 
the quadrupole distortion now precesses slowly in the 
direction of the c i rculati ng vortex flow (|Chevv et ali , 
\200d lHalian et aZ.I . l2001ab . as seen in Figs.^b) and (c). 
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IZambelli and Stringaril ( 19981 ) note that the magni- 
tude of the splitting provides a direct measure of 
the an gular momen t um o f the condensate. The ENS 
group (|Chevv et ali . l2000h use this approach to detect 
the first appearance of a vortex in the condensate and 
to verify that its angular momentum per particle is ~ h. 
It is remarkable that this macroscopic experiment with 
N ~ 10 6 atoms measures a single-particle fundamental 
constant; it demonstrates the coher ence of the ma cro- 
scopic wave function [see Sec. 4.7 of ILeggettI (|2006l) for 
a valuable per spective o n this and related experiments]. 
Section 2.4 of IDonnellvl ( 199ll ) discusses similar earlier 
macroscopic experiments on rotating superfluid He. 

The precession of the quadrupole mode provides a sen- 
sitive vortex detector (sometimes called "surface-wave 
spectroscopy") that complements t he direct visibility o f 
the vortex core. The JILA group ( Halian et ali . l2001af) 
use this method to track the large-amplitude tilting mo- 
tion of a vortex in an axisymmetric condensate. As 
a tilted vortex rotates about the axis of symmetry, its 
empty core initially moves away from the line of sight. 
Periodically, however, the vortex returns to its original 
orientation, when the optical visibility of the core and 
the quadrupole precession both attain their initial val- 
ues. The experiment indeed observes two recurrences, 
yielding g ood agreement with the predi cted precession 
frequency ( Svidzinskv and Fetterl . l2000bl ) . 



4. Bent vortex configurations 



When a disk-shaped condensate rotates faster than the 
critical velocity Q m for creation of a metastable vortex 
[see Eq. p. lip ], the initial equilibrium vortex is essen- 
tially straight. For condensates with larger aspect ra- 
tios R z /R±, however, the equilibrium vortex shape can 
be curved, depending on the angular velocity and the 
specific aspect ratio. Predictions of this phenomenon 
arise b oth from linearized stability analysis for a straight 
vortex (|Feder et ali . l200lk ISvidzinskv and Fetterl . l2000bD 
and from direct calculation of the state that minimizes 
the energy in the rotating frame (lAftalion and Rivierel 



120011: iGarria-Ripoll and Perez-Garrial l2001allbl ). ka 
long cigar-shaped condensate, for example, the initial 
vortex can reduce its total energy by forming a small 
symmetrical loop near the equator. As the angular veloc- 
ity fl increases, the vortex loop expands toward the axis 
of symmetry. Wh en these vortex lo o ps be come large, 
ENS experiments (Roscnb usch et ali . |2002| ) show that 
they can have both symmetric configurations (denoted 
"U") and antisymmetric configurations [denoted "(un- 
folde d) N" although "S" or "J" might be more appropri- 
ate]. iModugno et ali (2003) give a nice physical picture 
of how curvature can lower the total energy of the vortex. 



5. Kelvin (axial) waves on a vortex 

A long straight classical vortex line in an incompress- 
ible fluid represents a dynamical equilibrium between the 
kinetic-energy density and the pressure, as fol lows from 
Berno ulli's equation. In an influential article, iThomsonl 
(1880) (subsequently, he became Lord Kelvin) studied 
the small oscillations of such a vortex; he found a char- 
acteristic long- wavelength (|fc|£ <C 1) dispersion relation 



hk\ ( \ 

m — — 

2M \\k\£ 



(3.21) 



where k is the axial wave number, £ is the vortex core 
radius, and I assume a quantized circulation k — 2irK/M. 
For general k, these classical Kelvin modes are left- 
circularly polarized helical waves with velocity potential 
oc exp [i(kz — (/) — uJkt)] K\{kr), where (r, (j), z) are cylin- 
drical polar coordinates and K\(u) is a Bessel function 
of imaginary argument that vanishes for large u. Each 
element of the vortex core executes a circular orbit with 
a sense opposite to the circulating velocity of the vortex 
line. In a quantum-mechanical version, these modes have 
quantized angular momentum with m = —1. Note that 
this negative helicity (m — —1) is independent of the di- 
rection of the wave propagation (±fc). The overall helical 
behavior reflects the gyroscopic character of the vortex 
line, and the physics thus differs significantly from that 
of a stretched string. 7 

Unfortunately, Kelvin's treatment is rather intricate, 
for he allowed the unperturbed axisymmetric flow to have 
an arbitrary radial dependence (including both solid- 
body rotation and irrotational vortex flow). For ac- 
ces sible discu s sions of these Kelvin mod es, see S ec. 1. 7 
of IDonnellvl (Il99lh. pp. 337-344 of iFetterl (|l969ft . 
pp. 138-139 of iFetter and Waleckal (l2006h . Sec. 29 of 
iLifshitz and Pitaevskiil (|l980bl ) arid Sec. III.D of ISoninl 
(|1987D . 

In a modern quantum-mechanical context, I Pitaevskiil 
(|1961| ) used the Bogoliubov equations (j2"H"j) and (j2Tl2"|) 
to show that a singly quantized vortex in a dilute BEC 
also has this Kelvin mode with the same long-wavelength 
dispersion relation [Eq. (|3.21[) ]. Although the existence 
of Kelvin waves are widely accepted, direct evidence for 
their presence in low-te mperature s uperfl uids is relatively 
meager. Sectio n 6.1 of IDonnellvl (|l99lh and Sec. III.E 
of ISoninl ( 1987t l discuss the interpretation of relevant su- 
perfluid 4 He experiments in some detail. 

The evidence is much more direct in dilute BECs. 
iBretin et al. I (I2003h report a clear signature of such quan- 
tized Kelvin modes. They drive the condensate with a 
near-resonant quadrupole excitation (m — ±2) and ob- 
serve preferential damping of the mode with m = —2. 
Recall that the Kelvin modes have quantized angular 
momentum m = — 1 relative to the circulation around 



For a detailed comparison, see Fet ter and Harvevl dl97lT >. 
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the vortex, independent of the direction of k. Thus the 
authors suggest that this particular excited state with 
m = — 2 can decay by the creation of two Kelvin- wave 
excitations, each with m = — 1 and opposite momenta 
±hk along the axis determined by the conservation of 
energy. In contrast, the other quadrupole mode m = +2 
has no such open chan nel for decay. A theore tical anal- 
ysis of the line shape (jMizushima et all 120031 ) confirms 
this interpretation. Furthermore, transverse images of 
the expanded condensate excited with the quadrupole 
mode m = —2 exhibit the expected periodic spatial 
structure associated with the Kelvin-mode deformation, 
whereas t hose excited with m = +2 have no periodic 
structure (|Bretin et a/1 l2003h . 



6. Other ways to create vortices 

As seen throughout Sec. III.D, experimentalists have 
devised many procedures that successfully create one or 
more vortices in a dilute BEC. Section III.D. 1 discusses 
dynamical phase-imprinting in a two-component conden- 
sate, and Sec. III.D. 2 considers the more common ap- 
proach of rotating an anisotropic potential that spins 
up the condensate. A related technique (discussed be- 
low in Sec. IV. B in connection with vortex arrays) is 
to spin up the normal gas and the n cool into the su- 
perfluid state (H alian et all l2001b[ ). One-dimensional 
solitons provide an additional unusual alternative, for 
they generally decay to form three-dimensional vortex 
rings ([Anderson et all l200lt iDutton et all l200lh . 



a. transfer of orbital angular momentum from a laser A 
Gaussian laser beam can be decomposed into com- 
ponents with orbital angular momenta m relative to 
the propagation vector (the Laguerre-Gaussian modes). 
Thus a laser beam can, in principle, create a vor- 
tex by coherently transferring orbital angular momen- 
tum to the c o ndens ate. In a recent experiment, 
lAndersen et al. (200G) use a two-photon stimulated Ra- 
man process involving a counter-propagating linearly po- 
larized Laguerre-Gaussian mod e with m = 1 [see, for ex- 
ample (lAllen et aZ.l . l2003lll999h ] and a linearly polarized 
Gaussian mode with m = (the linear polarization en- 
sures that "spin" angular momentum does not cause the 
observed effect). Some of the atoms in an initial single 
23 Na condensate absorb a photon with one unit of an- 
gular momentum from the Laguerre-Gaussian beam and 
emit a photon to the Gaussian beam, acquiring both an- 
gular momentum and linear momentum. Consequently, 
these atoms form a new coherent moving condensate with 
a singly quantized vortex. After release from the trap, 
the resulting spatial separation between the new and the 
original condensates permits clear imaging of the vortex 

core. 

lAndersen et ali (2006) also use a sequence of such 
pulses to make vortices with higher quantized circula- 



tion; as noted below, however, such vortices are essen- 
tially unstable with respect to decay into singly quan- 
tized vortices. Recently, similar techniques create per- 
sistent quantiz ed flow wi t h one unit of circulation in a 
toroidal trap ( Rvu et all . l2007t k the ~ 10 s lifetime of 
this persistent flow arises from experimental factors and 
is not intrinsic. 



b. merging of multiple BECs Another interesting exper- 
iment creates a vortex by merging three 87 Rb conden - 
sates in a single disk-shaped trap. IScherer et ali ( 20071 ) 
apply a repulsive (blue-detuned) laser beam shaped as 
a three-fold optical barrier, dividing the noncondensed 
gas into three non-overlapping regions (see Sec. III.E for 
the effect of external laser beams). As the gas cools be- 
low T c , three independent and uncorrelated BECs form. 
The optical barrier is then slowly removed, and the three 
BECs merge. After a variable delay time, they remove 
the trap and image the expanded condensate. Depend- 
ing on the initial relative phases of the three condensates, 
a net circulating current can form, producing a clearly 
visible vortex core in the expanded image. For random 
phase differences, a vortex should appear with probabil- 
ity 0.25, in reasonable agreement with the experimental 
observations for repeated trials. 



c. imposition of topological phase As a final topic in this 
discussion of forming a single vortex, I consider the cre- 
ation of a doubly q uantized vortex through imposition of 
topological phase (|Leanhardt et all I2002D. The e xper- 
iment follows sugge stions by llsoshima et al. (2000) and 
lOgawa et all (|2002) . It relies on a special form of mag- 
netic confinement known as the Ioffe- Pritchard trap [for a 
detail ed description, see Sec. 4.1.3 of iPethick and Smitbl 
(2002)]. Tere, the important feature is that the x and 
y components of the magnetic field have a quadrupole 
form, whereas the z component can be taken as uniform. 

This MIT experiment studies 23 Na in the lower hyper- 
fine state \F = l,mj? = —1). Initially, the z component 
of magnetic field is large and positive, which aligns the 
hyperfine spin F vector with the direction of B w B z z . 
Thus this spin texture is uniform and oriented, with vari- 
able number density induced by the nonuniform trap po- 
tential. The experiment then slowly reverses the z com- 
ponent of magnetic field, and the hyperfine spin vector F 
adiabatically follows the local magnetic field. In particu- 
lar, when B z vanishes, the hyperfine spin vector F takes 
the (nonuniform) quadrupole form of the remaining x 
and y components of B. As B z then becomes negative, 
the topological winding of the spin texture continues adi- 
abatically, leading to a nontrivial final texture. Remark- 
ably, at the end of the B z inversion, a detailed analysis 
shows that the local spin texture has an additional factor 
e -i2m F 4> ^ wnere mp — _\ Consequently, this final state 
represent s a doubly quanti z ed vor tex with quantum num- 
ber two. iLeanhardt et al. (2002) use surface-wave spec- 
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troscopy (precession of the quadrupole oscillation mode, 
see Sec. III.D.3) to verify that the vortex indeed has an- 
gular momentum ss 2H per particle. 

The energy of a vortex arises principally from the total 
kinetic energy J dV ^Mv 2 \^>\ 2 , which is proportional to 
the square of the circulation. In a multiply connected 
geometry like an annulus, such flow represents a macro- 
scopic quantized circulation, and the presence of the in- 
ner wall ensures that the hydrodynamic flow generally re- 
mains stable (a persistent current) as long as the temper- 
ature is well below T c . The situation is very different for 
a free-standing multiply quantized vortex, which should 
be unstable with respect to splitting into the appropri- 
ate number of singly quantized vortice s. As mentioned at 
the end of Sec. II.C. |Pu et all (|l999h study the Bogoli- 
ubov equations for a doubly quantized vortex and find 
imaginar y eigenvalues that imply a dynamical i nsta bility 
[see also (|Castin and Duml . ll999l) ]. IShin et all (|2004h in- 
vestigate the stability of these doubly quantized vortices, 
with convincing images of the relatively large initial vor- 
tex core decayi n g into two smaller cores in « 75 ms. 
iMottonen et all ( 2003F) provide a theoretical analysis of 
the MIT experiments, emphasizing the unstable eigen- 
values of the Bogoliubov equations. 

An al ternative picture of the same topological tra nsfor- 
mation ( Leanhardt et all. 20021 ; lOgawa et all . 2002) relies 
on the idea of Berrvl l 1984 ) phase that arises from an adi- 
abatic alteration of the physical system (in this case, the 
overall direction of the magnetic field). In the context 
of quantum mechanics, many recent textbooks provide 
a careful de scripti o n of the Berry ph ase, for e x ample , 
Sec. 10.3 of lAberd (I2004D. Sec. 10 .2 of iGriffithd (|2005h . 
and pp. 592-607 of lShankarl (jl994l ). 



E. Laser trapping and optical lattices 

The response of an atom to elect romagnetic radia- 
tion is an ol d subject [see, for examp l e. I Allen and EberTvl 
(|l987h and ICohen-Tannoudii et all (l998|)], yet it has 
many important applications to the interaction between 
ultracold dilute atoms and laser beams. To understand 
the basic physics, place an electric dipole with moment 
p in an external electric field E. Standard electrostatics 
yields the interaction energy J7j n t = P • E, where the 
minus sign means that the lowest-energy configuration 
has the dipole moment oriented along E. 

This expression holds for a fixed dipole moment p, 
but an atom typically acquires its electric dipole mo- 
ment as a result of induced polarization. In such a case, 
the dipole moment is proportional to the applied field, 
with p = ae E in SI units, where e ~ 8.85 x 10~ 12 
farad/m is the permittivity of empty space and a is the 
atomic polarizability. As a very simple classical model, a 
grounded conducting sphere of radius b has a polarizabil- 
ity a = Airb 3 (three times its volume), and this picture 
provides a useful order-of-magnitude estimate for atomic 
polarizabilities. The total interaction energy of the atom 



with the external field then builds up from Ei 
integration yields the final expression 



\ae a E 2 



0, and 



(3.22) 



If the electric field E(r) is static and nonuniform, then 
Uint(r) = — ^ae E 2 (r) is also nonuniform. Thus an 
atom experiences a force that seeks to minimize the in- 
teraction energy, placing the atom at the maximum of E 2 
because the static polarizability is positive. This result 
also holds for low frequencies, but the situation eventu- 
ally changes because of induced transitions among atomic 
energy levels. The Drude model of a bound electron with 
resonant frequency Wq and damping time r subject to an 
oscillatory driving field yields the frequency-dependent 
polarizability 



a{u>) 



a{0)u 2 



iijj/r 



(3.23) 



where a(0) is the (positive) static polarizability. Note 
that the real part of a{ui) changes from positive to neg- 
ative when the applied frequency u> exceeds the reso- 
nant frequency u>q. The combined relation Ui n t(u) from 
Eqs. (|3.22p and (|3 . 23[) is sometimes called the ac Stark 
effect. 

These ideas apply directly to the optical trapping of 
atoms by a focused laser. If the frequency ui of the laser 
is less than that of the atomic resonance ("red detun- 
ing") with Rea(w) > 0, then Eqs. (j3~2"2"j) and (|3~2"3"]) 
indicate that the atoms preferentially move toward the 
squared-field maxima, whereas if the laser frequency is 
higher than the resonant frequency ("blue detuning") 
with Rea(w) < 0, then the atoms move toward the 
squared-field minima. In practice, it is necessary to 
use far-off-resonance lasers to avoid dissipation. For 
the alkali-metal atoms of interest here (Li, Na, K, Rb, 
Cs), the lowest s to p transition of the valence electron 
predominat es, as in the familiar yello w emission spec- 
trum of Na. IStamper-Kurn et~aJ\ ( 19981 ) successfully trap 
a pre-existing 23 Na BEC with a focused infrared laser 
beam. Unlike magnetic traps, which confine only cer- 
tain hyperfine states, these optical traps confine all such 
states. This experimental technique allows the study of 
"spinor" condensates (see Sec. VII. B. 3) that involve cou- 
pled BECs, one for each hyperfine component of a given 
F multiplet ((Hoi . Il998t Ibhmi and Machidall 19981) . 

A more subtle example is the formation of an optical 
standing wave to trap atoms in a periodic array. Con- 
sider a laser beam propagating along the z axis with 
electric field E(z,t) = Eocos(kz — ut). If the laser 
beam reflects normally off a mirror at z = 0, the re- 
sulting standing wave Eo [cos(kz — cot) — cos(fcz + cot)} = 
2Eq sin kz sin cot yields the corresponding squared electric 
field \E(z,t)\ 2 = 4|£ | 2 sin 2 kz sin 2 cut. This periodic lat- 
tice has a spatial period ^A, where A = 2jr/k is the wave- 
length of the original laser beam. For either sign of de- 
tuning, the atoms assume a spatially perio dic configura- 
tion in this optical lattice, as emphasized by Jaks ch et all 
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(|!998h . This seminal idea has wide applications in the 



study of ultracold dilute g discussed, for exam- 

ple, in Sec. 9.4 and C hap. 16 of Pitaevskii and Stringarl 
(|2003l ) and Sec. IV of iBloch et all (|2007l ). The strength 
of the laser field can be controlled externally, so that 
the height of the barrier for tunneling between adja- 
cent layers can vary from small to large. This pos- 
sibility leads to remarkable experiments that demon- 
strate the transition from superfluid to ins ulating behav- 
ior in three-dimensional optical lattices ( Greiner et all . 



2002) . In the present context o f one-dimensional lattices, 



Anderson and K ascvich ( 19981 ) are the first to demon- 



strate such atom trapping. Their vertical lattice expe- 
riences a gravitational potential difference MgX/2 be- 
tween adjacent layers, with associated interlayer tun- 
neling at a measured ac J osephson frequency wj w 
MgX/(2h) (|Josephsonlli1)62h . 

A particularly interesting situation arises when a ro- 
tating condensate is subsequently confined in a one- 
dimensional optical lattice oriented along the rotation 
axis. Consider an axisymmetric condensate with a sin- 
gle vortex on the axis of symmetry. Then adiabatically 
turn on an applied standing-wave laser field. If the laser 
beams propagate along the symmetry axis, the optical 
lattice slices the rotating condensate into many circular 
disks. The resulting periodic atomic array forms a set 
of effectively two-dimensional condensates that are cou- 
pled by tunneling between adjacent layers. Each two- 
dimensional condensate contains a segment of the previ- 
ously continuous vortex line. Thus it becomes a single 
pancake vortex, and the position of each vortex center 
can fluctuate more-or-less independently, depending on 
the strength of the la ser field and the height of r esult- 
ing potential barriers. iMartikainen and StooJ ( 2003T ) use 
a variational trial function to study the resulting axial 
modes in this discrete system, which are the analog of 
the Kelvin modes studied in Sec. III.D.5. At present, no 
direct experimental evidence for such behavior exists. 



F. Berezinskii-Kosterlitz-Thouless transition in dilute BECs 

Another remarkable application of an optical lat- 
tice is the study of the equilibrium state of a two- 
dimensional dilute Bose-Einstein gas. Although such 
a uniform two-dimensional system c annot have a true 
condensate with lon g -rang e order ( Hohenbergj 19671: 



M ermin and Wagnerl. 



Berezinskiil 



196 6]), neve r theles s „ 
(|l97lh and iKosterlitz and Thouless] (|l973t ) predict that 
it can indeed become superfluid below a certain critical 
temperature Tbkt- Above this temperature, free single 
vortices can form easily. The motion of such vor tices dis- 
rupts any organized long-range phase coherence (jLeggettl . 
l200ll ) and hence destroys the associated superfluid flow. 
Below Tbkt, in contrast, free vortices do not occur, al- 
though bound pairs of a vortex and an antivortex do exist 
[in the present context of di l ute ul tracold gases, Sec. 17.5 
of iPitaevskii and StringarJ (|2003f ) discusses this behav- 



ior in more detail]. As the temperature increases, so 
does the size of the bound pairs, and they eventually be- 
come unbound at Tbkt, producing f ree vortices. Experi- 
ment s with thin superfluid 4 He films ()Bishop and Reppvl . 
Il978l) and with Josephso n-coupled superconducting ar- 
rays (|Resnick et a/.l . ll98ll) have both confirmed these pre- 
dictions in considerable detail. 

In a recen t ex periment at ENS (Paris), 
lHadzibabic et at I (|2006l) apply a standing- wave laser field 
perpendicular to the symmetry axis of a cigar-shaped 
condensate, splitting the atoms into two essentially 
independent thin pancake strips with lateral dimensions 
ph 120 /Ltm x 10 («n. See Fig.[6]for the geometry (although 
there is no rotating laser beam in the present case): 
the laser beams propagate in the x direction and the 
pancakes are parallel to the yz plane. After these two 
isolated Bose-Einstein gases reach thermal equilibrium, 
the confining traps are turned off. Both pancakes then 
expand primarily in the x direction perpendicular to 
their planes and overlap spatially. 

Soon after the first observ ations of BECs in d ilute 
trapped gases, the MIT group ( Andrews et aLl . [l997h cut 
a cigar-shaped condensate in half with a blue-detuned 
laser, creating two separate three-dimensional conden- 
sates. When the trap is turned off, these two condensates 
overlap and produce clear straight interference fringes. 
Similarly, the ENS experiment here looks for and finds in- 
terference fringes between the two pancakes. At low tem- 
perature [see Fig fTUT a)]. the fringes are indeed straight, 
indicating constant relative phase. As the temperature 
increases, however, the fringes become wavy [Fig llOf b)]; 
furthermore, as seen in Fig llOf c). they sometimes contain 
one or more dislocations, which indic ates the presence of 
free vortices in one or bot h pancakes ( Simula and Blakid . 
120061 ; IStock et all , l2005bl ). They repeat the experiment 
at various temperatures and use the presence of free vor- 
tices to determine Tbkt- The res ulting value agrees well 
with the theoretical pred i ctions of lBerezinskiil <| 1971!) and 
IKosterlitz and Thoulessl ( 1973f) . although the finite size 
of the present samples yields a crossover instead of a 
sharp transition. 



(a) 



FIG. 10 Interference fringes between two pancake conden- 
sates (a) at low temperature, (b) at high temperature showing 
wavy interference fringes without vortices, (c) fringe disloca- 
tion indicat es the presence of a vortex in one or both conden- 
sates. From iKriieerl i|2007l ). 



ISchweikhard et all (|2007h report an intriguing study 
of vortex formation in a two-dimensional array of 
Josephson-coupled BECs. They start at T below T c with 
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a disk-shaped condensate along with a residual normal 
cloud and then apply a hexagonal optical lattice, trans- 
forming the system into an array of coupled BECs that 
act like Josephson junctions. Each condensate has a sin- 
gle effective phase, and the coupling energy between ad- 
jacent condensates has the form J(l — cos Acf>) where A<j) 
is their relative phase. Equivalently, the tunneling cur- 
rent between any pair of condensates is proport i onal t o 
JsinA(j > [for a careful discussio n, see Josephsonl ( 1962f ), 
iLeeeettl (120011) . and Chap. 6 of lTinkhaml (119961) ]. Vary- 
ing the amplitude of the optical lattice through the laser 
intensity tunes the coupling parameter J. 

The experiment starts by inferring the temperature T 
from a nondestructive image and then ramps down the 
o ptical lattice, a l lowin g the adjacent BECs to merge. As 
in lScherer et al. (2007), discussed in Sec. III.D.6, vortices 
can form depending on the relative phases, appearing as 
holes in the expanded gas clou d after the trap is t urned 
off. For a given temperature T, ISchweikhard et all ( 20071 ) 
measure the vortex density V near the center of the cloud 
as a function of the Josephson coupling parameter J. For 
two different temperatures, a plot D as a function of J/T 
shows nearly universal behavior: T> is negligibly small 
for J/T 3> 1 (low temperature), but it rises rapidly and 
saturates for J/T <C 1 (high temperature), as expected 
for the BKT vortex unbinding crossover. 



IV. VORTEX ARRAYS IN MEAN-FIELD 
THOMAS-FERMI (TF) REGIME 

A rapidly rotating condensate has a dense a r ray of vor- 
tices, with a uniform areal density ( Fevnmanl , 1 19551 ) 



MCI 
irH ' 



(4.1) 



as follows from these arguments: 

(1) Consider a steady classical incompressible flow with 
velocity v(r) in a rigid container that rotates with angu- 
lar velocity CI. If the fluid has number density n(r), then 
the relevant part of the energy in the rotating frame is 



E' = J dV (\Mv 2 - MCI -rxv)n. 



This result is equivalent to 



E' = ±M f dV [v - CI X r) 2 n - E. 



Sb : 



(4.2) 



(4.3) 



where E s h = \M J dV \Cl X r\ 2 n = ^/ s b^ 2 is the rota- 
tional energy for solid-body rotation with the correspond- 
ing moment of inertia I s h- The first term of Eq. (|4.3[) 
is non-negative, and hence the absolute minimum of 
E' occurs when v = CI X r, namely solid-body rota- 
tion = CI X r; in addition, the corresponding mini- 
mum value of E' min has the classical value —E^. Since 
V X v s b = 2CI, the irrotational superflow cannot assume 
such a solid-body rotation, but it seeks to mimic this flow 
through the nucleation of quantized vortices. 



The angular momentum L z 
thermodynamic relation 



(r xp) follows from the 



-(dE'/dCl) N , 



(4.4) 



[the Hellmann-Feynman theorem gives the same 

result. See, for e xamp le, pp. 1192-1193 of 

ICohen-Tannoudii et all ( 19771): for a brief history 
of this theorem, see p. 1 of lBrownl (|2000l )] . Since the 



actual E' exceeds the classical value — i/ s b^ 2 by a 
positive correction that depends on CI, the equilibrium 
value of the angular momentum necessarily differs from 
the classical value I s bCl. Evidently, this system exhibits 
a non classical moment of inertia [see Sec. 2.5 of ILeggettl 
(2006) for a discussion of these quantum-mechanical 
effects]. In addition, it also has transverse shear modes 
(the Tkachenko modes, see Sec. IV. B. 4), and it "melts" 
for sufficiently large angular velocities to form a "liquid" 
state (see Sec. VI). Whether this system should be 
considered a "supersolid" in the conyention al sense 
[see, for example, Sec. 8.3 of iLeggettl ( 20061 )] merits 
additional study. 

(2) For a superfluid, Eq. (|3 . 3[) shows that a singly quan- 
tized vortex at a position T"o has a localized vorticity 
V X v = (2-kK/M) (r — ro) z. In a rotating super- 
fluid, the vortices will be uniformly distributed to ap- 
proximate the uniform classical vorticity. Take a contour 
C containing J\f v such vortices in an area A; the total 
circulation around C is Tc = (2irh/M)N' v . Stokes's the- 
orem shows that the corresponding classical value would 
be r^ 1 = 2CIA, and comparison immediately yields the 
mean areal vortex density n v — Af v /A = MCI/(ttK). The 
inverse n~ x is the area per vortex irh/ (MCI) = nl 2 , which 
defines the radius 



I = 



MCI 



(4.5) 



of an equivalent circular cell. Thus the intervortex sep- 
aration is w 21. Note that I decreases with increasing CI 
(it is also the usual oscillator length for a frequency CI). 



A. Physics of BEC in axisymmetric harmonic traps for 
rapid rotation 

The Feynman relation (|4.ip shows that the number 
of vortices increases linearly with Ct assuming that the 
geometry remains fixed. Although this picture indeed 
applies to superfluid 4 He in a rotating bucket, the situa- 
tion is quite different for a dilute gas in a harmonic trap, 
because the centrifugal forces expand the condensate ra- 
dially (number conservation means that the condensate 
also shrinks along the axis of rotation). As a result, the 
total number of vortices in a given condensate 

Af v (Cl) w n v TTR 2 ± (Cl) = MClR\{Cl)/h = Rl(Cl)/l 2 (4.6) 

increases faster than linearly with CI. 
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To quantify the dependence of Rj_ on f2, it is conve- 
nient to focus on the TF limit of a large condensate, when 
the spatial variation of the condensate density is negligi- 
ble compared to the other terms in the energy functional. 
Thus = -ihV^/ = (WS)$ - ihe iS V\$>\ « Mof, 
and the GP energy functional in the rotating frame 
Eq. (|3.6p simplifies to 



Mft-rxv 

■ ■ 



(4.7) 



Like the TF limit E TF [$] in Eq. $Z7Ify for a stationary 
condensate, this equation depends only of \^f\ 2 and \^\ , 
but it now includes the hydrodynamic flow through v. 

In the typical case of rotation about z, simple manip- 
ulations yield the equivalent expression 



dV\M (v 
dV [(Vtr - 



Vsb) 



E 'v + ^TF i 



(4.8) 



where £Z£ and E' TF denote the contributions in the first 
and second lines. For a large condensate with many 
vortices, the actual superfluid velocity v closely approxi- 
mates the classical solid-body value v s ^, so that E' v [the 
term containing (v — v s ^) 2 ] vanishes to leading order. 
The remaining terms _E^ F (f2) have exactly the same form 
as Etf in Eq. (|2.24p . but the radial part of the con- 
fining potential is altered to \M{uj\ — Vt 2 )r 2 . Varia- 
tion with respect to |\I/| 2 gives an equation of the same 
form as Eq. (|2.25p , but with the reduced radial confining 
potential. Thus the TF density again has the form of 
Eq. (|2.27[) . but with modified TF condensate radii that 
depend on f2. The normalization condition in three di- 
mensions shows that 



Mtf(^) = /xtf(0)(1 - n 2 /wi) 2 / 5 , 



(4.9) 



so that the chemical potential decreases continuously and 
formally vanishes for Q — > luj_ , as does the central density 
fiTF (fi) / 9 (because of the reduced radial confinement) . A 
combination of these results leads to the specific SI depen- 
dence of the condensate dimensions dButts and Rokhsarl 
119991 : iFetterl . l200lt iRaman et all l200lh 



-3/10 



R z (0) \ - 2 



1/5 



R ± (n) _ / o 2 ■ 

(4.10) 

As anticipated, the external rotation expands the con- 
densate radially and shrinks it axially, approaching a two- 
dimensional configuration. The corresponding aspect ra- 
tio 



R Z (Q) _ ^uj'i - Q 2 



(4.11) 



provides a convenient diagnostic tool to in fer the actua l 
angular velocity of the rotating condensate ( Halian et all 
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FIG. 11 Side view of BECs in a trap, (a) Static BEC, with 
aspect ratio R z /R± = 1.57. Dramatically reduced aspect 
ratio is evident for (b) Q/u>± = 0.953 and (c) Q/uj± = 0.993. 
From lSchweikhard et all (|2004T ). 



l2001bt IRaman et all , l200lb ISchweikhard et all , I2004D . as 
seen in Fig. [TT1 

The two terms in Eq. (|4.8p have quite different magni- 
tudes in the mean-field TF regime. The dominant con- 
tribution E' T¥ (fi) is the TF energy of the rotating con- 
densate with the modified squared radial trap frequency 
u) 2 ^ — > iJ\_ — fl 2 . This quantity follows directly from 
Eqs. (j2~30l) and gU) 

^tf(^) - - 

The thermodynamic relation in Eq. (|4.4p yields the cor- 
responding TF angular momentum of the rotating con- 
densate, and some algebra gives the simple expression 



f n TF {n)Not (l-ft 2 A4) 2 / 5 . (4.12) 



Ltf = | MNR 2 ± (n) n = i sh (n) n, 



(4.13) 



where isb(^) = ^ MNRj_(Q.) is the solid-body moment 
of inertia for the deformed rotating condensate (this 
expected result follows from a separate calculation of 
(r 2 ) TF ). 

The remaining term E' v in Eq. (|4.8p gives a small non- 
classical contribution to the angular momentum of the ro- 
tating condensate. If the volume integral is approximated 
as a sum of integrals over circular cylindrical cells of ra- 
dius I centered at each vortex, a straightforward analysis 
yields 



E' v (n) w JVfifiln(Z/f). 



(4.14) 



With logarithmic accuracy (ignoring the f2 dependence 
inside the logarithm), the corresponding nonclassical an- 
gular momentum then becomes 



(4.15) 



Thus the total angular momentum of the rotating 
condensate indeed displays nonc lassical behavior [see 
Sees. 2.5 and 8.3 of lLeggettJ (|2006l )] 



I sh (n)n- Nhln(l/£) 



(4.16) 



Note that the correction reflects the inability of the su- 
perfluid to rotate like a solid body. It arises directly from 
the quantization of circulation and would vanish in the 
classical limit (h — > 0) when k — > and n v — 2Q/k — > oo. 
The correction L' v also becomes small in the limit of a 
dense vortex array (I — > £), as discussed in Sec. V.B.I. 
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B. Experimental and theoretical studies of vortex lattices 
in the TF regime 

The first images of one or a few vortices (Sec. III.D) 
rapidly led to the cre ation of large vortex ar rays. Ini- 
tially, the MIT group (|Abo-Shaeer et all l200ll ) use a ro- 
tating deformation (like the ENS technique — see Fig. |6|) 
to obtain remarkable images of triangular vortex lattices 
that are strikingly uniform, even near the outer edge of 
the condensate. These arrays closely resemble the clas- 
sic triangular flux-line la t tice in type-II superconductors 
predicted by lAbrikosovl ( 1957| ) and the triangular vor- 
tex lattice in rotating neutral in compressible superfluids 
predicted by iTkachenkol ( 1966al) . 



0>) 




FIG. 12 Images of expanded 87 Rb condensates showing (a) 
small vortex array for slow rotation and (b) large vortex array 
for rapid rota tion. Note the highly re gular triangular form. 
Adapted from lCoddington et al\ &004 ). 



The JILA group (|Engels et all 120031 : lHalian et all . 
l2001bD take a quite different approach to imparting an- 
gular momentum to the condensate. They start by 
spinning-up a normal cloud above T c with an elliptically 
deformed rotating disk-shaped trap. To cool into the 
BEC regime, they adiabatically alter their trap to cigar- 
shaped and then evaporate normal atoms near the ends 
of the axis of symmetry (these atoms have small angular 
momentum). This process yields a rotating condensate 
with relatively large angular momentum per particle and 
good triangular arrays, as seen in Fig. 1121 These rapidly 
rotating BECs facilitate a series of detailed experiments 
on the properties of the vortex lattice. 



1. Collective modes of rotating condensates 

The study of collective modes of a static condensate 
provides one of the crucial tests of the ti me-dependent 
GP eq uation, as discussed in Sec. IV of iDalfovo et al. 
(1999) and Chap. 12 of iPitaevskii and Stringari 



(2003). One particul arly effective theoretical ap- 
proach ( Stringari Il996l ) linearizes the hydrodynamic 
equations (|2.38[) and (|2.39[) in the small perturbations Sn 
and Sv = V(5$. Since this formalism explicitly assumes 
irrotational flow, it fails for a rapidly rotating TF 
condensate with a large vortex lattice, where the average 
s uperfluid velocity has uniform vorticity Vx»~ 2fi. 

ICozzini and Stringaril (|2003l ) propose a simple gener- 
alization, based on the Euler equation of nonviscous hy- 
drodynamics, retaining the full convective term dv/dt + 
(v ■ V)u. A familiar vector identity yields the equivalent 
form of the Euler equation in the laboratory frame 

M—+V(-Mv 2 + V tI +gn) = Mv X (V Xv), (4.17) 

which includes both the trap potential Vt r and the 
Hartree potential gn, but omits the quantum kinetic en- 
ergy from Eq. (|2.39j) . as appropriate in the TF limit. 
For a nonrotating irrotational condensate, V X v van- 
ishes, and Eq. (|4.17[) then omits the right-hand side, re- 
producing the gradient of Eq. (|2.39p in the TF limit. In 
contrast, a rotating TF c ondensate has many v o rtices , 
with V X v m 2fl. Thus ICozzini and Stringaril (|2003[) 
replace the right-hand side of Eq. (|4.17p with 2Mv X ft, 
which they call the limit of "diffused vorticity," although 
"mean-field" or coarse-grained" vorticity may be more 
descriptive. A detailed analysis gives the frequencies of 
the lowest m — ±2 quadrupole modes in the laboratory 
frame 



uj (m = ±2) = aJ: 



2uj 2 , -n 2 ±n. 



(4.18) 



In particular, the difference of the two frequencies is 20, 
whereas the sum of the squar ed frequencies is 4a; 2 , inde- 
pendent of the rotation rate. lHalian et all ( 2001bl ) fully 
confirm both these detailed predictions. 

As an extension of these experiments, lEngels et al\ 
(2002) study the dynamics of the (originally tri- 
angular) vortex lattice in a rapidly rotating con- 
densate. They can image in all three per- 
pendicular directions, which allows them t o con- 
firm theoretical predictions (Feder and Clarkl l200ll : 
iGarria-Ripoll and Perez-Garciall2001b ) that the vortices 
are nearly straight in such dense arrays. Furthermore, 
Eq. (|4.18p shows that the m — ±2 quadrupole modes 
have several i ntriguing fea t ures i n the limit Q — ► uj±. 
In particular, lEngels et all ( 20021 ) note that the result- 
ing m = —2 mode is small and hence nearly stationary 
in the laboratory frame, but the corresponding m = 2 
mode is nearly stationary in the rotating frame because 
uj(m = 2)' = uj(m = 2) — 2fi then becomes small [see 

In a particularly interesting experiment, lEngels et all 
(2002) form a rotating vortex lattice with an angular ve- 
locity fl/oj± ~ 0.95. They then apply a relatively small 
trap distortion with (uj 2 — lu 2 )/(oj 2 + uj 2 ) ~ 0.036 that is 
stationary in the laboratory. Such a fixed perturbation is 
nearly resonant with the slow m = —2 quadrupole mode. 
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As a result, the rotating condensate acquires a stationary 
quadrupole distortion that grows to ~ 0.4 over a time of 
~ 300 ms. Hence the rapidly rotating triangular vortex 
lattice experiences an oscillating shear distortion with a 
period that is g times the overall rotation period 2tt/£1 
because of the six-fold symmetry. Every time a basis vec- 
tor of the rotating lattice lines up with the minor axis of 
the fixed elliptical confining potential, the vortex lattice 
deforms significantly, forming what are effectively paral- 
lel stripes of closely spaced vortices. 

To understand this behavior in detail, recall that the 
stationary elliptical distortion of the rotating conden- 
sate induces a quadrupolar irrotational flow (see the 
discussion at the beginning of Sec. III.D.2) that com- 
bines with the overall rotation v s \, = X r to pro- 
duce the total flow velocity v. In the laboratory frame, 
the vortices follow stream li nes o f this total local veloc- 
ity. ICozzini and Stringaril (|2003l ) integrate the relevant 
dynamical equations of motion for the vortices. They in- 
deed find stripe-like patterns whenever one of the lattice 
vectors lies along the minor axis of the distorted conden- 
sate. 

In a related approach, iMueller and Hoi ( 20031 ) use a 
variational two-dimensional wave function. Although 
their formalism strictly applies only in the extreme limit 
of Q — > u>± when the interaction energy is small com- 
pared to the single-particle energy (the "lowest Landau- 
level" regime, discussed below in Sec. V), their trial 
wave function can describe a two-dimensional lattice of 
vortices with general symmetry, using Jacobian elliptic 
functions th at are closely relat ed to those in the ear- 
lier work of iTkachenkol (fl966al) . In their simulation, a 
small elliptical distortion of the trap potential at t = 
induces an increasing distortion of the condensate's as- 
pect ratio. The initial triangular vortex array remains 
a well-defined lattice but deforms following the com- 
bined velocity field of the rotational and irrotational 
flow. As time evolves, stripes appear periodically, pro- 
ducing im ages that clo s ely re semble those found experi- 
mentally (|Engels et aLl . l2002h . 



2. Uniformity of the vortex array 



Nearly 30 years ago, ICampbell and Zif3 ( 1979f ) stud- 
ied numerically the various equilibrium arrangements of 
two-dimensional vortex arrays in uniform incompressible 
ideal fluids inside a fixed circular boundary (as a model 
for superfluid 4 He in a rotating bucket). They found pre- 
dominantly distorted triangular arrangements with con- 
centric circl es of vortices and a dep leted region near the 
outer edge. lYarmchuk et all ( 1979T ) see similar patterns 
in rotating superfluid 4 He (see Fig. [8J, although the small 
number of vortices precludes a definitive comparison. In 
contrast, images of vortex lattices in rotating BECs are 
strikingly triangular and regular, even out to the edge of 
the visible condensate, as seen in Fig. [T2j 

The treatment of the density profile in the mean-field 



Thomas-Fermi regime omits the term (v — v s ^) 2 entirely 
[see Eq. (|4.8[) ]. which assumes a uniform vortex distri- 
bution a nd ignores the fine-grain s t ructure of individual 
vortices. ISheehv and Radzihovskvl ( 2004al lbl) investigate 
these effects by considering a two-dimensional TF con- 
densate with number density n(r) = n(0) (1 — r 2 /i?^J, 
where R± grows with increasing f2. This rotating con- 
densate contains an array of vortices at regular two- 
dimensional lattice sites Tj with mean density n v = 
Mfl/(nh) [the Feynman value in Eq. P~Tj) ]. They then 
subject each vortex to a small displacement field u(r), so 
that the original site moves to rj + u(rj), and the vortex 
density becomes 



n v (r) =n v [1 - V • u(r)} 



(4.19) 



This distorted vortex density in turn induces a change in 
the mean background fluid flow based on the following 
integral 



v(r) « / d 2 r'n v (r') Z * ^ ^ = ftxr-2ft Xu(r), 



h 

M J ' ' \r — r' 

(4.20) 

where $1 = flz and the second term is a correction to the 
solid-body value. 

A detailed analysis with Eqs. ([4~8]) . (j4~19l) and (|4~20|) 
expresses the deformation part of the energy E' in terms 
of it, and the corresponding Euler-Lagrange equation 
leads to the general expression 



I 2 fl 2 ' 
u ( r ) ~ -g- in ( ^ ) Vlnn(r) 



(4.21) 



where n(r) is the condensate number density, / is the 
mean circular cell radius from Eq. (|4. 5[) . and £ is eval- 
uated with the central density. 8 In the particular case 
of the two-dimensional TF density profile with n(r) oc 
1 — r 2 /R\, this equation yields 



u(r) 



I 2 



1 -r 2 /R\' 



(4.22) 



Note that the deformation of the regular vortex lattice 
is purely radial, as might be expected from the circu- 
lar symmetry, and it increases with increasing distance 
from the center. Furthermore, the quantity R\jl 2 is the 
number M v of vortices in the rotating TF condensate [see 
Eq. ()4.6|) ] , so that the nonuniform distortion is of order 
M~ x (at most a few %), even though the TF density 
decreases significantly near the edge of the condensate. 



8 As a very simple application, assume a uniform particle den- 
sity, so that u vanishes. Equation 114.191 1 then indicates that the 
vortex density is also uniform. Note that this formalism relies 
on the TF energy functional and has little connection to the 
nu merical studies of vo r tex a rrays in uniform incompressible flu- 
ids l|Campbell and Zifil. Il979h . where a vortex-free region occurs 
near the outer wall. 
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FIG. 13 Lattice spacing as a function of the scaled radial 
position r/R±. The solid curve is the prediction (|4,24[) from 
ISheehv and Radzihovsky (2004a b), with the rigid-body value 
as a dashed line. Plots (c) and (d) show that the fractional 
amplitude decreases with increasing Q, and (e) restores the 
suppressed zero in (c), indicating that the overall effect is 
small but measurable. The intervortex spacing increases by 
less than 2% ove r a region where the atom ic density changes 
by ~ 35%. From lCoddington et al\ (120041 1. 



Finally, Eq. (|4.19p gives the corresponding nonuniform 
axisymmetric vortex density 



1 



n v {r) = n v + — In I —z V lnn(r) 



8tt 



I 2 



which becomes 



n v (r) 



2txR\ 



In 



(1 -r 2 /R 2 ± 



2 ^2 ' 



(4.23) 



(4.24) 



for the two-dimensional TF co ndensate. The cor r ection 
is again small, of order N^ 1 . ICoddington et al. I (12004TI 
study this behavior in detail, as shown in Fig. 1131 Com- 
parison of images in Fig. [TBT c) and [TBT d) show that the 
distortion decreases with increasing number of vortices 
(namely with increasing fl), and Fig. ITBT e) restores the 
suppressed zero in (c), indicating that the overall effect 
is indeed small. 



3. Vortex core size for large rotation speeds 

The principal effect of rapid rotation is the altered as- 
pect ratio seen in Fig. [TTJ but the rotation also affects 
the vortex core radius £. In the TF limit of a large con- 
densate, the basic definition in Eq. (|2.3ip shows that £ 2 



increases with increasing Q because the radial expansion 
reduces the central density n(0) and thus the chemical 
potential /i = gn(0). To quantify this behavior, it is help- 
ful to consider the dimensionless ratio £, 2 /l 2 — Ml/(2fi), 
which is roughly the fraction of each vortex cell that 
the core occupies. A combination of Eqs. (12.281) . ()2.29D . 
(123T|) . and flO} yields the TF expression (jFetterl . 2001) 



e 



Ml 
2^ 



it 



(1 - 57 2 ) 2 /5 \ uj z i5Na J 



2/5 



(4.25) 



where H = H/lu± is a scaled dimensionless rotation speed. 

Since d±/(Na) is small for a typical TF condensate, 
the ratio £ 2 /I 2 also remains small until SI approaches one, 
and it t hen grows rapidly becaus e of the small denomina- 
tor. As lFischer and Bav m (2003) emphasize, this limit of 
rapid rotation requires a more careful treatment. Specif- 
ically, they do not assume that £ 2 is proportional to jti -1 . 
Instead, they treat £ 2 as a variational parameter and find 
that the ratio £, 2 /l 2 grows only until it reaches w i, when 
the vortex core occupies a significant fraction of each unit 
cell in the vortex lattice. Strictly, the TF limit no longer 
applies for this large fi < 1 because the density variation 
now m akes a significant contribut ion to the total kinetic 
energy. iBavm and Pethickl ( 2004[ ) construct a more gen- 
eral variational theory that includes both the TF limit 
and the lowest-Landau-level regime described below in 
Sec. V. 
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FIG. 14 Fractional vortex core area as function of 2hQ/fi. 
The dashed line is the TF theory and the dotted line is 
the l owest-Landau-level theor y discussed below in Sec. V. 
From lCoddington et all (|2004r i. 



ICoddington et al] ( 20041 ) study the size of the vortex 
core from both experimental and numerical perspectives. 
Because of different fits to the vortex-core density profiles 
and different definitions of the fractional core area, their 
numer ical values differ somewhat fro m those in Eq . (14.251) 
and inlFischer and Bavm (2003) and lBavm and Pethickl 
(2004). Figure [14] shows typical experimental and nu- 
merical values of the fractional core area compared with 
the ratio 2Ml/fi. For small values of f2, Eq. (|4.25[) in- 
dicates that £ 2 // 2 should vary linearly with Q,/fx in the 
TF regime. As f2 increases toward the critical value u±, 
however, the fractional core area sa turates at a constant 
value. Detailed numerical studies ( Cozzini et all |2006[ ) 
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of a rapidly rotating condensate yield a very similar pic- 
ture of this behavior. 



4. Tkachenko oscillations of the vortex lattice 



In 1966, iTkachenk o (1966a.b) published two remark- 
able papers on the behavior of arrays of straight vor- 
tices in an unbounded rotating incompressible irrota- 
tional fluid (as a model for superfluid 4 He). The first 
paper shows that a triangular lattice has the lowest en- 
ergy of all simple lattices (those with one vortex per unit 
cell). The second paper studies small perturbations of 
a general lattice, showing that the square lattice is un- 
stable for waves in certain directions, but that the tri- 
angular lattice is stable for all normal modes. For the 
latter structure, he also determines the dispersion rela- 
tion of the small-amplitude normal modes of a vortex 
lattice, along with the corresponding eigenvectors for a 
given wave vector k lying in the xy plane perpendicular 
to fi. The calculation is a tour de force of analytic func- 
tion theory, involving several elliptic-type functions, es- 
pecially the Weierstr ass C function that has a s i mple pole 
at ea ch lattice site ( Chandrasekharanl Il985t IWolframl 
120071 ). Ultimately, however, the final result is simple. At 
long wavelengths fci < 1, where I — y/fr/(MQ) is the 
vortex-cell radius, the wave is predominantly transverse 
with a linear dispersion relation 



uj(k) « crk, 



where 



ct — 




(4.26) 
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is the speed of propagation. This long-wavelength motion 
is effectively a transverse phonon in the vortex lattice. 9 
The significant radial expansion of rapidly rotating 
BECs means that the resulting vortices are essentially 
two-dimensional. Thus bending modes are irrelevant, 
but the nonzero compressibility of th ese atomic gases r e- 
quires a si gnifican t modi ficatio n to Tkachenkol ( 1966bf )'s 
analysis. ISoninl (|l987f) and iBavml (|2003f T generalize 
Tkachenko's result in Eq. (14. 26[) to find the new long- 
wavelength expression 



uj{kf 



„2 7 2 



„2j.2 



4Sl 2 + s 2 fc 2: 



(4.28) 



where s is the speed of sound. This express i ons as - 
sumes an infinite unifo rm system, b ut lAnglinl (|2003 ) ; 
iBaksmatv etai\ (|2004h ; IBavml (|2003l ); ISoninl (l2005l Talso 
include the nonuniformity of the trapped BEC. If sk ^ SI 



9 Unlike Newtonian particles in a two-dimensional lattice, however, 
vortices obey first-order dynamical differential equations. Hence, 
there is only one normal mode for each wave vector fc. 



(the short- wavelength or incompressible limit) , then this 
expression reduces to Tkachenko's uj{k) ss orfc, but if 
sk <C SI (the long- wavelength or compressible limit), 
then the mode becomes soft with uj(k) oc fc 2 . Simi- 
lar softening of the collective-mode spectrum has impor- 
tant con sequen c es for stability of the vo rtex lattice at 
large SI (|Bavml 120051: ISinova et all [2002), as discussed 
in Sec. VI. 




FIG. 15 Lowest Tkachenko mode of the vortex lattice ex- 
cited by atom removal taken (A) 500 ms after the end of 
the blasting pulse and (B) 1650 ms after the blasting pulse. 
Lines are sine fits to the di stortion of the vortex lattice. 
From ICoddington et all (|2003ft . 



Coddin gton et al. I (|2003l) observe these Tkachenko os- 
cillations in considerable detail. They start with an es- 
sentially uniform vortex lattice and then apply a weak 
perturbation (along with other approaches, they use a 
"blasting pulse" that removes a small fraction of the 
atoms near the center of the trap). At a scaled rota- 
tion frequency SI = 0.95, they wait a variable time and 
then turn off the trap, yielding an image of the result- 
ing distorted vortex lattice. Figure [15] shows the vortex 
lattice at one quarter and three quarters of the oscilla- 
tion period, with the lines as a sine fit to the distor- 
tion of the vortex lattice. In practice, the normal mode 
has the correct shape, but the measured frequency is 
significantly smaller than that predicted in Eq. (|4.28D . 
lAnglin and Crescimannol (2002) suggest that this dis- 
crepancy arises from reliance on a continuum theory that 
ig nores the increase d vortex core radius for large SI. 

ITkachenko! (|l969h also evaluates the shear modulus of 
the triangular vortex lattice in an incompressible fluid, 
which is directly related to th e long-wavelength o scilla- 
tion spectrum. Subsequently, ICozzini et al. (2006) use 
numerical methods to study the more general situation 
of a vortex lattice in a rotating BEC, which also includes 
both the inc ompressible (TF) limit and the c ompressible 
(LLP limit ()Sinova et all l2002HSoninl . l2005h . 



5. Rotating two-component BECs 

The experimental creation of relatively large 
vorte x arrays in s ingle- components BECs in 
23 Na (lAbo-Shaeer et all l200lh and 87 Rb (lHalian et all 
l2001bl : lHodbv et a/.l . 12002 ; iMadison et aLl . l2000af, rapidly 



led to p roposals for corresponding two-component vorte x 
arrays (jKasamatsu et all 120031 : Mueller and~Hcl l2002h . 
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As discussed in Sec. III.D.l, the interspecies interac- 
tion parameter g 12 significantly affects the behavior, 
especially for trapped condensates (the intraspecies 
p arameters g\ and gi are assumed positive). 

iMueller and Hoi ( 2002f ) initiated the theory of these 
two-component condensates, assuming rapid rotation 
(Q < u)±). This limit allows a considerable simplification 
(see Sec. V) because the condensate wave functions "Ji 
and ^ 2 can then be constructed from the lowest-Landau- 
lev el (LLL) single-part i cle sta tes [the same as those used 
bv iButts and Rokhsarl (|l999h in Sec. HI.C]. Apart from 
an overall Gaussian factor, the resulting trial solutions 
are analytic functions of the complex variable £ = x + iy; 
they involve quasiperiodic Jacobi theta functions that 
are closely related to other elliptic functions appearing 
in the treatment of normal modes in a rotating v ortex 
lattice ( Chandrasekhararl Il985t iTkachenkol . Il966al lbl). as 
mentioned in Sec. IV. B. 4. This analytical approach al- 
lows a relatively simple treatment of the change in the 
vortex lattice as the parameter a = 9x2/ \Jg\92 varies 
from negative (when the interspecies interaction is at- 
tractive) to positive. Specifically, for negative a, the two 
components prefer to overlap and form a single triangular 
lattice. As a becomes increasingly positive, however, the 
overall repulsive interaction energy leads the two com- 
ponents to separate spatially, with interlaced triangular 
vortex arrays to minimize the overall density variation. 
Eventually, the triangular lattices should distort, forming 
square or rectangular arrays when a ~ 1. 

Subsequently, iKasamatsu et all (2003) study the be- 
havior at lower angular velocity in the mean-field TF 
regime. Their extensive numerical analysis of the equi- 
librium vortex-lattice configuration as a function of the 
scaled dimensionless rotation rate f2 = Q/u>± and the 
relative interaction strength a gives a qualitatively sim- 
ilar picture: For any reasonable f2, the interpenetrating 
vortex arrays assume four-fold square symmetry as a ap- 
proaches 1~ from below. 

ISchweikhard et all ( 2004h perform a series of experi- 
ments on these rotating two-component systems (Cl ss 
0.75). They start with a single-component triangular vor- 
tex array in the state |1) = \F = l,mp = —I) (compare 
Sec. III.D.l). A short pulse with the two-photon elec- 
tromagnetic coupling fields transfers ~ 80% of the pop- 
ulation to the state |2) = \F = 2, nip = I). The initial 
triangular vortex lattice is now unstable, and after a tur- 
bulent stage (~ 2 s), it transforms to a square structure 
that persists for another ~ 4 s. Ultimately, however, the 
state 1 2) decays through a spin-flip transition. As a re- 
sult, the remaining single-component condensate in state 
|1) again transforms back to a triangular array, similar 
to the initial state. 



6. Rotating gas of paired fermions in tightly bound (BEC) limit 



fermionic atoms such as 6 Li and 40 K. Although a de- 
tailed discussion is inappropri ate here, Sees. I and V III of 



iBloch et all (l2007h . Sec. Ill of lGiorgini efaU (120071) . and 

3tt ( 



Sec. 8.4 of lLeggetd (|2006h provide excellent treatments of 
the behavior near a Feshbach resonance. Here, I simply 
review the basic elements of scattering theory to explain 
how the s-wave scattering length can exhibit a resonance 
associated with the appearance of a new bound state. 

The simplest description of a dilute atomic BEC as- 
sumes a weak repulsive effective interaction g > with 
a positive s-wave scattering length a > (this picture 
holds for the bosonic atoms 87 Rb and 23 Na). Such a sys- 
tem is stable both in bulk and in a trap, with a variety 
of associated collective modes. In contrast, the bosonic 
atom 7 Li has an attractive effective interaction g < 
with a negative a < 0. In bulk, the negative interac- 
tion would yield an imaginary speed of sound, implying a 
long-wavelength instability. In a trap, however, the pos- 
itive kinetic and potential energy can stabilize the con- 
densa te for sufficient l y small N, but the sy s tem w ill col- 
lapse ( Bradley et all , Tl997l : [Ruprccht et all , Il995h above 



a critical N c . As an alternative picture of this stabi- 
lization, note that the quantized wavenumbers in a trap 
provide a lower bound k mm ~ tt/R, where R depends on 
N. If the discrete analog of the Bogoliubov energy spec- 
trum is real for all k > /c m i n , then the trapped condensate 
is stable even for negative a. 

The direct relation between the sign of g and the sign 
of a holds only for weak potentials that have no bound 
states. The situation becomes more interesting when the 
potential can support one (or more) s-wave bound states. 
In this case, the s-wave scattering length exhibits diver- 
gent resonant behavior. To be very specific, let us con- 
sider a spherical square- well potential V(r) = Vo6(b — r) 
with height Vo and range 6, where 9 denotes the unit 
positive step function. 

1. Bound states occur only for negative Vb = — \Vo\. 
A standard analysis [see, f or example, Sec. 33 of 
[Landau and L ifshitz (|l965f) and pp. 164 and 349 
of IShankarl (jl994l )] shows that there are no s-wave 
bound states if \V \ < n 2 h 2 / (8Mb 2 ); as \V \ in- 
creases, the nth s-wave bound state appears when 



(2n-l) 2 TT 2 h 2 
8Mb 2 



\Vo\ 



(4.29) 



2. Consider now the scattering of a particle with pos- 
itive energy E = h 2 k 2 /(2M). The radial I = 
wave function tp(r) = u(r)/r has the asymp- 
totic form u(r) = sin(fcr + 8), where 5(k) is the 
s-wave phase shift. A detailed analysis for the 
same squa re-well potential V(r) [see , for example 
Sec. 130 of Landau and Lifshitz! (|!965h and pp. 549- 
554 o flSJiankad |1994)] yields the explicit expression 



One of the many remarkable new results from the study 
of ultracold gases is the superfluidity of bound pairs of 



5 = -kb- 



aril an ( — tan k'b 



(4.30) 
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where k' 2 = 2M(E — V Q )/H 2 determines the 
wavenumber k' inside the square well. 10 It is ev- 
ident physically that S is negative for positive Vo 
because the repulsive potential pushes out the wave 
function, whereas 5 is positive for negative Vo be- 
cause the attractive potential pulls in the wave 
function. This phase shift completely characterizes 
the s-wave scattering by this potential. 

It is conventional and convenient to introduce the 
s-wave scattering length a = — lim^o 8{k)/k] for 
the square-well potential, this result gives 



Magnetic Field [G| 

vn 



tan k'b 
~~k^ 1 



(4.31) 



where k' 2 — —2MVo/h 2 is now evaluated for k = 
(and hence E — 0). 11 If Vo is positive and large, 
then a — > b, but the dependence can be very differ- 
ent in other situations. For example, a vanishes if 
Vq = (as expected). More significant is the behav- 
ior for negative Vo, especially near the appearance 
of a bound state. Equation (|4.29|) shows that k'b 
increases through (n — ^)tt when the nth s-wave 
bound state first forms. Consequently, Eq. (|4.3ip 
diverges to negative infinity when the state is just 
unbound and then decreases from positive infinity 
when the trapped state is weakly bound. Such res- 
onant dependence is a typical and generic feature 
of scattering by an attractive potential that can 
support one or more bound states. 

In the context of dilute ultracold gases, the attrac- 
tive long-range van der Waals potential acts to bind 
the two fermionic atoms to form a molecule. In gen- 
eral, there are several molecular states, each with differ- 
ent magnetic moments. Thus an applied magnetic field 
can shift a bound state in a closed channel across the 
asymptotic continuum in an open chann el, leading to a 
"Feshbach resonance" [see, for exam ple, Moerdiik et all 
(ll995lh lKohler et all (120061). Se c. 2 of ICastinl (l2006h and 
Sees. I and VIII of iBloch et all (|2007h for more detailed 
discussions]. In essence, sweeping an external magnetic 
field through a Feshbach resonance transforms a tightly 
bound bosonic molecule of the two fermions into a weakly 
attracting fermionic pair. The bosonic molecules can 
form a BEC that displays essentially all the proper- 
ties of an atomic BEC, whereas the weakly attracting 
fermionic pairs are analogous to the Cooper pairs in the 
BCS the ory of conventional low-temperat ure supercon- 
ductors (de Genned . [1961 iTinkha m. 1996). As a result, 
this transition through a Feshbach resonance is called the 
"BEC-BCS crossover." 
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FIG. 16 Vortex lattice in a rotating gas of Li fermions. The 
left figure is on the BEC side, the central figure is on reso- 
nance, and the right figure is on the BCS side [note that the 
horizontal scale is — (fcFa) -1 , which i s effectively — (n 1 ^ 3 a)~ 1 
apart from a numerical factor]. From [Zwierleinl (|2007l). 



Particula rly relevant in t h e pre sent context is the ex- 
periment of IZwierlein et all ( 20051 ) who rotate a BEC of 
tightly bound bosonic molecules with a deformed opti- 
cal trap (the magnetically induced Feshbach resonance 
precludes the use of the usual magnetic trap). They 
then sweep the magnetic field from the positive scatter- 
ing length (the BEC side of bound molecules) through 
the resonance to the negative scattering length (the BCS 
side of attracting fermionic pairs) . Figure [16] shows the 
vortex lattice at three different values of the magnetic 
field. Note that the scattering length a diverges in the 
resonance region. Consequently, it is conventional to 
use (kpa)~ l as the independent variable, which varies 
smoothly from positive in the BEC limit through zero 
near the crossover region to negative in the BCS limit. 
Here, fcj? is the Fermi wavenumber, but, apart from a 
numerical factor of order unity, it is the same as n 1 / 3 , 
with n the boson density. The explicit persistence of 
quantized vortices through the BEC-BCS crossover pro- 
vides strong evidence for superfluidity on the fermionic 
side, although the experimen t did not reach the w eakly 
bound BCS limit [Sec. VIII of lGiorgini et all (j2007f ) gives 
a detailed discussion of rotating ultracold Fermi gases] . 



10 If Vo > E, then k' is written as in and tan/c'fe/fc' becomes 
tanh Kb/a. 

11 If Vq is positive, then this relation is replaced by ft 2 = 2MVo/h 2 ; 
correspondingly, Eq. 114.311 1 becomes a = b — tanh Kb/ ' k. 



ISensarma et d\ (120061) use the fermionic Bogoliubov- 
de Gennes equa tions [see, for e xample, Chap. 5 of 
Ide Gennesl (|l966h and Sec. 10.1 of iTinkhaml (|l996f) ] to 
study the structure of a single quantized vortex in a uni- 
form three-dimensional atomic Fermi gas at zero tem- 
perature. Here I consider only the BEC side, since a full 
summary of ultracold Fermi gases is beyond the scope of 
this review. 

Let us focus on the vortex core radius (healing length) 
£ defined in Eq. (|2.31[) for a trapped condensate. As 
noted in Sec. IV. B. 3, £ grows with increasing rotation 
speed because of the reduced central density [and hence 
reduced chemical potential /i(fi)], but Eq. (|4.25p also ex- 
hibits the explicit dependence on the scattering length 
a. Here, it is preferable to consider the slightly different 
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dimensionless ratio 



d 



(1 



J ( d \ 

n 2 y/s \l5Na) 



1/5 



(4.32) 



which involves only the scaled rotation speed f2 = Q/u/±, 
the s-wave scattering length a, and the trap geometry 
t hrough dp and the rad ial rap frequency luj_. 

ISensarma et all (2006) remark that the healing length 
in a uniform cold BEC Eq. (|2.21[) is proportional to 
(na) -1 / 2 and decreases with increasing positive a for 
fixed n. Consequently, the vortex core size should shrink 
on approaching the resonance regime from the BEC side. 
Note, however, that inclusion of the trap in Eq. (|4.32[) 
alters the fractional power from a -1 / 2 to a -1 / 5 . Thus 
the dependence on a may be difficult to measure in cur- 
rent experiments. It is also clear from Eq. (|2.22[) that the 
GP picture fails as £ approaches the interparticle spac- 
ing, which is expected to be the minimum value for £. 
On the BCS side, the qualitative relation £ ~ hvp/A 
shows that the vortex core again grows on moving away 
from the resonance regi me because the BCS e nergy gap 
A then be comes small (|Sensarma et all , l2006t ) [see also 
Sec. 8.4 of lLeeeetd (12006b ]. 

On the BEC side, recall that the condensate density 
n(r) around a vortex decreases rapidly from its bulk value 
n for r < £ inside the vortex core. Consequently, the 
circulating mass current j(r) = Mn(r)v(r) = n(r)h/r 
reaches a maximum "critical" value j c ~ nh/£ for r ~ £ 
[compare the discussion below Eq. (|3 .4(1 ] . As a increases 
on the BEC side toward the Feshbach resonance, the de- 
creased £ implies that this critical mass current j c grows 
and sa turates when $ ~ tt," 1 / 3 . A more detailed anal- 
ysis by ISensarma et all ( 20061 ) shows that j c then again 
decreases on moving from the resonance regime to the 
BCS side. Thus both £ and j c display non-monotonic 
dependence in traversing the BEC-BCS crossover region. 



C. Effect of additional quartic confinement 

The discussion of rotating condensates in the mean- 
field Thomas-Fermi regime (Sec. IV. A) makes it clear 
the a purely harmonic trap with radial frequency tu±_ will 
cease to confine the system if the rotation speed f2 reaches 
or exceeds u>±. Indeed, the limit fl — » u>± is singular, 
because the condensate radius and the total angular mo- 
mentum both diverge. In the regime lo±_ — Q <C u>±, the 
system is expected to undergo a series of crossovers and 
complicated phase transitions (Sees. V and VI) whose 
exact nature remains under debate. 

Thus it is interesting to consider the effect of adding 
an additional strongly confining potential that rises more 
rapidly than r 2 . I here focus on the most common 



(quartic) model oc r 4 dFetterl. l200lt iKasamatsu et all. 



l2002t iKavoulakis and Bavml . 120031 iLundhl . I2002T) . Such 
an anharmonic trap will confine the condensate even for 
51 > u)_\_, thus allowing a more controlled study of possible 



new states (they are not generally the same as those pre- 
dicted in a harmonic trap in the limit Q — > to±). In terms 
of the usual dimensional quantities, the anharmonic ra- 
dial trap potential has the form 



V tr (r) 



hut i 



r 

7W 



+ \ 



rf 4 



(4.33) 



where dj_ = \J h/(Muj±) is the radial oscillator length 
and A is a dimensionless parameter that characterizes 
the admixture of quartic contribution. 

For simplicity, I focus on the limit of rapid ro- 
tation and consider only a two-dimensional conden- 
sat e that is uniform in the z direction over a length 
Z ([Fetter et all . 120051 : IKavoulakis . and Bavml . l2003h . with 
^(r, z) = ip(r)/y/Z. It is also convenient to measure en- 
ergies in units of fkij± and lengths in units of d±, and 
to normalize the dimensionless condensate wave function 
ifi with Jd 2 r\4>\ 2 = N. In the mean-field TF regime, 
the dimensionless energy functional in the rotating frame 
[compare Eq. (|4.7p ] becomes 



E'[tP] = J d 2 r [i (v 2 +r 2 + Ar 4 ) |?/>| 2 



(4.34) 



where g = Aira/Z is a dimensionless two-dimensional in- 
teraction parameter. An alternative model for the axial 
confinement uses a tight harmonic trap, in which case 



the interaction parameter is g = V87ra/d z , but the sub- 
sequent description is unaffected. 

When the rotating condensate has many vortices, the 
total velocity v closely approximates the solid-body form 
f2 X r, in which case Eq. (|4.34p reduces to 

£'[</>] = J d 2 r{i[(l-fi 2 )r 2 +Ar 4 ] |Vf + \g W) . 

(4.35) 

As in Eq. (|2.25[) . variation with respect to |-0| 2 at fixed 
normalization yields two-dimensional density 



5|<Kr)| 2 =/i+i[(fl 2 -l) r 2 -Ar 4 



(4.36) 



For Q < 1, the density has a local maximum at the center, 
but it changes to a local minimum for £1 > 1. In addition, 
the central density [■(/'(O)l 2 is proportional to the chemical 
potential [i. As in the previous Eq. (|4.9|1 . [i depends on 
f2 through the normalization and decreases continuously 
with increasing f2, reflecting the radial expansion of the 
condensate. 

In contrast to the mean-field TF picture for a har- 
monic radial trap, however, the chemical potential here 
can vanish at a critical angular velocity flh > 1, when the 
central density also vanishes. This behavior indicates the 
formation of a central hole; for fl > flh, the chemical po- 
tential fi becomes neg ative and the cond ensate assumes 
an annular form [see ( Fetter et al. I. l2005h for a detailed 
description of this unusual structure] . 

As n continues to increase, the system can, in prin- 
ciple, make a second transition to a different annular 
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state with purely irrotational (vortex-free) circulating 
flow, known as a "giant vortex." In this case, the cir- 
culating velocity v(r) has the quantized form Af v /r, re- 
flecting the presence of M v phantom vortices in the hole 
(equivalently, the flow simply represents A/",, -fold quan- 
tized circulation in a multip ly connected condensate). 
IKavoulakis and Bavml ( 20031 ) suggest the possibility of 



such irrotat ional flow in a quadratic p lus quartic poten- 
tial [earlier, iFischer and Bavml ([2003) propose a similar 
structure in a rotating rigid cylinder]. As noted at the 
start of Sec. IV, solid-body rotation [the model used in 
obtaining Eq. (|4.36[) ] gives the absolute minimum energy 
for a rotating system. Consequently, any irrotational gi- 
ant vortex necessarily depends on the finite quantization 
of circulation that introduces grainincss in the superflow. 
In practice, estimates of the critical angular velocity for 

formation of a giant vortex require considerable numeri- 
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cal work [see (Fetter et ai, 2005; Fu and Zaremba, 2006; 
IKavoulakis and Bavml . 120031 : iKim and Fetterl . 120051 ) for 
more details]. 

In a n elegant series of e xperiments, the E NS (Paris) 
group (|Bretin et all 12004 IStock et all \2004) study the 
effect of a weak quartic confinement in ad dition to the 
usual quadratic trap potential [see also ( Stock et all 
2005a) for a more general review of rotating BECs]. 
Specifically, they use a blue-detuned (repulsive) Gaus- 
sian laser beam propagating along the symmetry axis of 
their elongated condensate, with a resulting dipole po- 
tential 
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(4.37) 



The positive constant term Vq simply shifts to zero of en- 
ergy, the quadratic term reduces the original radial trap 
frequency, and the quartic term provides weak positive 
confinement. The measured oscillation frequency of the 
center of mass of the condensate give s the effectiv e radia l 
trap frequency uj c s/2it = 64.8 Hz (jStock et all 120041) . 
and the other trap parameters yield the quartic coupling 
constant A « 10~ 3 , indicating that the added confine- 
ment in Eq. (|4.33p is indeed weak. 
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FIG. 17 Pictures of the rotating gas taken along the rotation 
axis. The number below each picture is the rotation frequency 
in Hz (note that u c s /2n = 64.8 Hz is t he eff ective radial trap 
frequency). Taken from l|Bretin et aUl2004h . 

For the dimensionless rotation speed U w 1, the 
rotational deformation produces a nearly spherical con- 
densa te that remains stable up to SI < 1.05 ( Bretin et all 
l2004h . Figure [17] shows the expanded cloud for various 
rotation speeds both below and above the critical value 
uj c ff. Initially, a regular vortex array appears, but it 
becomes irregular in Figs. [T7TdV(f). In addition, for 



fl > w e fi, the condensate exhibits a clear local mini- 
mum in the density near the center, as expected from 
Eq. (|4.36p in the mean-field TF picture. Throughout the 
analysis, they fit the measured condensate shape to the 
TF prediction, which provides a direct determination 
of SI. They also use s urface-wave spectroscopy as 
discussed in Sec. III. D. 3 dCozzini and Stringaril 12003 ; 
Svidzinskv and Fetterl . 1 1991 IZambelli and Stringaril . 



19981 ) to obtain an independent measure of Q (even 
though the number of visible vortices appears to be too 
small for SI > Lo e s). 

The absence of visible vortices in Fig. [T7T g) is puz- 
zling. One possible reason for the discrepancy is the 
three-dimensional character of the condensate, which al- 
lows the vortices to bend. Also, the low temperature may 
make it difficult for the system to eq uilibrate. Numerical 
studies ( Aftalion and Danailal . I20CM) suggest that much 
longer times are needed to attain a well-ordered vortex 
lattice in these regimes of large SI. In any case, Fig. [T7T h) 
implies that the condensate simply collapses, although 
the reason remains unknown. 



V. VORTEX ARRAYS IN MEAN-FIELD 
LOWEST-LANDAU-LEVEL (LLL) REGIME 

In the mean-field Thomas- Fermi description of a rotat- 
ing BEC, the basic approximation is the neglect of the 
"kinetic energy" associated with the density variations. 
Specifically, Eq. |g2J omits the term h 2 (V\^\) 2 / (2M) 
whereas the usual kinetic energy from the phase varia- 
tion iMu 2 | , F| 2 remains. As seen in Eq. (|4.25|) , the vortex 
cores are small for moderate values of the dimensionless 
rotation speed S7 = tt/uj±, so that this approximation is 
valid. As SI increases toward 1, however, it fails when 
1 — SI ~ d±/ (15Na) <C 1 because the vortex cores then 
become comparable with the intervortex separation ~ 21. 
In this case, it becomes essential to return to the full GP 
energy functional E'[^f] in Eq. (|3.6p . 

In thinking about this new regime, IHoI (|200lh ob- 
serves that the interaction energy per particle [of order 
/i(Sl) = gn(0)] is small because the condensate expands 
radially [see Eq. (|4.9[) ] . In essence, the condensate be- 
comes two-dimensional, and it is convenient to study a 
two-dimensional condensate that is uniform in the z di- 
rection over a length Z . Thus the original condensate 
wave function ^(r, z) can be written as ip(r)/\/Z, where 
ip(f) is a two-dimensional wave function with normal- 
ization J d 2 r\ip\ 2 — N. The GP energy functional then 
becomes 



d 2 r ijj* 



\2M 



1 



Mlo 2 , r 2 - VLL 7 



+ ^2d\4>\ 2 )4>, (5.1) 



where g2D = g/Z- As in Sec. IV. C, one can instead 
assume a tight Gaussian axial confinement, in which case 
.92D = g/(V2nd x ), with d z = y/h/{Mu z ). 
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A. Physics of lowest-Landau-level (LLL) one-body states 



The one-body Hamiltonian in Eq. (|5.ip has the crucial 
advantage that it is exactly soluble with very si mple en 



ergyeigenvalues [see, for exampl e, Sec. 12.14 oflMessiab 
1961) and pp. 727-741 of iCohen-Tannoudii et al 



19771 )]. Consider the (nonrotating) two-dimensional 



isotropic harmonic oscillator with Hamiltonian 



P i 1 71 /r 2 2 



2M 



(5.2) 



Classically, the motion is that of two independent oscil- 
lators (along x and y) that can combine to give, for ex- 
ample, circular motion with frequency u>± in the positive 
(counterclockwise) or negative (clockwise) sense, depend- 
ing on the relative phases. When viewed from a frame 
that rotates with angular velocity fi in the positive sense, 
the two circular orbits with originally degenerate frequen- 
cies now have different frequencies 



u>± = lu± =F fi, 



(5.3) 



since the positive (negative) mode has a reduced (in- 
creased) frequency in the rotating frame. 

The quantum-mechanical description relies on the fa- 
miliar creation and annihilation operators 



1 / x 



. Pxdi 



al = 



1 ( x p x d A 
~~ y/2 \d± h 



(5.4) 

along with similar operators a y and al. In the present 
case, it is more appropriate to use circularly polarized 
states with operators 



a± = 



a x ^f i a 



y 



V2 



t 4 ± i a l 

V2 



(5.5) 



that obey the usual commutation relations [a_|_,at_] = 

[o_,oL] = 1, with all other commutato rs vanish i ng. A 
straightforward ana lysis [Sec. 12.14 of Messiah! ( 196lh 
and pp. 727-741 of ICohen-Tannoudii et all (|1977T) ] ex- 
presses Hq in terms of these new operators 



Hq = fkv_i 



1 



(5.6) 



Similarly the angular momentum L z = xp y — yp x be- 
comes 



(5.7) 



The number operators a±a± have non- negative integer 
eigenvalues n±; the corresponding operators a\\_ and a+ 
(al_ and a_ ) create and destroy one quantum with posi- 
tive (negative) circular polarization and one unit of posi- 
tive (negative) angular momentum. For the nonrotating 
system with Hamiltonian Hq, each quantum simply has 
an energy tuv±. 



The principal advantage of this particular circularly 
polarized basis appears for the rotating system, with 
Hamiltonian H' a = Hq — flL z , because Eqs. (|5.6[) and 
(|5.7p together imply 



H' Q = huj± + fku + a\\_a + + /kj_aLa_, 



(5.8) 



where uj± = lo±_ =p Q are the frequencies in Eq. (|5.3[) for 
the positive and negative classical orbits as seen in the 
rotating frame. Correspondingly, the energy eigenvalues 
are labeled with two non-negative integers n± 

e(n_|_, ri-) = n + h(ujj_ — fi) + n^h(ujj_ + f2), (5-9) 

where I omit the zero-point energy fko± for simplicity. 
In the limit of rapid rotations (f2 — > luj_), these eigen- 
values are essentially independent of n+, which implies 
a large degeneracy. The other integer n_ then becomes 
the Landau-lcvcl index, with different Landau levels sep- 
arated by an energy gap of ~ 2hu>±. 

For small angular velocity it is often convenient to 
introdu ce a different set of quantum numbers [pp. 727- 
741 of ICohen-Tannoudii et all (|l977t) ] n = n + + n_ for 
the energy and m — n + — n_ for the angular momentum, 
so that the energy eigenvalues have the equivalent form 



e = e rl 



h(n + to), \(n — to)] = nhujj_ — nihil. (5.10) 



Figure [TH] uses this new basis to illustrate the energy 
spectrum e nm for different values of Q. 

1. On the left side of the figure (fi = 0), the exci- 
tation energy is simply nhuj±, independent of to, 
forming an inverted pyramid of states: for each 
non-negative integer n, there are n + 1 degener- 
ate angular-momentum states ranging from — n to 
n in steps of 2. For small fi 1, the energy quan- 
tum number n determines the large splitting and 
the angular-momentum quantum number to lifts 
the remaining degeneracy (physically, this Coriolis- 
induced splitting is like the weak-held Zeeman split- 
ting of magnetic sublevels) . 

2. The central figure shows the situation for moder- 
ate fi, where the (now signficant) Coriolis effect 
eliminates the degeneracy of the nonrotating spec- 
trum. States with to = are unshifted, states with 
negative m shift up by |m|?lfi, whereas those with 
positive to shift down by the same amount. 

3. The right side shows the extreme situation for fi — > 
ui±, when the states again become nearly degener- 
ate, forming essentially horizontal rows (the Lan- 
dau levels). The lowest Landau level has n = to, 
which means n_ = and n + = to; the first excited 
Landau level has n — m — 1 (namely rt_ — 1 and 
n + = to + 1), etc. In this limit (fi < 1), the quan- 
tum numbers n+ and n_ of the rotating basis are 
most convenient. 
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FIG. 18 Excitation energy levels e nm = nhuij_ — mQ. from 
Eq. (|5.10[) for various angular velocity Q expressed in terms 
of quantum numbers n = n++n_ for energy and m = n+ — n_ 
for angular momentum. Left side is for £1 = 0, showing an 
inverted pyramid of states: for each non-negative integer n, 
there are n + 1 degenerate angular-momentum states ranging 
from — n to n in steps of 2. Center shows the situation for 
moderate Q. States with m = are unshifted, states with 
negative m shift up whereas those with positive m shift down. 
Right side shows the extreme situation for Q — > ui± , when the 
states form essentially horizontal rows (the Landau levels). 
The lowest Landau level has n — m, which means n_ =0, 
the first e xcited Landau le vel has n — m — 1 (namely n_ = 1), 
etc. From lDalibardl (|2007T) . 



Since the same LLL wave functions appear in the 
fractional quantum Hall effect for two-dimensional elec- 
trons in a strong magnetic field, this mean-field LLL 
regime is so metimes called the "m ean-field quantum 
Hall" regime ( Watanabe et all 12004 ). In fact, related 
many-body states from the fractional quantum Hall ef- 
fect play an important role for a rotating dilute Bose gas 
at still higher angular velocity, where the physics will 
turn out to be quite different (Sec. VI). Thus it seems 
preferable to reserve the name "quantum Hall regime" 
for this more unusual (nonsuperfluid) limit. 



B. Direct treatment of LLL states 

Before using the LLL wave function to construct a GP 
condensate wave function, I briefly discuss an alternative 
approach to the LLL states. It is natural to rewrite the 
single-particle Hamiltonian in the rotating frame Hq = 
p 2 / (2M) + |Ma; j_r 2 — £1 ■ r X p by completing the square 
to obtain two equivalent forms 



For rapid rotation (O < 1), it is appropriate to focus 
on states in the lowest Landau level (those with n_ = 0, 
denoted LLL), when the single quantum number n + = m 
suffices to specify the one-body states. The relevant en- 
ergies are the gap 2hoj± between adjacent Landau levels 
and the mean interaction energy gn(0) = fi, which is 
assumed small in the present limit. Thus the basic ex- 
pansion parameter of the theory is n/(2hu>±), although 
curr ent experiments have only achieved values as small as 
0.6 (ISchweikhard et all \2004) . The corresponding two- 
dimensional eigenfunctions i\) m (r) of both H' and L z 
have a very simple form in plane-polar coordinates [note 
that these are precisely the trial states from Eq. (|3.19|) 
used to study vortex arrays i n a w eakly interacting con- 
densate ( Butts and Rokhsarll 19991 )] 



ip m (r) oc r "V"^ e -<- 2 /(2di\ 



(5-11) 



where the normalization is unimportant for the present 
purpose. 

In particular, the ground state ipo is just a two- 
dimensional isotropic Gaussian. It represents the vacuum 
for both circularly polarized modes 



a± ip = 0. 



(5.12) 



The higher states ip m for positive m are proportional to 
(a^)' Tl ^o; they can be written in terms of £ = (x + 

iy)/d± = (r/d±)e i<)> as ip m cx ( m e^^ 24 ^, where the 
first factor is just a non-negative integral power of the 
dimensionlcss complex variable £. For a given m 3> 1, 
the density |i/'m('")| 2 is large only in a circular strip of 
radius ~ ^/rnd± and width ~ dj_/^/m. In addition, it 
is easy to see that the mean-square radius for any non- 
negative m is given by 



H 



# o = 



(p-Amxrf 1 , 



2M 



(p - Mwi x r ) 



2M 



+ -M(ujI - ft 2 ) r 2 (5.14) 



(o>x -«)•£, (5.15) 



where = luj_z. Either form is reminiscent of the non- 
relativistic Hamiltonian of a particle with charge q in a 
magnetic field B, where the momentum p and vector 
potential A appear in the gauge-invariant combination 
p - qA, with B = V X A. 

Landau solved the quantum problem of an electron in 
a uniform magnetic fiel d B [s ee, for example, Chap. XV 
of lLandau and Lifshitzl ( 19651 ). but the vector potential 
has a different (asy mmetric) form A = —By x]. I n con- 
trast, pp. 742-764 of lCohen-Tannoudii et all (|l977l ) solve 
the same problem in the symmetric gauge A — ^B X r. 
It is of course possible simply to quote those results to 
obtain the LLL wave functions ip m , but such an ap- 
proach requires a charge q, whose sign (positive or neg- 
ative) determines some crucial conventions in the choice 
of quantum-mechanical states. As written, Eq. (|5.L5|) 
suggests a unit positive charge with fictitious vector po- 
tential A = Mojj_ X r; the resulting uniform magnetic 
field is B = 2Mu)±. Unfortunately, the usual LLL pic- 
ture treats electrons with negative charge, and the cor- 
responding one-body states involve reinterpretations of 
some o f the quantum numbers , as di scussed in pp. 750- 
757 of ICohen-Tannoudii et al\ ( 1977| ). For this reason, 
I choose to work directly from the original Hamiltonian 



Hq = Hq — fl ■ L, because the role of : 



and 



is par- 



(r 2 



(m+l)d\. 



(5.13) 



ticularly transparent. 

Another interesting feature of the transformation to a 
quadratic form in p — qA is the possibility of topologi- 
cal gauge potentials. In the context of optical lattices, 
iJaksch and Zollerl (|2003h propose inducing various terms 
in the single-particle Hamiltonian that mimic the effect of 
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an external magnetic field. Specifically, when a two-level 
atom hops around a closed path in the lattice, external 
lasers can generate a net phase change that is essentially 
equivalent to the magnetic flux of an applied field [see , 
(iMuellerl . 12004 ISatiia et all . 120061 : ISorensen et oil . 120051 ) 
for alternative approaches, including the possibility of 
"non-Abelian" gauge potentials]. 

Simi lar ideas also apply t o a simple trapped conden- 
sate. lJuzeliunas et all (2005) propose using a three-level 
system in the "A" configuration [two nearly degenerate 
lower levels, each coupled to a high er third state by exter- 
nal laser beams ( Arimondol . Il996[ ); for a "dressed-atom 
approach" to t hese three-level systems , see, fo r example, 
pp. 451-454 of [Cohen-Tannoudii et all (|l998t) ]. Such an 
interacting three-level system has a "dark" state that ef- 
fectively uncouples from the remaining states. If one of 
the external lasers has orbital angular mom entum [see, 
for example , Sec. 52 of iGottfriedl ( 19661 ) and I Allen et all 
(|2003Lll999h ]. ;hen the dark state experiences a nontriv- 
ial topological gauge potential that in principle can act 
like an essentially arbitrary magnetic field (and hence an 
arbitrary applied vorticity for the trapped condensate). 
Whether such methods can indeed generate angular mo- 
mentum in the condensate and the associated vortices 
remains uncertain. More work (both theory and experi- 
ment) is clearly desirable. 



C. LLL condensate wave functions 



IHqI 1)20011 ) proposes a linear combination of the LLL 
single-particle states in Eq. (|5.11[) as the GP condensate 
wave function for the rapidly rotating two-dimensional 
BEC 

4>LLL{r) = c m ^m(r) = f(C)e- r2/{2dl \ (5.16) 

m>0 

where /(£) = J2 m >o c mC m 1S an analytic function of the 
compl ex variab l e C~ [c ompare Sec. III.C; see also pp. 587- 
592 o f lShankari(ll994h ]. Specifically, for a truncated basis 
set, the analytic function /(£) is a complex polynomial 
and thus has a factorizcd form 



/(Oocncc-o), 



(5.17) 



apart from an overall constant factor. Evidently, /(£) 
vanishes at each of the points Q , which are the positions 
of the nodes of the condensate wave function ipLLL(r). 
In addition, the phase of this wave function increases by 
2ir whenever £ moves in the positive sense around any of 
these zeros. Thus the points Cj are precisely the positions 
of the vortices in the trial state ipLLL, and minimization 
with respect to the constants c m is effectively the same 
as mi nimization with respect to the position of the vor- 



tices (Butts and Rokhsa: 



spect to 
r, 1999). 



1. Unrestricted minimization 

This trial state has some very interesting and un- 
usual properties. Apart from the overall Gaussian en- 
velope, the spacing of the vortices completely determines 
the spatial variation of the number density nhLhir) = 
\ipLLL(r)\ 2 — |/(C)| 2e ~ r ' l d '^. The number density varies 
smoothly between adjacent vortices because |/(C)| 2 con- 
sists of harmonic functions. Hence the vortex core size is 
comparable with I = \J h/(MQ), which is essentially d± 
in the rapidly rotating limit. In contrast to the Thomas- 
Fermi limit at lower angular velocity, the present trial 
function ipLLL automatically incorporates all the kinetic 
energy. 

It is important to emphasize that the one-body state 
tJjlll arises from the GP energy functional and that the 
corresponding many-body state 



N 



(5.18) 



is just a Hartree product of these LLL states. For any 
such ground state, the system exhibits "off-di a gonal" 
long- range order (jPenrose and Onsagen . Il956l lYangl . 
Il962l ). and i^lll is the corresponding macroscopically 
occupied condensate wave function. Thus the mean-field 
LLL regime still has a BEC and, in general, exhibits su- 
perfluidity. 

In the limit of rapid rotation, this LLL wave func- 
tion can serve as a possible variational trial solution 
for the GP energy functional. Let denote 
an expectation value evaluated with tpLLL- Equations 
flS/F} and (I5.12|) show that (L z ) LLL /h = (a ] + a+) LLL be- 
cause (a_a,-)LLL vanishes for any LLL state. Further- 
more, the quantity r 2 = x 2 + y 2 is easily expressed 
as (x + iy)(x — iy) = d\{a- + a^)(a+ + al_), whose 

expectation value gives (r 2 ) lll/ d 2 ^ = (a+a+)ii,i + 1. 
Comparison yields the remarkably simple relation be- 
tween the_aiigulBi^mon^ mean-square ra- 
dius (|Aftalion et all , 120051 IHoI . l200ll IWatanabe et all . 
2004) ' 



(L Z ) LLL = Muj ± (r 2 ) LLL - h 



(5.19) 



for any linear combination of single-particle states in the 
lowest Landau level [as discussed in Sec. III.C, a slightly 
di fferent form of this re l ation is important in the analysis 
of iButts and Rokhsarl ( 1999) )]. Note that this relation 
also follows directly from Eq. (|5.13[) . 

A combination of Eqs. (|5.2|) and (|5.19p expresses the 
full GP energy functional (|5.ip in the very intuitive and 
familiar form 



e'[^ lll } = nn + jd 2 r [Muj\ (i - n) r 2 \i; LLL \ 



\92D\lp. 



LLL 



(5.20) 



where the integrand again involves only r 2 | | 2 and 
IV'lllI 4 [compare Eq. (|2.24[) ]. Unrestricted minimization 
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with respect to \iPlll\ 2 at fixed normalization readily 
yields the number density 



|V>min(r)| 2 = n min (r) = n min (0) 1 



7?2 ; ' 



(5.21) 



where the (two-dimensional) central density is n m i n (0) = 
Mmin/ff2D, and the condensate radius is given by 



R 2 = 



/Anil 



(5.22) 



Remarkably, this expression has the same parabolic form 
as the usual TF result in Eq. (|2.27p . even though the 
present functional explicitly includes all the kinetic en- 
ergy. It does, however, rely on the restriction to the 
lowest Landau level. 

The normalization condition J d 2 r |'0min| 2 = N yields 
N = 7j7ri?Qn m j n (0), which is equivalent to 



R 



2Ng 2D 

TT/^min 



(5.23) 



The product of Eqs. (|5.22|) and (|5.23|) gives an explicit 
expression for Rq 



Rq 



2Ng 2D 



&Nad 4 , 



(5.24) 



that explicitly shows the radial expansion for Q — > 1. 
Similarly, the ratio of the same two equations yields the 
chemical potential 



2Xu-> D M^ 2 (1 - <>) _ ^ ;8Na(l - V.) 



hw i 



(5.25) 

The first form is similar to Eq. (|4.9p in the mean-field 
TF regime (the different fractional powers arise from 
the different dimensions for the normalization integrals), 
whereas the second form gives the condition /i m i n < 2Scjj_ 
for the validity of the mean-field LLL regime 



1 - < 



2Na 



(5.26) 



If Z/2a ~ 100 and N ~ 10 4 , then a dimensionless rota- 
tion speed f2 ~ 0.99 just barely approaches the mean-field 
LLL regime. 

The energy of this approximate minimizing mean-field 
LLL density |?/> m in| 2 follows immediately by evaluating 
Eq. IfOO)) 



E' min - nn 2^ 



Nhiu 



= 4 j2JVo(l-n) 
3Huj± 3 



(5.27) 



As emphasized bv lAftalion et all ( 20051 ). the approximate 
density in Eq. (|5.2ip arises from an unrestricted varia- 
tion and is not within the LLL. Thus the actual energy 
will be higher, because the minimization then includes 
the vortices through the zeros of the analytic function in 



Eq. (|5.17p . As seen below, the main effect of the vor- 
tices is to increase the numerical coefficient in Eq. (|5.27p , 
keeping the same dependence on the various parameters. 

The derivative of Eq. (|5.27[) with respect to £1 yields 
the angular momentum of the minimizing mean-field LLL 
state 



(£*)min = lMNRlu ± -h, 



(5.28) 



which also follows from Eq. (|5.19j) . Here, the first term 
is just the solid-body value (since direct integration with 
the parabolic IV'minl 2 yields (r 2 ) m i n = jNRq). Unlike the 
mean-field TF regime discussed in Sec. IV. A, the LLL 
reduction from the classical value is only of order N^ 1 
and thus negligible in the thermodynamic limit. 



2. Inclusion of the vortices 



The minimizing parabolic density profile in Eq. (|5.21|) 
cannot provide a full description of the rotating LLL con- 
densate, for it omits the fine-grained structure associated 
with the vortices. 



a. renormalization of g2D In the present mean-field LLL 
regime, the intervortex separation ~ d± is small com- 
pared to the size of the condensate. Thus it is 
possible to write the number density obtained from 
Eq. (|5.16p in the form \ip L LL{r)\ 2 « \4>LLL{r)\ 2 \m(r)\ 2 
slowly varying envelope func- 
a rapidly varying modulation 
vanishes at the center o f each 
2005; Bavm and PethickL [2004 



as the product of a 
tion \4iLLh{r)\ 2 and 
functio n \m(r)\ 2 that 
vortex (lAftalion et all. 

i Fetterl . l200lt iFischer and Bavml . 120031 : IWatanabe et all , 
2004). For a large vortex lattice, this modulating func- 
tion is essentially periodic, and it is convenient to nor- 
malize it so that J ceU d 2 r \m(r)\ 2 = 1 over each unit cell. 
Substitution of this approximate density into the LLL 
energy functional (|5.20p yields products of slowly vary- 
ing functions involving powers of \^lll{t)\ 2 and rapidly 
varying periodic modulating functions involving powers 
of |m(r)| 2 . 

The integral over the condensate separates into a sum 
of integrals over each unit cell, with the slowly varying 
quantities acting as locally constant factors. Thus the 
only effect of the vortex lattice is that the interaction 
term acquires a numerical factor (3 = f n d 2 r |m(r)| 4 , 
which effectiv ely renormalizes the interaction constant 



92D — * flg2D Aftalion et all (|2005h discuss these steps 



carefully]. For an unbounded triangular vortex lattice, 
the resulting (3 w 1.1596 is the numer ical value for the 
corresponding Abrik osov vortex lattice ( Abrikoso"vlll957t 
iKleiner et al. , Il964h . Consequently, the radius Rq in 
Eq. ((OH) expands by the factor (3 1 / 4 w 1.0377. 



b. nonuniformity of the vortex a rray To s tudy the effect 
of the vortices in more detail, iHol (|200lh considers the 



3G 



logarithm of the particle density for a LLL condensate 
wave function of the form in Eq. (|5.16p , obtaining 



ln|-0Lii(r)| 2 = \nn LLL (r) = --^ + 2 ^ In \r - rj\, 

j 

(5.29) 

where rj is the position of the jth vortex. The two- 
dimensional Laplacian of this equation yields 

V 2 lnn LLL (r) = -4- + 4tt^ 6™ (r - r 3 ) (5.30) 



because V 2 In \r — Tj\ — 2tt8^ (r — rj). The sum over 
the delta functions is precisely the two-dimensional vor- 
tex density n v (r), which gives a striking relation be- 
tween the particle density fiLLL if) in the mean- field LLL 
regime and the vortex densit y ( Aftalion et ali . 120051 : \~ 
120011: IWatanabe et al. 1. 12004) 



n v (r) 



4tt 



V 2 h\n LLL (r). 



(5.31) 



The first term is an effective vortex density n c ff = 
Muj±/ (nh) for £1 = and the second is a (small) cor- 
rection. 

To appreciate the implications of this elegant equa- 
tion, first assume that the vortex density is uniform 
with the actual Feynman value im v — mQ/h. Equa- 
tion (|5.3ip then has the Gaussian solution riLLLif) = 
it-lllIO) exp(— r 2 /<7 2 ), where o~ 2 = 7r(n ff — n v ); equiva- 
lently, the effective squared condensate radius is 



rf 2 



m (oj± - n) 1 - o ' 



(5.32) 



As expected from the radial expansion, the length a di- 
verges as fl — > u)±. 

In fact. lCooper et al. (2004) use numerical methods to 
show that a small distortion of the vortex lattice can both 
lower the energy and significantly change the mean den- 
sity profile from Gaussian to one closely resembling the 
parabolic shape in Eq. (|5.2ip . Thus it is better to substi- 
tute this profile into Eq. (|5.31[) to find t he approximate 
and s lightly nonuniform vor tex density ( Aftalion et ali , 
2005; Watanabe et al. I. I2004T ) 



n v (r) 



1 



% {l-r 2 /RlY 



(5.33) 



This expression is similar to the result in the mean-field 
Thomas- Fermi limit in Eq. (|4.24|) . and the correction is 
again small, of order 1/N V . 



3. Rapidly rotating anisotropic traps 

In practice, most traps are not exactly axisymmetric 
(namely, lj x 7^ uj y , and I take u> x < uj y for definite- 
ness). The rotational properties of such realistic traps 



differ in certain significant ways from those of symmetric 
traps. Similar probl ems arise in the study of deformed 
nuclei ( Valatinl . ll956f ). If the anistropic trap is stationary, 



the noninteracting single-particle ground state is simply 
a product of two one-dimensional Gaussians in x and y 
with oscillator lengths d x and d y . With increasing rota- 
tion speed f2, however, the coupling term —flL z mixes 
the x and y components. As a result, the ground-state 
density elongates along the direction of w eak confine- 
ment (here along x) and diverges as fl — > ui x ( Linn et ~al\ . 
I2001D . 

In complete analogy with the operators a± and aj_ for 
the symmetric trap [see Eq. f|5 . 5[) ] , a rotating anisotropic 
trap has a similar set of creation and destruct ion op- 
erators, but their fo rm is rather more intricate (|Fetterl . 
120071 : IValatml I 1956). Nevertheless, in the limit Q — » ui x , 
the lowest single-particle states constitute a nearly de- 
generate set that again acts like the lowest Landau level. 
The corresponding LLL condensate wave function again 
involves a polynomial in a scaled complex variable (a lin- 
ear combination of x and iy). Consequently, many of 
the results derived for the mean-field LLL regime in an 
isotropic trap continue to apply, at least qualitatively, for 
the rapidly rotating anisotropic trap. 

One important difference is the significan t anisotropy 
of th e condens ate as Q — > u> x (Fetter . 2007 : iLinn et ali 
2001). Indeed, Sanchez-Lotero and Palaciosl ( 20051 ) and 



Sinha and Shlvapnikovl (|2005f ) both model this situation 



with a one-dimensional effectively infinite strip of finite 
width, containing one or more rows of vortices. As dis- 
cussed in Sec. VI below, however, this behavior may well 
be unobservable. Before reaching such a one-dimensional 
limit, there may well be a quantum phase transition to 
a correlated many-body state (sometimes called a "vor- 
tex liquid") that is not superfluid and has no BEC. Such 
questions currently remain unanswered and merit more 
study. 



VI. QUANTUM PHASE TRANSITION TO HIGHLY 
CORRELATED STATES 

The mean-field LLL regime relies on the assumption 
that the interaction energy per particle ~ /i is small com- 
pared to the gap 2hu>±, leading to the lower bound for Ct 
[compare Eq. (|5.26|) ] 



1 - 



Z 
2N~d 



< n 



(6.1) 



where a = /3a is the renormalized scattering length that 
includes the vortex-induced density variation (note that 
(3 pa 1.1596 is simply a numerical factor, see Sec. V.C.2 
for more details). In practice, typical experiments require 
f2 > 0.99 to reach this regime. For the present purpose, 
it is important to emphasize that the mean-field LLL 
regime has a BEC with macroscopic occupation and is 
indeed a superfluid. Specifically, the associated many- 
body GP ground state is simply the Hartree product in 
Eq. HETgp . 
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For still larger values of Cl, however, the mean-field 
LLL regime should eventually disappear through a quan- 
tum phase transition, leading to a wholly different highly 
correlated many-body ground state ^corr- Unlike the 
mean-field LLL ground state, the new ground state ^/corr 
does not have a BEC and is not a superfluid. As ev- 
idence for such a transition, the exact two-dimensional 
ground state for bosons with small N and relatively 
large fixed angular momentum L 7 indeed has a corre 



lated form (|Wilkin and Gunnl . 120001 : IWilkin et all 1 1998ft . 
Furth ermore, detailed numerical studies ( Cooper et all 
l200ll ) explicitly exhibit such a phase transition. 

To understand the physics of this quantum phase tran- 
sition, it is helpful to focus on the number of vortices Af v . 
In the mean-field LLL regime, the two-dimensional con- 
densate has a radius Rq, and Eq. (|4.6p yields the number 
of vortices Af v = R 2 Q /l 2 w R%/d 2 ± . Use of Eq. (I5T24)) and 
the discussion in Sec. V.C.2 give 



R o 

dl 



8Na 



Z(l - ft)' 



(6.2) 



As a n alternative derivation of this result (jBloch et all 
2007), note that the mth LLL state Vm has a mean- 
square radius ~ md\ for m ^> 1 [Eq. (|5. 13[) ] . If the 
sum in Eq. (|5 . 16[) extends to m max , then the conden- 
sate has an effective squared radius ("r 2 )m max ~ TOmaxdj_- 
The number of vortices (the degree of the polynomial) is 
Af v = TOmax, which is the same as Eq. (|6.2[) found directly 
from the Feynman relation for the vortex density. The 
parameters used below Eq. ()5.26|) yield Af v rs 200. 

It is valuable t o con s ider the ratio v = 

N/Af v (ICooper et all l200lh IWilkin and Gunnl . l2000ft 
known as the "filling factor" or "filling fraction" because 
of a similar quantity in the fracti onal quantum Hall 
effect [see, for example, Chap. 1 of iPrange and Girvml 
( 1987ft ]. In the mean-field LLL regime, it has the simple 
form 



N 



(6.3) 



but it also remains well defined for general many- 
body ground states of rotating bosons [and for charged 
fermio ns in a magnetic fiel d, which is its original 
role ( Prange and Girvml Il987ft ]. 

In the mean-field LLL approach for weakly interact- 
ing bosons in a ro t ating trap, typified by the work of 
iButts and Rokhsarl ( 1999ft and discussed in Sees. III.C 
and V.C, the equilibrium state is a vortex array 
that breaks rotational symmetry and is not an eigen - 
state of L z . As an alternative, ICooper et all (2001) 
use exact diagonalization to study the ground states 
for Af v = 4, 6, and 8 and for many values of N 
[see also (ICooper and Wilkinl . Il999l : IViefers et al] , l2000t 
IWilkin and Gunnl . l2000ft for less extensive studies]. This 
numerical analysis uses a toroidal geometry with periodic 
boundary conditions in both directions, allowing a study 
of bosonic systems with relatively large integral values of 



both N and Af v . For a particular aspect ratio of the torus 
that can accommodate a triangular array with Af v = 8, 
they investigate filling fractions v ranging through sev- 
eral rational values from ^ up to 9. They also study the 
mean-field GP ground state containing a vortex lattice 
with Af v = 8 for the same values of v. Comparison of 
these two classes of states shows that the GP vortex lat- 
tice is the ground state for v c > 6, where the specific 
numerical value depends on the aspect ratio of the torus 
(and hence the number of vortices that can fit into the 
periodic structure). Other choices of this aspect ratio 
yield critical values in the range v c ~ 6-10. 

A combination of Eq. (|6.ip for the validity of the 
mean-field LLL regime and the lower bound v > 1 
for the mean-field GP ground state yields the inequal- 
ities t hat define th e allow ed range of the mean-field LLL 
state (jBloch et all l2007ft 



2Na ~ ~ ZN 



^7 (6-4) 



The right-hand inequality for ft illustrates the experi- 
mental difficulty of studying the quantum phase transi- 
tion. Since Z/a typically exceeds 100, Eq. (|6.3|) makes 
it clear that reaching even v ~ 10 will require signifi- 
cant reduction of N (note that the high rotation speed 
ft ~ 0.999 presents a daunting experimental challenge). 

In contrast to the vortex lattices for v > v Cl the ground 
states for smaller v < v c are rotationally symmetric in- 
compressible vortex liquids that are eigenstates of L z . 
They have close similarities to the bosonic analogs of the 
I J ami (|l989ft sequence of fractional quantum Hall states 



that describe a two-dimensional electron gas in a strong 
magnetic field. The simplest of these many-body gro und 
states is the bosonic Laughlin state (jLaughlinl . Il983ft 

*Laughlin(T , l,r 2 , ■•• , V N ) 



N 



II( 

n<n' 



where z n = (x n + iy n )/d± refers to the nth particle. Here 
the power 2 ensures that this state is symmetric under 
the interchange of any two particles (in contrast to th e 
fermionic case wit h powe r 31. ICooper and Wilkinl (Il999ft: 
IWilkin and Gunnl (|2000ft ; IWilkin et all (|l998ft note that 
this state occurs for v = i. 

It is worth emphasizing the difference between this 
Laughlin ground state and the GP Hartree ground state 
in Eq. ([5TT8]) . For the Laughlin state in Eq. (|6T5]) . 
the first (product) factor means that there is no off- 
diagonal lo ng-range order ( Penrose and Onsagerl . Il956t 
lYand. 11962ft and hence no BEC. In addition, the Laugh- 
lin state vanishes whenever two particles come to- 
gether, enforcing the many-body correlations. Indeed, 
the usual short-range two-body potential V(r — r') = 
g5(r — r') has zero expectation va lue in this correlated 
state ( Trugman and Kivelsonl . [i985ft . In contrast, the GP 
ground state has every particle in the same single-particle 
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Hartree state ipLLL j thus particles in the GP condensate 
te nd to overlap where ver \ipLLL( r )\ is large. 

ISinova et all ( 20021 ) provide an appealing physical pic- 
ture of the onset of this quantum phase transition. Us- 
ing various approximate descriptions, they find that the 
long-wavelength spectrum of collective modes becomes 
quadratic in the limit of rapid rotations, in contrast to the 
familiar linear Bogoliubov spectrum in Eq. (|2.44|1 . This 
softened spectrum implies that the vortex lattice melts 
for a critical filling fraction v r 8, in re a sonab le agree- 
ment with the estimates of ICooper et all (|200lh . An al- 
terna tive description of the vortex-lattice melting ( Bayml . 
I2005D relies on the softening of the Tkachenko vortex- 
lattice modes, as seen in Eq. (|4.28|) . The connection 
between these apparently distinct pictures remains un- 
certain. 

The periodic toroidal geometry of lCooper et all ( 200lh 
has the strong advantage of avoidin g boundary effects; 
iRegnault and Jolicoueurl (|2003l[2003 ) use a spherical ge- 
ometry to perform a similar numerical analysis. Both 
groups propose specific sequences of known quantum Hall 
states for the rotating dilute Bose gas when i < v < v c , 
although there is currently no experimental e viden ce for 
such sequences. Nevertheless, ICooper et al. I (|200ll ) find 
strong overlap between their results from exact diagonal- 
ization and the Laughlin state {y = 4) and other similar 
but more complicated quantum Hall states for the partic- 
ular filling fractions v = 1, |, 2, |, 6. For a more detailed 
discussion of these various quantum Hall states, see, for 
exam ple. (ICooper et all , l200ll; IRegnault and Jolicoueurl 
2003, 200f and Sec. VII. C of lBloch et all (1200 7l) . 



What is the physics of this ground-state transition 
from a vortex lattice to a vortex liquid? Why is the crit- 
ical value of the filling fraction v c ~ 5 — 10, as opposed 
to v c ~ 1 or v c ~ 0.1? To think about this question, 
it is helpful to consider N bosonic particles in a plane, 
with 2N degrees of freedom. Vortices appear as the sys- 
tem rotates, and the corresponding vortex coordinates 
provide Af v collective degrees of freedom. 12 This relation 
between particle and collective degrees of freedom is fa- 
miliar in many areas of many-body theory: For example, 
the Debye theory of specific heat uses the phonon normal 
modes instead of the particle positions in the crystal, and 
the original number of particle degrees of freedom (3./V) 
also determines the cutoff frequenc y for the phonon spec- 
trum [see, for example, Chap. 23 of lAshcroft and Merminl 
( 19761 )]. Similarly, the theory of a charged electron gas 
often relies on collective density modes ( plasma oscil la- 
tions) instead of the particle coordinates (|Pinesl [l96lh . 

For a slowly rotating Bose-Einstein gas in a plane, the 
2N particle coordinates provide a convenient description. 
In principle, the Af v collective vortex degrees of freedom 



12 Note that vortices obey first-order dynamical equations, so that 
the associated x and y coordinates are not independent; in par- 
ticula r, they can serve as a pair of Hamilt onian canonical vari- 
ables jFetten , ll967tlHaldane and Wul . ll985t) . 



should reduce the original total 2N degrees of freedom 
to 2N — AT V , but the effect is unimportant as long as 
N v <C N . When M v becomes comparable with N, how- 
ever, the depletion of the particle degrees of freedom be- 
comes crucial, ultimately producing a phase transition 
to a wholly different ground state. With this picture, 
the critical value v c ~ 5 — 10 is not surprising, since 
N v /N = 1/v is still relatively small. 

If experiments can indeed create these correlated 
states, how might they be detected? Section VII. C of 
iBloch etall ( 2007h proposes several techniques, includ- 
ing (1) a reduced three-body recombination rate [be- 
cause of the anti-correlation between the bosons inherent 
in the first factor in Eq. (|6.5p ] and (2) the characteris- 
tic density profile of these various quantum Hall states. 
Another possibility would be to study th e ensemble av- 
eraged density-density noise co rrelations ( Altman et all , 
12004 iBrown and Twissl . Il956f ). which should exhibit a 
different structure in the correlated state from that in 
the mean-field LLL regime. 



VII. OUTLOOK 

As this review demonstrates in detail, experimental 
and theoretical studies of rotating Bose-Einstein conden- 
sates generally agree well for scaled angular velocities 
= fl/uj± up to - 0.995. 



A. Successes 

These accomplishments include several different 
regimes. 

1. For slow rotation f2 = Cl/ui± -c 1 (Sec. Ill), the 
condensate includes only a few vortices, and cen- 
trifugal effects are negligible. Thus the overall den- 
sity profile is essentially that of the nonrotating 
condensate, apart from the small holes associated 
with the vortex cores. Experimental observations 
of vortex precession in one- and two-component su- 
perfluids (see Sec. III.D.l) provide the only direct 
measurements of vortex dynamics in dilute BECs. 

2. As the rotation increases, the centrifugal effects be- 
come important, leading to significant changes in 
the condensate's aspect ratio (Sec. IV). For not too 
large f2, however, the intervortex spacing ~ / — 

1/2 

[H/(Mu)] remains large compared to the vor- 
tex core size £. In this mean-field Thomas-Fermi 
regime, the gradient energy associated with den- 
sity variation remains negligible, and the hydrody- 
namic flow velocity v = KVS/M predominates in 
the energy balance. Typically this case holds for 
0-75 < Q < 0.99. 

3. Finally, for very large rotation speed 0.99 < fi < 
0.999, the vortex cores expand, filling much of the 
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available space. The energy associated with spa- 
tial variation of the density now becomes compa- 
rable to the hydrodynamic flow energy. In this 
mean- field lowest-Landau-level regime (Sec. V), the 
single-particle part of the energy in the rotating 
frame is exactly soluble, leading to the rotational 
analog of the Landau levels for an electron in a 
uniform magnetic field. In particular, there is an 
energy gap ~ 2hu>± between the nearly degenerate 
lowest Landau level and the first excited Landau 
level. In addition, the condensate expands radi- 
ally, decreasing both the central particle density 
and the mean interaction energy [i, which thus be- 
comes smaller than the energy gap 2ftwj_. The set 
of lowest-Landau-level states then provides a con- 
venient and flexible description of the equilibrium 
of the rapidly rotating condensate. Although cur- 
rent experiments with fi / (2huj j_) « 0.6 only just 
reach this mean- field LLL regime (sec Sec. V.A), it 
is reasonably well understood. 



B. Challenges 

Despite the many remarkable achievements of the past 
decade, several theoretical predictions of considerable in- 
terest remain unverified experimentally. 



1. Rapid rotation of both Bose and Fermi gases in the 
quantum Hall regime at small filling fraction 

For dilute ultracold trapped single-component Bose 
gases of the sort considered in this review, one obvious 
challenge for future experiments is the quantum phase 
transition from a superfluid with a Bose-Einstein con- 
densate to a nonsuperfluid correlated state. As discussed 
in Sec. VI, formation of this new state requires large an- 
gular momentum and smaller particle number N than 
is typical in most current experiments. Whether such 
states can indeed be realized is uncertain. In principle, 
one might start with an anharmonic trap of the sort con- 
sidered in Sec. IV. C. Once the condensate rotates with 
fl > 1, it might be feasible adiabatically to turn off the 
external laser beam (and hence the anharmonic compo- 
nent of the trap potential), raising the effective oj± and 
leaving the system with high angular velocity fi < 1. 

Another possibility is be to rely on topological gauge 
potentials (Sec. V.B), in which applied laser fields with 
orbital angular momentum can produce effective nonuni- 
form topological "magnetic" fields that are more general 
than the usual constant B. To my knowledge, all studies 
of rapidly rotating trapped bosons rely on the detailed 
form of the operator — ft ■ L = —A p, with A = fl X r as 
the equivalent vector potential. The existence of a non- 
superfluid correlated many-body state reflects the high 
degeneracy of the lowest Landau level, which in turn de- 
pends on the specific form of the effective vector poten- 
tial (A = fi X r, apart from a possible gauge transfor- 



mation). It is not obvious that such nonsuperfluid cor- 
related many-body states occur for more general topo- 
logical gauge potentials. This question certainly merits 
additional detailed study. 

A related area for future activity is the study of 
rapidly rotating Fermi gases. At present, this the- 
oretical topic is only partia l ly explored [see, for ex- 
ample. dBhongale et all l2007t iMoller and Cooped . l2007t 
IZhai and Ho . 2006 )]. and much work remains. The 
only r elevant experiments are those of IZwierlein et all 
(120051) on the vortex lattice in rotating fermionic 6 Li 
(see Sec. IV. B. 6). The experiment starts with a BEC of 
tightly bound bosonic molecules and then sweeps across 
the Feshbach resonance toward the more weakly bound 
fermionic regime. Unfortunately, any move far to the 
BCS side of the transition at fixed low temperature nec- 
essarily encounters the severely reduced critical temper- 
ature of the weak-coupling BCS limit. 

As a pure theoretical argument, let us consider the 
weak-coupling fermionic ground state as a function of 
f2 = n/u±. For in t erme diate f2, a vortex lattice forms, as 
in lZwierlein et al. (2005). What happens for larger f2 ? Is 
there a transition to a normal state similar to that at H C 2 
in a type-II sup erconductor [see, for example, Sees. 4.7, 
4.8, and 5.3.3 of lTinkhaml (|l996l )]? Or does a transition 
to a correlated quantum Hall state intervene? If so, for 
reasonable experimental parameters, how does the fill- 
ing fraction v relate to Q ? These intriguing questions 
deserve careful study. 



2. Rotation of ultracold dilute gases with long-range 
dipole-dipole interactions 

Dilute quantum gases with long-range anisotropic 
dipole-dipole interactions constitute a fascinating gen- 
eralization of the usual case with short-range isotropic 
interactions characterized by a scattering length a. Al- 
though the general dipole-dipole potential is rather 
complicated, experiments usually use a strong uniform 
field to orient the dipoles (either magnetic or electric) 
along the z direction. Such fields lead to the sim- 
ple r interaction potential [s e e, for example, the review 
by iMenotti and Lewensteml (12007) and the experimen- 



tal studies of 52 Cr 
120071 : lLahave et all 

V dd (r - r') = 



Griesmaier et all 120051 iKoch et all 
2007h ] 



r — r 



r\3 



1 - 3(£ ■ n)' 



(7.1) 



where d is the relevant dipole moment and n is a unit 
vector along the direction r — r' joining the two atoms. 

The dipole-dipole potential in Eq. (|7.ip is repulsive if 
z ■ n vanishes (namely if the oriented dipoles lie in the xy 
plane), whereas it is attractive if (i • n) 2 > | (in partic- 
ular, for oriented dipoles located along the z axis). Since 
a BEC with attractive interactions tends to collapse, a 
long cigar-shaped condensate with dipole-dipole interac- 
tions should be unstable, but a flat disk-shaped conden- 
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sate should be stable. iKoch et al. I (|2007l) use a Feshbach 
resonance to decrease the s-wave scattering length, thus 
enhancing the relative role of Vdd- These experiments 
confirm the qualitative expectations in great detail. 

Recent theoretical studies have predicted many prop- 
erties of the ground states of rotating BEC s with dipolar 
intera ctions. In t he weak-interaction limi t, ICooper et al\ 
(|2005h [see also (|Zhang and Zhail . I2005T) ] propose that 
an increased dipolar coupling strength induces a series 
of transitions in the mean-field (large filling fraction v) 
condensate involving vortex lattices of varying symmetry 
(triangular, square, stripe, and "bubble" phases). They 
also study the behavior for small v . when the system is no 
longer superfluid. Subsequently, iKomineas and Cooperl 
(|2007h analyze the more general situation as a function 
of both the dipolar interaction and the chemical potential 
relative to the s-wave contact interaction. 

Another system of great interest is the rapidly 
rotati ng dipolar Fermi ga s, whose predicted ground 
state ( Baranov et all [20051 ) has a simple incompressible 
form for filling fraction v = |. For smaller u, however, 
theory predicts that the quantum Hall states ultimately 
become unstable to the form ation of crystalline Wigner 
states (|Osterloh et all l2007t h Experimental studies of 
such nonsuperfluid fermionic states involve many techni- 
cal issues that remain for future investigation. 



3. Rotating spinor condensates 

The development of o ptical laser traps that confin e all 
the magnetic substates ( Stamper-Kurn et all . Il998l ) al- 
lows a direct study of spin-one BECs (as well as higher- 
spin cases). These experiments use a focused (red- 
detuned) infrared laser that attracts the BEC to the 
waist of the beam. For a general theo retical introduction, 
see Se c. III.E, alo ng with Sec. 1 2.2 of iPethick and Smithl 
(l2002l). Se c. 3.4 oflFetterl (l2002h. ari d the o riginal papers 
iHol (ll998| ) and lOhmi and Machidal (jl998h . In contrast 
to the two-component mixtures of the different hyperfine 
states F = 1 and F = 2 (see Sees. III.D.l and IV.B.5), 
the rotational invariance of the interaction between two 
atoms with the same F = 1 now imposes special restric- 
tions on the allowed interaction potential. The macro- 
scopic order parameter becomes a spin-one object with 
three components 




(7.2) 



labeled by the mp value. In the low-energy limit where 
only s-wave scattering is relevant, the effective short- 
range interaction has the form 

Mnt(n - r 2 ) = (.g + g 2 F 1 ■ F 2 ) <f< 3 >(n - r 2 ), (7.3) 

where Fj is the hyperfine spin of the jth atom. The cou- 
pling constant go is generally positive, but g 2 is propor- 
tional to a,2 — ao (the difference of the scattering lengths 



in the two symmetric channels F = F\ + F 2 — or 2). 
Thus g 2 can be either positive or negative, depending on 
the details of the two-body scattering. 

A mean-field description (|Hol . Il998l ) uses an effective 
energy functional and writes the spinor order parame- 
ter in a factored form *f? a {r) = y/n(r) ( a (r). Here, a 
runs over the three allowed values a = (1, 0, —1), n(r) is 
the assumed common density for all three spin compo- 
nents, and Ca( r ) is a normalized spinor with C f • C = 1- 
Apart from the gradient term in the effective energy, 
the ground-state spinor follows by minimizing the spin- 
dependent part of the energy density ^n 2 g 2 (F} 2 , where 

If g 2 is positive, then the minimum energy occurs for 
| = 0. Such states are known as "polar" or "anti- 
ferromagnetic," and they follow by spatial rotations of 
the h yperfine state | mp = 0) [this situation applies for 
23 Na (|Stenger et all Il998| )]. If, however, g 2 is negative, 
then the minimum energy is for a "ferromagnetic" state 
with | (F) |=1, obtained by spatial rotatio ns of the state 
\mp = 1) [this situation applie s for 87 Rb ( Sadler et all 
l2006t IVengalattore et ^H . l2007h ]. 

The multicomponent structure of the spinor order pa- 
rameter in Eq. (|7.2p is reminiscent of the order param- 
eter for superfluid 3 He, where each atom has a nuclear 
spin 1/2 and zero net electronic spin. Low-temperature 
experiments show that 3 He forms a p-wave paired super- 
fluid, with the two atoms in a triplet nuclear spin state 
and un it relative orbital angu l ar mo mentum [see, for ex- 
ample, ^nhardLandW61flg (199(3)]. In particular, the 
3 He order parameter has several components, analogous 
to those for a spinor BEC. Many groups have studied 
the structure of vortices in rotating 3 He, with various 
proposed candidates, both for a single vortex a nd fo r 
vortex lattices [for d e tailed reviews, see iFetterl (|l98fih : 
Salo maa and Volovikl (| 19871 )] . 

Unfortunately, the semi-phenomenological Ginzburg- 
Landau theory for superfluid 3 He is often inadequate to 
predict the detailed experimental phase diagram for var- 
ious vortices as a function of pressure and temperature. 
Specifically, the energy difference between many of the 
vortex states is small, and the relevant phenomenologi- 
cal Ginzburg-Landau parameters have considerable un- 
certainties their ex perimental values [see , for e xample, 
Sees. 7.5 and 7.6 of IVollhardt and Wolfld (|l990l) ]. Con- 
sequently, experimental studie s provide ess e ntial guid- 
ance for theoretic al analyses. lEltsov et al .] (|2005h and 
IFinne et all (2006) review recent experimental work on 
vortices and other related structures in superfuid 3 He. 

In an unbounded dilute trapped BEC with a single- 
component order parameter, the vortex core size and 
structure arise from the balance of the kinetic energy 
and the interaction energy (see Sec. III. A). For a trapped 
BEC with a single component, the new feature is the ad- 
ditional energy of the confining trap, and various regimes 
occur depending on the dimensionless interaction param- 
eter Na/d (see Sees. II.B and III.B.l). If Na/d Q < 1, 
then the vortex core is comparable with the condensate 
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radius (both of order do, the mean oscillator length) or 
the intervortex separation (see Sec. III.C). In the more 
typical Thomas- Fermi regime with Na/do ^> 1, however, 
the vortex core radius £ (the healing length) is much 
smaller than both c?o an d the still larger mean condensate 
radius Rq. 

For a three-component spinor condensate, the spatial 
variation of each component provides additional freedom 
in minimizing the overall free energy. Since the resulting 
coupled GP field equations are nonlinear, several solu- 
tions can exist, at least in principle, each corresponding 
to a local minimum of the GP energy [see Sec. 7.5.7 of 
IVollhardt and Wolfld (|l990T i for the similar situation in- 
volving the Ginzburg-Landau theory of superfluid 3 He]. 
In practic e, it turns out that |<?2| <C go f° r both 23 Na 
and 87 Rb (fYhl I1999T) . which means that gradient ener- 
gies and other interaction energies play a significant role 
in determining the condensate wave function. This near 
degeneracy is especially important in understanding the 
core of a quantized vortex in a spinor condensate, be- 
cause these various energies co mpete with the energy of 
the circulating flow. lYipl (|1999n proposes several different 
vortex structures, and many other authors have ampli- 
fied these considerations, both for single vortices [see, for 
example. dBulgakov and SadreevL 120031 : llsoshima et all . 
l200lllMizushima et all . 2002h and refer ences therein] and 
for v ortex lattices [see, for example, ( Mizushima et all . 
l2004f ) and references therein]. One experimental chal- 
lenge in studying a rotating spinor BEC is ensuring both 
the full axisymmetry of the optical trap and its alignment 
with the axis of rotation; otherwise the lifetime decreases 
significantly. These difficulties appear in the experiments 
on a rotating paired Fermi gas, where the magnetic field 
associated with a Feshbach res onance also requires a n 
optical trap [see Sec. IV. B. 6 and lZwierlein et al\ ( 20051 )]. 
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